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We review the progress and main challenges in implementing large-scale quantum computing by optical con- 
trol of electron spins in quantum dots (QDs). Relevant systems include self-assembled QDs of III-V or II- VI 
compound semiconductors (such as InGaAs and CdSe), monolayer fluctuation QDs in compound semiconduc- 
tor quantum wells, and impurity centers in solids such as P-donors in silicon and nitrogen- vacancy centers in 
diamond. The decoherence of the electron spin qubits is discussed and various schemes for countering the de- 
coherence problem are reviewed. We put forward designs of local nodes consisting of a few qubits which can 
be individually addressed and controlled. Remotely separated local nodes are connected by photonic structures 
(microcavities and waveguides) to form a large-scale distributed quantum system or a quantum network. The 
operation of the quantum network consists of optical control of a single electron spin, coupling of two spins in 
a local nodes, optically controlled quantum interfacing between stationary spin qubits in QDs and flying photon 
qubits in waveguides, rapid initialization of spin qubits, and qubit-specific single-shot non-demolition quantum 
measurement. The rapid qubit initialization may be realized by selectively enhancing certain entropy dump- 
ing channels via phonon or photon baths. The single- shot quantum measurement may be in- situ implemented 
through the integrated photonic network. The relevance of quantum non-demolition measurement to large-scale 
quantum computation is discussed. To illustrate the feasibility and demand, the resources are estimated for the 
benchmark problem of factorizing 15 with Shor's algorithm. 
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I. INTRODUCTION 



Based on the quantum parallelism rooted in the superposition principle of quantum mechanics, quantum computers are ex- 
pected to dramatically outperform their classical counterpart, particularly with exponential speedup in solving some hard prob- 
lems such as factoring 1 1 1. Toward the ambitious realization of practical quantum computation such as factoring a million-bit 
number, enormous eff'orts are still to be put on both the design and invention of software (quantum algorithms) and the hardware 
development (physical implementation). Here we are mostly interested in the latter part. In principle, all quantum systems 
(which arguably amount to all physical systems) could be considered for the physical realization. But certain qualifications 
(such as the DiVincenzo criteria |2|) are to be fulfilled for them to be brought into consideration. Still there is a vast range 
of systems in the candidate pool, including nuclear spins in liquids |3 -7|, trapped ions or atoms |8 - 18], atoms in optical lat- 
tices (T9H24I, photons I25H341 , superconducting circuits ||35H55]| , electrons suspended over liquid helium surfaces 156-591, 
molecular magnets f60l, nuclear spins in solids tSlllS^I, electron spins in semiconductor quantum dots (QDs) |63-79|, hole 
spins in QDs |80, 81], electron spins in impurity centers in semiconductors such as phosphorus donors in silicon [82- 86] and 
nitrogen- vacancy (NV) centers in diamonds (STUlOOL and non-Abelian any on excitations in quantum matters with topological 
orders 1 101;-. 105 J . Here we concentrate on solid-state systems, and in particular on electron spins in semiconductors, while 
nuclear spins therein are also considered either as an adverse noise source or as a beneficial information storage 11061 11071 . 
Electrical, magnetic and/or optical means may be employed to access and/or manipulate the spins. In this review, we discuss 
the optical operations which may be applied to electron spins in III-V compound semiconductors where the direct band-gaps 
facilitate controllable optical transitions. 

We shall not give an overall review of diff'erent schemes under current investigation for quantum computation. A compre- 
hensive review of progresses and challenges in research of diff'erent systems may be found, e.g., in the Quantum Computation 
Roadmap published by Quantum Institute, LANL. It, however, would be useful to have a perspective of the position and con- 
nections of the systems under review in the global picture of quantum computation research. As compared with their "soft" 
counterparts such as trapped atoms or ions, cavity-trapped or flying photons, electrons floating on liquid helium surfaces, and 
nuclear spins in liquids, the "hard" solid-state systems as candidates for quantum computers have the advantages of stability and 
integratability, but have the disadvantages of relatively short coherence time due to interactions with complex environments in 
solids. The solid-state systems under current investigation include superconducting circuits, nuclear spins in solids, and electron 
spins in semiconductors. While the qubits in superconducting circuits are made of excitations with macroscopic coherence in 
superconductors under designed confinement, nuclear or electron spins are natural qubit carriers since the information can be 
encoded in an intrinsic degree of freedom of elementary particles and thus are very stable. For example, while a superconduct- 
ing qubit may be lost during the measurement and control processes, a spin qubit would always exist unless the hosting particle 
like a nucleus or an electron disappears (by decay, ionization, thermal activation, etc.). Also, a spin does not feel an electrical 
field directly, which makes spins less vulnerable than superconducting qubits to charge or current noises from environments or 
operating devices. The decoherence of spin qubits may be caused by coupling to other environmental spins, local magnetic field 
fluctuations, or phonon scattering via spin-orbital interaction, all of which are usually rather weak in semiconductors. So the 
coherence time of spins are usually very long at low temperature and under a moderately strong magnetic field, varying from 
microseconds to milliseconds for electron spins (excluding the inhomogeneous broadening eff'ect) L108rflll 1 . and longer than 
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seconds for nuclear spins r il2H115l . The weak coupling of spins to environments, of course, has also detrimental effects - the 
control, initialization and measurement of spins are all challenging tasks. In this regard, electron spins are more tractable than 
nuclear spins. Nuclear spins have longer coherence memory time but slow operating rate and low detection efficiency. Electron 
spins are relatively more controllable but less resilient to decoherence. Schemes have been pursued to combine the advantages 
of the two kinds of systems by using nuclei as storage l l06lllQ7llll6i and electrons as operating units and interfaces 1 1 17 1. The 
coupling between nuclear and electron spins is the hyperfine interaction. The hyperfine interaction is a main mechanism causing 
electron spin decoherence, but it can also be utilized to realize coupled qubit systems. 

Electron spin qubits can be formed in various structures, such as doped electrons in QDs and impurities in solids. The fab- 
rication of such systems with designed patterns and structures is possible because of the advances in modem semiconductor 
technologies and nano-technologies - compatibility with the modern semiconductor industry is an extra advantage of using elec- 
tron spins as qubits in quantum computation. The direct pump, control and probe of electron spins may be done by electron spin 
resonance techniques with microwave pulses |88, iTOlllllllllTl . For faster operation clock as desired, the electron spin states 
may be converted to other degrees of freedom through quantum interfacing (the same as one does from nuclear spin states to 
electron spin states). Then electron spins may be accessed indirectly by control of agents such as excitons in semiconductors and 
photons generated by recombination of excitons [64l [6814791 II 1 8 -1201. Quantum interfacing between spins and photons |119l 
makes also possible quantum communication between distributed quantum nodes which is required for scalable quantum com- 
putation. In this sense, the direct-gap semiconductors such as In As and GaAs, where ultrafast optical control and interfacing 
are possible, have some advantages over the indirect gap materials such as silicon, where electrical gating and microwave pulses 
may be the only possible means of control. As compared with silicon, a main concern with the III-V materials is the much 
shorter spin coherence time due to the abundance of nuclear spins as noisy environments (in GaAs, e.g., the electron spin coher- 
ence time is in the order of microseconds, while in silicon, it is in milliseconds, excluding the effects of phonon scattering and 
inhomogeneous broadening) 11121141321 . Fortunately, there already exist various schemes to elongate the spin coherence time by 
orders of magnitude, via dynamical control 1 1301 II 33^1421 . or nuclear state preparation lll3H115lll43H151l . 

All these said, we would like to remark that at this point, it would be premature to discourage effort in exploring different 
physical systems, existing or emerging. It is conceivable that the future quantum computers will be realized by combination of 
innovative technologies, ideas, concepts, and synthesis of materials and systems. For instance, the idea of topological quantum 
computation may be implemented with trapped atoms in optical lattice 1 152|; the systems under the focus of this review involve 
both stationary electron spins and flying photons I119L which may also be applied to coupled systems of photons and supercon- 
ducting qubits 147114914511 : semiconductor chips may provide micro-trap for ions; photon-based quantum computation may use 
quantum lights from trapped ions 1 153 1, atoms 11541 [T55 K or QDs [ 156 , 157|; and so on. In the present initial stage of quantum 
computation technology, it would be highly risky to exclude certain candidates just because of difficulties encountered in the be- 
ginning of the adventure, since different systems may have their bottleneck problems at different stages. In particular, solid-state 
systems, while promising with their large-scale stability and integratability in the future, are still facing severe obstacles of envi- 
ronmental noise and control errors for one or a few qubits. In this review, based on many experimentally demonstrated elements 
and theoretically proposed schemes, we would like to put forward blue prints of relatively large-scale quantum computing via 
optical control of electron spins in QDs. We should point out that such targets are by no means easy and still require significant 
advances of technologies and concepts. Also, although our discussions, to be specific, will be based on electron spins in InAs or 
GaAs QDs, the schemes, with certain modifications, can be applied to a few emerging novel systems such as hole spins in QDs 
and NV centers in diamond, where the physics is similar to electron spins in QDs. 

There are two main concerns with solid-state systems for quantum computation, when compared with their atomic peers like 
trapped ions or atoms. One is how to fabricate and construct a large array of reproducible or identical qubit units (such as QDs). 
And the other one is the many-body problem. 

The fabrication issue is even worse for the systems to be discussed in this review, namely, QDs under optical control. Most 
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FIG. 1. Relevant energy scales for a typical GaAs semiconductor QD, in terms of frequency, energy, and temperature. The hyperfine constant 
refers to the interaction between a single nucleus and an electron. 
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likely, the QDs are formed by molecular beam epitaxy (MBE) growth process such as the self-assembled QDs and the fluctuation 
QDs. Efl'ort is being made on growth process so well controlled that all the QDs are almost identical and regularly located. 
There are also such promising systems as impurities in semiconductors such as phosphorus donors in silicon 1 83 1 and the NV 
centers in diamond 115811159 1 where the nanometer-precision ion implantation technology may allow a patterned array of qubits 
represented by the electron spins in the impurity centers. In the foreseeable future, however, we may have to live with the problem 
of the irregularity. On the positive side, we could take advantage of the irregularity produced by the system fabrication. In small 
systems for demonstration purposes, the varying size and position of QDs may be used as fingerprints by which diff'erent qubits 
may be addressed and individually controlled by light beams which usually have resolution no better than half the wavelength. 

The many-body problem raises two related questions: Can an isolated qubit be properly defined at all among so many particles, 
and will a coherent superposition state of a qubit last long enough for quantum gates before it collapses due to decoherence in 
the noisy environment? 

It turns out 1 160] that large gaps between valence bands and conduction bands in semiconductors protect elementary excita- 
tions such as excitons and electrons so well that an extra electron doped in a QD is well-defined as a single particle moving with 
a renormalized eff'ective mass and coupling constants such as the g-factor, in analogy to a low-energy electron in the Dirac sea. 
The protective gap sets a fundamental limit on the operation speed of quantum computation in solid-systems. Such a limit is 
far from being approached in current experiments, whether in typical semiconductors where the gap is about 1 eV (for fs-order 
operation time) or a few meV in superconductors. 

In optical control of spins in QDs, decoherence is caused mainly by three mechanisms, namely, coupling to nuclear spins of the 
host lattice 1 108, 121-132, 161], phonon scattering [162 - 168 1, and spontaneous photon emission during the optical control |68j 
11691 . For some pseudo-spin qubits with orbital-state dependence (such as the singlet-triplet qubits I1Q8I ), the charge fluctuation 
may also contribute to the environment noise 11701 [TTl I. Besides the standard quantum error correction protocols 1 172], there 
are specific strategies to deal with various decoherence mechanisms, which are needed anyway to achieve fidelity of quantum 
gates above the threshold for quantum error correction. The spontaneous photon emission could be suppressed by completing the 
optical control rapidly or via ofl'-resonance (virtual) excitation L68l[72l[76lll73i The phonon scattering may be quenched simply 
by lowering the temperature to a few kelvins I162H1651 117414177II or by using light-element materials such as diamond 1871 - 
[100] or organic materials | f78l 11791 where the spin-orbital coupling is weak. The nuclear spins, being a slow bath, may have 
their decoherence controlled by certain dynamical decoupling or disentanglement control [ 130II133I4142L Again, the normally 
harmful noise sources could be made useful by design. The photon and phonon baths are rapid entropy dumping pools when 
certain quantum channels are selectively enhanced |120|. The photon emission, when enhanced by cavities and guided by 
quantum channels, is an important basis for quantum communication between remotely separated qubits | 76, 119|. The nuclear 
spins, having very slow dynamics, are considered as good local quantum memories 11061 11071 11161 with an electron spin in 
contact acting as a mediator for quantum information operation and transfer. 

To achieve large-scale quantum computation, consensus has been reached on several criteria to be fulfilled, known as the 
"DiVincenzo criteria" |2|. We quote these criteria below as the guidelines for reviewing the progresses and main challenges 
toward the realization of quantum computation by optical control of electron spins in QDs: 

7. A scalable physical system with well characterized qubits; 

2. The ability to initialize the state of the qubits to a simple fiducial state; 

3. Long relevant decoherence times, much longer than the gate operation time; 

4. A universal set of quantum gates; 

5. A qubit-specific measurement capability; 

6. The ability to interconvert stationary and flying qubits; 

7. The ability faithfully to transmit flying qubits between specified locations. 

The stationary qubits under our focus are well defined by electron spins in QDs, and the flying qubits carrying quantum 
information between distributed nodes are photons flying in waveguides. We will discuss in more details the decoherence of 
the spin qubits and show that the decoherence time (~ 10"^ sec in a typical GaAs QD with the inhomogeneous broadening ex- 
cluded CMIIOSIIIIIIIlTlIIlOl^^ is indeed much longer than the quantum gate operation time (~ lO'^^ sec 11791 fTsTlfT^ ). 
The one- and two-qubit gates, which form a universal set 1 17T,'183","1841, are realized by optical excitation of charged excitons. 
Some fundamental physics issues with the initialization and quantum measurement of qubits will be reviewed. Measurement 
and initialization are put together because they are related to the same physical process. Initialization disposes of entropy to the 
environment while in measurement the environment acts as part of a readout device. For a large-scale quantum computation 
blueprint, we use designs of local nodes of a few qubits and structures of distributed nodes connected by quantum channels 
which may be realized by photonic elements such as waveguides and microcavities. Control schemes of quantum interfacing 
are also an important topic to be covered. To illustrate the feasibility and demands of the quantum computation in the discussed 
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systems, the resources, in terms of the number of optical pulses and operation time (compared with the spin decoherence time), 
will be estimated for the benchmark problem of factoring 15 with Shor's algorithm. 



II. SPIN QUBITS IN QUANTUM DOTS 



In this section, we begin (in Sec. |IIA|) with a brief review of the confinement of single electrons in optically controllable 
semiconductor QDs, followed (in Sec.|IIB|) by discussions of QD energy level structures and optical properties. In Sec.|IIC[ we 



briefly outline recent theoretical and experimental results on the spin coherence properties of single electrons confined in QDs. 
Both theories and experiments show that, as phonon mechanisms are suppressed at low temperature (~ 4 K and below), lattice 
nuclear spins become the dominant cause for the electron spin decoherence. In Sec. II D we review the theory of electron spin 
decoherence by interacting nuclear spins in a QD. Coherence protection of electron spin in the interacting nuclear spin bath is 
possible by applying a sequence of n pulses to the electron, as discussed in Sec. HE In Sec. IIP an overview of QD electron 
spins as qubits is given from the perspective of fault-tolerance requirement for scalable quantum computation, and two other 
promising spin qubit systems, namely, hole spins and NV center spins, are also discussed. 



A. Confinement of a single electron in a quantum dot 



Two types of MBE grown QDs formed in direct bandgap III-V compounds ofl'er a great deal of controllability by ultra-fast 
optics and are being investigated as building blocks for optically manipulated quantum computers. 

The first is referred to as interface fluctuation formed QDs in GaAs/AlGaAs quantum well structures [185-1931. We will 
refer them in short as GaAs fluctuation QDs. As illustrated in Fig.|2ja) and [3] an electron in such a structure is confined in 
the growth direction along the z axis in a low bandgap GaAs layer, with a thickness of tens of A, between two higher bandgap 
PAxGdii-xA.^ layers. In the III-V materials, the conduction band minimum occurs at the F point in the momentum space and 
the hetero structure wavefunction of the electron is constructed from the conduction band Bloch functions in the vicinity of the 
r point. Thus the diff'erence in energy of the conduction band minimum in the high and low bandgap materials form a square 
well potential for electron in the growth direction. The band discontinuity in the interface of the low and high bandgap material 
is typically hundreds of meV so that the vertical confinement is strong. Growth interruption leads to roughness at the material 
interfaces, usually a thickness fluctuation of one monolayer. Electron confinement within the plane of the quantum well (in 
the X and y directions) is caused by the quantum well thickness fluctuation. The energy lift by the quantum confinement in the 
z direction is roughly ^(f )^ for the lowest energy state, where m* is the eff'ective mass of conduction electron and d is the 
quantum well thickness. This energy lift is larger where the quantum well is thinner. Therefore lateral confinement is formed 
where the quantum well has an island. Fig. [3] shows schematically how monolayer- size fluctuation in a quantum well gives rise 
to the localized envelope function in the plane. The energy scale of this lateral confinement is typically of several to tens of meV 
in GaAs fluctuation dot. A GaAs fluctuation dot with lateral size ~ 40 nm can hold several localized energy levels with level 
spacing of several meV. 

The second type of QDs, referred as InAs self-assembled dots, are formed using the Stranski-Krastanow growth mode which 
utilizes the strain caused by the lattice mismatch between InAs layers and GaAs substrates. InAs self-assembles into islands 
which are primarily in the shape of a pyramid [see Fig. [2](b)], with a height of tens of Aand a base size of tens of nanome- 
ters 1 196 - 199 1 . In an InAs self-assembled QD, the lateral confinement is much stronger due to the pyramid structure, and the 
smaller dot size leads to level spacing of ~ 10 meV or larger. 




FIG. 2. Illustration of two types of optically controllable QDs. The arrows indicate the growth direction. The ellipse regions show schematically 
the confinement of electrons in the dots, (a) A GaAs fluctuation QD. (b) An InAs self-assembled QD. 
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For a GaAs fluctuation QD, a single conduction band electron can be incorporated in the dot by modulation Si doping in 
the barriers |200|. For an InAs self-assembled dot, gate voltage tuning in an n^-intrinsic-Schottky diode structure is a more 
controllable way to charge and discharge the QD with a single electron |194[|195l (see the schematic illustration in Fig.|4]). The 
qubit is typically encoded in the spin subspace of the lowest energy level of the single electron in the QD. 

There are other notable systems where single electrons are localized in nanoscale regions in semiconductors. These include the 
confinement by electric gates on top of two-dimensional electron gas in GaAs (63] 201 , 202j, and the localization by impurities 
such as phosphorus donors in silicon | 82- 86| or NV centers in diamond 1881 [92H98nl00inT7ll203l [2041 . The spins of single 
electrons localized in these systems are also under extensive investigation as qubit carriers. 



B. Energy levels in a charged quantum dot 

Control of a spin qubit makes use of a larger Hilbert space in a QD, involving optical transitions from valance bands to 
conduction bands. In this subsection, we briefly describe the relevant energy level structures and the corresponding optical 
transition selection rules. 

A direct interband transition creates an additional electron in the conduction band by leaving a hole in the valance band. At the 
r point, the fourfold degeneracy of the bulk Fg valance band is lifted by the quantum well confinement eff'ect. The top valance 
subband is a doublet derived from the = ±3/2 bulk band, which is also denoted as the "heavy hole" band (see Appendix. [A|. 
Other valance bands are irrelevant in our control schemes as they are always far off'-resonance. Electrons can be excited from the 
valance band states of angular momentum to conduction band states of spin S^hy absorbing a photon, with the selection rule 
= Jz + cr where cr = ±1 corresponds to the circular polarization of light. Therefore, the single electron states are optically 
coupled to the charged exciton states (also known as trion states) composed of two conduction electrons and one heavy hole. 
Fig. [5] shows schematically two diff'erent types of trion states. The two electrons in Fig.[5ja) occupy diff'erent electronic levels 
of the QD and the two electrons in Fig.[5jb) are in a spin singlet configuration on the same electronic level. The energy spacing 
between diff'erent trion configurations is ~ 1-10 meV depending on the QD size. The bandwidth and Rabi frequency of the 
optical field is much smaller than this energy spacing. So when the frequency of the optical field is near the resonance of one 
type of trion configuration, the remaining trion configurations can be neglected. 

Most manipulation schemes use the ground state trions [Fig. |5jb)] to mediate optical control of the spin. In such cases, 
the relevant Hilbert space is composed of the two single electron spin states and the two ground state trions (see Fig. [6]). For 
simplicity, we use below the e\ to denote the creation of a conduction band electron of 5"^ = ±1/2 in the lowest energy level 
of the QD and similarly h\ to denote the creation of a heavy hole with = ±3/2 (annihilation of an electron with = =f3/2). 



(a) (d) 

AlGaAs GaAs AlGaAs V^{x) V^{x) 




■> z > z 



FIG. 3. Schematics of the three-dimensional confinement of electrons and holes in a GaAs fluctuation QD. (a) Interface fluctuation, typically 
of one monolayer, forms a fluctuation dot. The bright region is the GaAs material and the dark regions are AlGaAs. (b) Vertical confinement 
for electrons and holes in the growth direction (the z direction) at x = Xi where the quantum well is thinner d - d\. (c) Vertical confinement 
in the growth direction dX x - X2 where the quantum well is one monolayer thicker d - d2. (d) The difference of the energy lift by vertical 
confinement in regions of different thickness forms a lateral confinement for the electrons and holes. 
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The relevant Hilbert space for optical control of spin is: | t> = e\\G), \ i) = el\G), Its) = elelhl\G} and \t_3} = e\elhl\G) 
where |G) denotes the configuration with empty conduction bands and full valence bands. The transition selection rule between 
these four states is shown in Fig.|6ja). For diff'erent control schemes, the selection rules represented in other basis sets are also 
useful. In Fig. [6|b) and (c), we have changed the basis for the spin states and the trion states so that the electron spin states 
are eigenstates along x-direction and the trion states are |r±) = (1^3/2) ± k-3/2)) / V2. The selection rule can also be represented 
with the linearly polarized basis for the optical field as shown in Fig.[6jc) 11181 . The selective coupling of the QD transitions 
with light of diff'erent polarizations off'ers sufficient freedom for optical control of the spin. Depending on the polarization and 
frequency of the optical field, various optical pathways can be established to realize the control in the spin subspace via a second 
order process through the trion. For example, if the optical field is (T+ circularly polarized [Fig.|6jb)], the \t_3) trion state is 
decoupled from optical field and the relevant dynamics is in a A-type three-level system. Raman processes in such three-level 
systems are central to the optical control of spin dynamics as will be discussed in details later. 



C. Spin relaxation and decoherence in quantum dots 



Coherence properties of qubits are crucial to quantum information processing. Spin decoherence of a single electron in a solid 
results from the coupling to various environmental modes. Typical optical manipulations of spin qubits in QDs are performed in 
the Voigt or Faraday geometry with a strong magnetic field (-1-10 T), and at low temperature (~ or < K) to suppress thermal 
excitations in the environment 1205 [ I206i] . In this subsection, we single out the environmental eff'ects that dominate the spin 
decoherence under these experimental conditions. 

We first show the quantum mechanics of decoherence due to the coupling to a general environment I207H2091 . The initial 
state of the electron spin, |^^(0)) = C+\+) + C_|-), is prepared as a coherent superposition of the spin up and down states 
1+) in an external magnetic field. The state of the total system of the spin plus bath at that instant forms an unentangled 
state, 1^(0)) = l^'(O)) \J). It evolves over time t to |^(0) = C+(OI+) ® \J^{t)) + C-{t)\-) ® \J-{t)) where the bath states 
\J'^{t)) and \J'~{t)) are generally diff'erent. The mixed state of the electron spin is determined by its reduced density matrix 
obtained by tracing over the environment states p^^/(0 = C^/(OCo-(0(^^ (01^^(0)- The diagonal element of the reduced 
density matrix p^^, gives the probability of finding the spin in state |(t). Either off'-diagonal element is a measure of the 
spin phase coherence. The environment-driven transfer of the probability between the spin states is known as the longitudinal 
relaxation, while the loss of the off'-diagonal element is known as the transverse decoherence. The longitudinal relaxation 
also contributes to the transverse decoherence. The spin decoherence without longitudinal relaxation is called pure dephasing, 
characterized by the quantity {J'^ {t)\J'~ {t)) . Pure dephasing is thus a consequence of system-bath entanglement when the 
bath evolution \ J'-(t)) is conditioned on the system states |±). The above discussion can be generalized to the situation where 
the environment is initially in a mixed state. The single electron plus its environment is then described at the initial time by 
the density matrix p(0) = \(p\Qi)){(p\0)\ Y^j Pj\J'){J'\, where the reduced density matrix for the electron at any moment is 

p;,,.(0 = Y.jPjC':^Xt)CAt){J''\t)\J''{t)). 

Novel experimental techniques have been developed in the past several years for measurement of the spin relaxation and de- 
coherence in QDs. In GaAs fluctuation dots and InAs self-assembled QDs, optical techniques are commonly used for measuring 
the spin Ti time. These include the time domain measurement by optical generation and detection of non-equilibrium spin pop- 
ulation L176.L and the frequency domain approach of the coherent phase-modulation spectroscopy 11 77 J . In gate-defined QDs, 
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FIG. 4. Gate controlled charging of an InAs self- assembled QD in an n+-intrinsic-Schottky (NIS) diode structure, (a) Schematic illustration of 
the NIS diode structure given in II194II195I . (b) Band diagram of the NIS diode structure where the InAs QD is uncharged, (c) An InAs QD 
charged with a single electron. 
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spin-to-charge conversion can be implemented to trace the time evolution of the initially created spin population 1174111751 . In 
strong magnetic field of ~ 1 - 10 T and at low temperature (< 1 K), the experimentally measured value Ti ~ 10""^ - 10"^ s II 1741 - 
11771 are in good agreement with the theoretical predictions of the spin longitudinal relaxation induced by phonon ll62H165ll . In 
a recent experiment by Amasha et al | 210 |, spin relaxation rate of a single electron in a lateral QD is studied when the orbital 
wavefunction is manipulated using gate voltages. The measured dependence of of Ti on orbital confinement and magnetic field 
is in excellent agreement with theory ll63M166t which confirms that phonon scattering in the presence of spin-orbit coupling 
is the dominant cause for spin relaxation in magnetic field down to 1 T. In low magnetic field and when spin-orbit coupling 
strength is weak, Ti as long as 1 sec is observed | 210 |. 

Transverse decoherence may be measured in optically accessible QDs by the frequency domain approach utilizing the Hanle 
eflfect 1 112], or by the time domain pump-probe measurement. In the later approach, a circularly polarized pump pulse initiates 
spin polarization in the growth direction, whose precession about an in-plane magnetic field is tracked by the diff'erential trans- 
mission of a circularly polarized probe beam |205 | or the Faraday rotation angle of a linearly polarized probe beam Il211i . In a 
gate-defined double-dot structure, the spin-to-charge conversion process can also be implemented to probe the relative coherence 
between states of a coupled spin pair, which provides information about the decoherence of a single spin I2Q2II212|| . In these 
experiments, transverse decoherence times of ~ 1 - 10 ns were obtained either from measurements on spatial ensembles 
of QDs Ii205l 12111 12131 or from time-ensemble measurements of single dots II1121 12021 [2121 . Spin echo type of measurement 
was also performed on the gate-defined single GaAs dot IIO8III8OII , which showed an echo decay time of T// ~ yus. The sharp 
diff'erence between Th and T* suggests that the ensemble dephasing is mainly aff'ected by the inhomogeneous broadening of the 
local environment of the QDs, which is removed in the spin echo measurement. Greilich et al have shown that by using a peri- 
odic train of circularly polarized light pulses to excite an ensemble of InAs self-assembled QDs, spin polarization is amplified 
only in a subset of the ensemble where the electron spin Zeeman frequency has to be an integer multiple of the pulse repetition 
frequency 11091 . These quasi-discrete spectra leads to constructive interference of the spin precession at each pulse arrival time. 
As dephasing by inhomogeneous broadening is thus removed, single spin T2 ~ //s is obtained 1 109, ll24[|146[|150lll51l . Very 
recently, Hahn echo measurement was also made possible for impurity spins in GaAs where rotations of spins were achieved by 
ultrafast optical pulses I181II . The measured Hahn echo decay time is 6 yus, consistent with theories and other experiments. 

These experiments all show that transverse decoherence times (r*, Th and T2) are orders of magnitude faster than the lon- 
gitudinal relaxation time Ti . On the other hand, theoretical analysis of the phonon mechanisms concludes that pure dephasing 
due to phonon is well suppressed at the temperature where these experiments are performed (< 1 K) I166M167I . As phonon is 
unlikely to be responsible for the observed fast transverse decoherence, the remaining possibility is then the nuclear spins sitting 
on the lattice sites which are coupled to the electrons through the hyperfine interaction 1. 1 08 . . 1 21. - . 1 27 . . 1 29,-. 1 32. . 1 34. . 1 6 1 . .2 1 4 J . 



D. Decoherence by an interacting nuclear spin bath 

In this subsection, we discuss the eff'ects of the lattice nuclear spins on the electron spin coherence. In the relevant III-V 
materials, all stable isotopes have non-zero nuclear spins. The nuclear magneton is about 3 orders of magnitude smaller than 
the electron Bohr magneton. A strong magnetic field of 10 T only results in a nuclear zeeman energy of ~ mK, much smaller 
than the experimentally achievable temperatures in cryostats (~ 4 K) or even in dilution refrigerators (> 50 mK) (see Fig.[T]). 
Therefore, coupling to a thermalized nuclear spin bath will be an inevitable source of decoherence for quantum computation in 
III-V materials. As the longitudinal spin relaxation of the electron is found to be much slower than the transverse decoherence, 
we will focus on the pure dephasing of electron spins. In this Subsection, we set the z-direction along the direction of the external 
magnetic field. 



(a) (b) 

~y i^r ~y ^ 

vAAAAAAA/^ vAAAAAAAA^ 



FIG. 5. Schematic illustration of the creation of trion states in a QD. The QD initially holds a single conduction band electron with spin up 
in the lowest state of the confinement potential, (a) By the Pauli exclusion principle, a cr- circularly polarized light can create an additional 
exciton to a higher excited state of the QD. (b) cr+ polarized light can create the lowest energy trion state with two electrons forming a singlet 
and a hole in the spin up state. 
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1. Single electron in a mesoscopic bath of interacting nuclear spins 



In III-V semiconductors, the electron spin is coupled to the lattice nuclear spin through the contact hyperfine interaction. The 
averaged magnitude of the coupHng to a nucleus is inversely proportional to the total number N of nuclei in the QD. For QDs of 
all practical sizes (A/^ ~ 10^ - 10^), this hyperfine coupling is much stronger (~ MHz) than the mutual interactions between the 
nuclear spins (< kHz). Therefore, a mesoscopic bath consisting of all nuclear spins within the QD (i.e. in direct contact with the 
single electron) can be identified (see Fig. [7]) 1 130|. The coupling between the mesoscopic bath and nuclei outside the boundary 
can be neglected since such dynamics occurs in a much slower timescale as compared with the electron spin decoherence caused 
by the mesoscopic bath. The decoherence problem can then be solved by considering the quantum dynamics of the coupled 
mesoscopic system of electron and nuclear spins. The assumption that this mesoscopic system is well isolated from the back 
ground lattice has been confirmed in numerical studies where the boundary of the mesoscopic bath is systematically extended 
and the electron spin coherence shows fast convergence 1 130|. The nuclear spin bath is typically of a randomized configuration 
as schematically illustrated in Fig. [7] since the experimentally achievable temperature is always much higher than the nuclear 
Zeeman energy. 

We briefly describe below the key ingredients for electron spin decoherence in the high field limit with details given in 
Appendix |B] By the diagonal part of the electron-nuclear hyperfine interaction, the electron Zeeman energy is conditioned on 
the nuclear spin states (see FigjSja)). In a QD ensemble, the diff'erent nuclear spin configurations for each ensemble member 
lead to inhomogeneous broadening of the Overhauser field for the electron spin. This is the dominant cause of ensemble 
dephasing in the timescale of which is inversely proportional to the inhomogeneous broadening. Nuclear-nuclear interactions 
become relevant due to the non-uniform hyperfine coupling strength between the electron and diff'erent nuclear spins. Pair- wise 
nuclear flip-flops can then lead to dynamical fluctuation of the nuclear Overhauser field (see FigjSjb)). This is the cause of 
single electron spin decoherence in the nuclear environment which begins with a pure state (referred as single-system dynamics 
hereafter). Electron-nuclear coupling also has an off'-diagonal part which tends to cause flip-flop between the electron and a 
nuclear spin. Because the electron Zeeman energy is much larger than the strength of the hyperfine interaction, the real process 
of electron nuclear flip-flop is suppressed fl24 127 ] . However, a second order process which consists of two virtual flips of the 
single electron ends up as a flip-flop between two nuclear spins (see FigjSjc)). This efl'ective nuclear interaction due to the single 
electron is designated as the extrinsic interaction II1271I1301|134L as opposed to the intrinsic nuclear interactions that exist in the 
semiconductor matrix, e.g. the dipole-dipole coupling and the indirect coupling mediated by virtual interband transitions via the 
hyperfine interaction |[215fl219|. 

The extrinsic nuclear interaction couples any two nuclear spins within the mesoscopic bath and is therefore infinitely-ranged, 
i.e., throughout the entire mesoscopic region. By contrast, the intrinsic nuclear interaction is finite-ranged. For near neighbors, 
the intrinsic one is much stronger than the extrinsic one for the field strength under consideration. In addition, for the extrinsic 
nuclear interaction, the magnitude is inversely proportional to the external magnetic field and the sign is conditioned on the 
electron spin states. 



(b) (c) 

1^3/2) 1^-3/2) |r-) |r+) 

K ixi 

FIG. 6. Optical transition selection rules in a QD illustrated in various basis sets, (a) The basis is the eigenstates of J^. The solid two-headed 
arrow denotes cr+ polarized light and the hollow two-headed one cr_ polarized light, (b) The two electron spin states are transformed to the 
basis in the x direction: \x±) = (| t> ± I V2. The solid two-headed arrow denotes cr+ polarized light, (c) The two trion states are also 
transformed to the basis in the x direction: \T±) = (|?3 > + lk_3»/ V2. Here, the hollow two-headed arrow denotes X-polarized light and the 
solid one 7-polarized light. 
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2. Nuclear spin pair-flip excitations and pair correlation approximation 



From the quantum mechanical picture of decoherence (see section [nC]), pure dephasing of a single quantum system is caused 
by the bifurcation of environmental evolution under the drive by different system states, or system-bath entanglement. Thus, the 
nuclear bath evolutions conditioned on different spin states of the electron are key to the solution of electron spin decoherence. 
The elementary excitations in the nuclear spin bath are pair-flip excitations as shown in Fig. [Sjb) and Fig. [Sjc). The flip-pairs 
are independent of each other if they are well separated. For a typical QD with A/^ ~ 10^ nuclear spins, the number of pairs 
that can flip-flop in a random configuration is large. We have 0{N) local flip-pairs whose dynamics are dominated by the finite 
range intrinsic nuclear interactions, and 0(N^) non-local flip-pairs whose dynamics are dominated by the infinite range extrinsic 
nuclear interactions. On the other hand, the number of pair-flips that can occur on the timescale of electron spin decoherence is 
negligibly small as compared with N. This is due to the slowness in the nuclear spin interacting dynamics. Thus, the probability 
of two pair-flips occurring in the neighborhood of each other is negligibly small and pair-flips as elementary excitations can be 
treated as independent of each other (see Appendix [B] for details). This approximation is further confirmed by the linked cluster 
expansion approach |129|. Independent pair-excitation approximation corresponds to keeping the lowest (2nd) order linked 
diagrams in the exponential factor. Higher order linked diagrams contain more nuclear interaction lines and are negligible in 
the relevant timescale because of the weakness of the nuclear interactions as compared with the electron-nuclear coupling. For 
long time evolution (such as for decoherence under pulse control) and for relatively small spin baths, higher order correlations 
would be important. The linked cluster expansion would be increasingly inefficient for calculating higher order correlations. 
The density matrix cluster expansion is an alternative method which requires no evaluation of higher order Feynman diagrams 



and is thus quite convenient I22QI . For small spin baths, however, it has been shown that the cluster expansion may not converge 
to the exact results |221 1. For a systematic and accurate account of the higher order correlations in spin baths in the qubit 
decoherence problem, a cluster-correlation expansion method has been developed 1221. 222 L which covers the valid ranges 
of all the methods mentioned above and in particularly produces the exact results even for a relatively small bath where the 
standard cluster expansion fails. The cluster correlation expansion is based on the factorization of the bath dynamics into non- 
factorizable correlations of certain groups of bath spins. The lowest order of the cluster correlation expansion coincides with 
the independent pair excitation approximation. We also note that recently the higher order effects of the hyperfine interaction 
(beyond the pair-excitation considered here) has also been considered for a relatively weak external magnetic field while the 
pure dephasing condition is still satisfied I I132I I214II 

The evolution of independent pair-correlations can be described using a geometric representation. A pair-flip k can be mapped 
to a Bloch vector which precesses about a pseudo-field (see Fig.|9]) 

ht^(±2Ak + 2Bk.0,±Ekl (1) 
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FIG. 7. Schematics of an electron (the shadow) and one layer of lattice nuclear spins in a QD. The two boxes in dotted lines indicate two 
possible choices of boundary of the mesoscopic nuclear spin bath, which are relatively arbitrary due to the interaction between nuclei within 
and without the boundary. When the hyperfine interaction dominates over the nuclear spin interaction, such arbitrariness has negligible effects 
on calculation of the electron spin decoherence as long as all the nuclei in direct contact with the electron spin have been enclosed. 
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where, for the electron spin state |±), ±Ak and are the pair-flip transition amplitudes contributed by the extrinsic nuclear inter- 
action and the intrinsic nuclear interaction, respectively, and ±Ek is the energy cost of the pair-flip contributed by the hyperfine 
interaction (see Appendix [B]). At the initial time when the electron spin coherence is prepared (see the general formulation of 



decoherence process in Section II C), the Bloch vector for each pair-flip points in the pseudo -\-z direction. As the pseudo-field 
direction depends on the electron spin state |±), the pair evolution takes diff'erent trajectories conditioned on the electron spin 
state, and the distance 6k (i.e., distinguishability) between the two separated trajectories turns out to be a measure of the electron 
spin decoherence, 

k 

where X+ _(0 is the ratio between the electron spin coherence at time t and at time (see the general formulation of decoherence 
in section II Q . 6k also quantifies the amount of entanglement between the electron spin and the ^th flip-pair. The electron spin 
decoherence is thus the consequence of the entanglement with the pair-flip excitations in the interacting nuclear spin bath. 

With single system dynamics solved for an arbitrary random configuration, ensemble dynamics is simply the statistical average 
of the single system dynamics with the nuclear bath initially in all possible configurations. As the number of flip-pairs is large 
(0(N) for local pairs and O(N^) non-local pairs), the central limit theorem of statistics leads to a factorized form for ensemble 
spin coherence, 

i;+,-(0 = x^,-(0x4'-(0, (3) 

where X+ _(0 is the single-system decoherence in a typical configuration of the nuclear bath, and 

£fl(t) = J]Pje-'^'^'^^ (4) 
J 

where (pj(t) = (Q^ + 8j)t in free-induction decay, and the summation runs over all possible nuclear configurations J^. fie 
and 8 J are the electron Zeeman energy resulting from the external magnetic field and from the Overhauser field, respectively, 
with the latter dependent on the nuclear configuration J'. The ensemble eff'ect resides entirely in the factor x!+^_(t), which may 
be read as the inhomogeneous broadening of the Overhauser field 8j with a distribution function Pj. The inhomogeneous 
broadening eff'ect dominates the free-induction decay (FID) in the ensemble dynamics in the form of X!+^_(t) = ^-«^e?-(?/7^2)\ with 
the dephasing time T* ~ ^/NJ^-^ ~ 10 ns as measured ifTTIl [ml fT24l l20l l2Q5l [2T^ where Ji ~THz is the hyperfine 
constant of the material. 

A sequence of tt pulses can be applied to the electron spin to eliminate the effects of the inhomogeneous broadening I iTOl 
11321 1180[ I224II . In a general scenario where the electron spin is flipped at time ti, T2, . . . , and t„ respectively, we have 



(a) 



e-n interaction (diagonal) 



(b) 



n-n interaction (intrinsic) 
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FIG. 8. Nuclear spin processes relevant for electron spin decoherence. (a) By the diagonal part of electron nuclear hyperfine interaction 
(which involves only the spin vector components along the field z direction), the electron Zeeman energy depends on the nuclear configuration, 
(b) Nuclear pair- wise flip-flop by the intrinsic nuclear interactions, (c) Nuclear pair- wise flip-flop mediated by two virtual flips of the electron 
spin, which results in an eff'ective extrinsic nuclear interaction. is the electron Zeeman energy in the external magnetic field, and a„ is the 
hyperfine interaction strength between the electron and the ;ith nuclear spin. 
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0^(0 = (Q, + 6j) [Ti - (t2 - Ti) + • • • + (-ir(r - Tn)l When n - (t2 - rO + • • • + (-ir(r - r J = is satisfied, £^^^0) = 1 
and a spin echo is expected. The echo magnitude will be determined by the dynamical part X+ Under the simplest scenario, 
a single 7r-pulse is applied at r and a spin echo is expected at ^ = 2t, known as Hahn echo 12251 . The spin echo profile, i.e., the 
echo magnitude X+ _(2t) as a function of the echo delay time 2t, reveals the dynamical processes that leads to decoherence. 

It is worth noting that the factorized form of the ensemble spin coherence, Eq. ([3]), allows direct observation of single-system 
dynamics behavior £+ _(t) from a spatial ensemble measurement when dephasing by inhomogeneous broadening is removed at 
a general time t. For example, the mode locking experiment reported in 11091 1211 1 opens up such possibilities as discussed in 
section II C where the single spin T2 (defined here as the FID timescale of has been extracted from the experimental 

data. 



3. Timescales of single spin decoherence and ensemble spin echo decay 

In FID, the conjugate Bloch vectors precess along opposite directions for non-local pairs {k e Ka), and symmetrically with 
respect to the pseudo y-z plane for the near-neighbor pairs (k e Kb) [Fig.|9j. The decoherence can be readily grouped by the two 
diff'erent mechanisms as 

keKB keKA 

where = |h*|. We can see that the extrinsic hyperfine-mediated and the intrinsic couplings lead to the ^"(^/^24)^ and the 
^-{tiT2,Bf behavior respectively in time shorter than the inverse pair-flip energy cost (which corresponds to the width of the 
excitation spectrum), 

T2,B ~ b-^'^J{-^'^N^'^\ T2,A ~ ^e^r^N, (6) 

where b is the typical value of near neighbor intrinsic nuclear coupling strength Bk (see Appendix [B]). The super-exponential 
decay behaviour of the spin coherence indicates the strong non-Markovian characteristic of the bath dynamics in the short time 
limit. In the long-time limit, the super-exponential decay will change to an exponential decay time first which indicates the 
onset of the Markovian dynamics 11301 11311 1214II . The dynamics in the even longer time limit (which could occur, e.g., in a 
highly polarized spin bath), determined by the complex structure of the collective modes of the bath, is rather complicated, and 
power-law decays have been predicted 1128112141 . 

Fig. [To] shows the FID in single-system dynamics for a typical dot under various field strengths ^ext- The strong field depen- 
dence of T2 demonstrates the significance of the extrinsic hyperfine mediated nuclear coupling up to a strong field (~ 10 T). 
The short time behavior of decoherence by extrinsic nuclear interactions and the behavior by intrinsic nuclear inter- 
action hold well within the relevant timescale for single spin FID. The field and dot-size dependence shown in Fig. [TT] by the 
intrinsic and extrinsic mechanisms, agrees well with the simple form given in Eqns. ([5| and ([6]). For a small QD or in a small 
magnetic field, the extrinsic nuclear interaction dominates, while the intrinsic nuclear interaction dominates otherwise. When 




FIG. 9. A geometric picture for understanding the free-induction decay, (a) The evolution of pair-excitation illustrated by the rotation of a 
Bloch vector and the projected trajectory on the pseudo x-y plane. Direction of the effective pseudo-field are indicated respectively for 
the set of non-local pairs Ka and the set of local pairs Kb. (b) The projection of the Bloch vector trajectories to the pseudo x-y plane for 
pair-excitation driven by extrinsic nuclear interaction. The solid (dashed) line denotes the pair evolution conditioned on the electron spin state 
|+> (|-». As the rotation angle 6 oc t, the distance between the conjugated vectors 6k oc t slI short time, (c) The projection of the Bloch vector 
trajectories to the pseudo x-y plane for pair-excitations driven by intrinsic nuclear interactions. The distance between the conjugated vectors 
6k oc f at short time. 
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the two mechanisms are comparable, the single-system FID begins with e behavior and may cross over towards e decay as 



time increases (e.g., see the curve at ^ext = 12 T in Fig.[TO[b)). The timescale of single spin FID ranges from 0.1 yus to 10 //s 
depending on the dot size and external magnetic field, which agrees well with the experimental observation of 3//s by Greilich 
etal fT09|. 

In the spin echo scenario, as the electron spin is flipped by the 7r-pulse, the transition amplitude by extrinsic nuclear interaction 
A]^ and the hyperfine energy cost E]^ for each pair-flip will change sign after the pulse. Thus, the pair-excitations by the extrinsic 
hyperfine-mediated nuclear coupling will reverse their precession after the pulse and return to the origin at ^ = 2t, disentangling 



the electron spin and the pair excitations (see Fig. 12 a)). So the decoherence driven by the extrinsic hyperfine-mediated coupling 



is largely eliminated in the spin-echo configuration as shown by the calculation in Fig.p3llll25lll27L For the pair-excitations 



driven by the intrinsic coupling, the conjugate Bloch vectors will switch their precession axes which also reverse the entangle- 



ment to some extent but no full recovery is obtained at the echo time (see Fig.[T2[b)). Finally, the electron spin coherence at the 
echo time is derived as, 

i;,,_(2T) = W e-^-'m^^-\K-l^), (7) 

Similar to the analysis for single system FID, the spin echo signal begins with the short-time behavior as ^-(^t/t'^)'^ 

The ensemble spin echo profile is numerically calculated and compared with the single-system FID for a typical QD in Fig. 13 
While it has been a common practice to equal the spin echo decay time Th to the single-system FID time T2 11231 1220ll , the 
two timescales can in fact be significantly diff'erent since the bath dynamics is modified by the pulse control of the electron spin. 
And the spin echo decay and the single-system FID follow diff'erent temporal behavior [Fig. pj]. Th e spin echo decay time of //s 
from calculation 1127111301 is in agreement with the Hahn echo measurements by Clark et al II8III for impurity spins in GaAs 
and by Petta et al 11081 for gate-defined dot in GaAs. 



E. Coherence restoration and protection in the nuclear spin bath 



Protection of the electron spin coherence by active physical control is desired which can result in a better physical qubit before 
the informatic approaches of quantum error correction may be implemented. This is indeed possible for single electron spins in 
interacting nuclear spin baths. 

By a sequence of 7r-rotations of the electron spin, the ensemble dephasing by inhomogeneous broadening as well as decoher- 
ence by extrinsically driven nuclear pair-dynamics is efficiently removed at the classical spin echo time, as shown by the previous 
discussions. Nuclear pair-dynamics driven by the intrinsic interactions is also affected by such control as shown in Fig.[T2jb): 
the two separated trajectories meet again sometime after the 7r-pulse. This intersection signals the disentanglement of the elec- 
tron from the pair-excitation. Surprisingly, even though different pair-excitations have very different precession frequencies hk, 
the trajectory separation 6k is eliminated for all local pairs in the leading order of Bkt dXt = V2t. This leads to a recovery of 
the electron spin coherence as illustrated by numerical evaluation shown in Fig. [14] Remarkably, even when the electron spin 
is flipped after the coherence has completely vanished in single-system dynamics, the coherence may be well recovered at time 
V2t whereas no coherence is visible at the conventional spin echo time 2t 11341 . Thus, in this context, the decay of Hahn echo 
does not mean the irreversible lost of coherence due to the nuclear interacting dynamics. It is simply because the classical spin 
echo time for phase refocusing in ensemble does not respect the quantum behavior of the interacting dynamics in a mesoscopic 




0.0 0.5 1.0 0.01 0.1 1 

t (|Lls) t (|Lls) 



FIG. 10. (Extracted from Ref. [130]) (a) Electron spin coherence as functions of time for various field strengths, (b) The logarithm plot of (a), 
in which the curve for 5ext = 12 tesla is compared with the contribution by the extrinsic nuclear interaction (the dashed line) and that by the 
intrinsic nuclear interaction (the dotted line), respectively. The size of the InAs dot is 33x33x3 nm^ and the nuclear- spin initial state \J') is 
randomly selected from an ensemble at temperature 1 K. The field strength is indicated by the numbers for each curve. 
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FIG. 11. (Extracted from Ref. 1 130|) (a) Field dependence of decoherence times with the inhomogeneous broadening effect excluded. The 
QD is the same as in figure[TO|and the temperature is 1 K. (b) Dot-size dependence of decoherence times with the inhomogeneous broadening 
excluded. The QD size is varied with fixed width : depth : height ratio 33 : 33 : 6. The field strength is 10 telsa and the temperature is 1 K. 
Tile (solid line with circle symbols) - time for FID to 1/^ of its initial value, T2,a (solid line with triangle symbols) - FID decoherence time 
resulting solely from extrinsic hyperfine-mediated pair-flips, T2,b (solid line with diamond symbols) - FID decoherence time resulting solely 
from the intrinsic nuclear spin interaction, and Th (circle symbols) - decay time of the Hahn echo signal. The ^/2T2,b (dashed curve) is plotted 
to compare with the Hahn echo decay time. 
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FIG. 12. (a)and (b) Evolution of the pair-excitations under the single-pulse control, driven by the extrinsic and intrinsic nuclear spin interac- 
tions, respectively. The red solid (blue dashed) trajectories denote the pair evolutions conditioned on the electron spin state |+> (|-». 



bath. When the extrinsic nuclear interaction is significant, the V2t echo can also be weakened by the non-local pair-dynamics 
(see Fig. [14]). 

The V2t coherence echo is observable when the ensemble factor has a timescale longer or comparable to the single spin 
T2 ti me. This is possible with the narrowing of inhomogeneous distribution by nuclear state preparation 1 1131 11 15lll43H15ll (see 
Sec. II F). Furthermore, the unusual echo at V2t turns out to have a more general occurrence in other echo processes. Months 



after its first prediction I13QI 11341 , such echo behavior is observed in NMR experiments in a ^ He gas undergoing Brownian 
motion in a magnetic field gradient 12261 . 

Disentanglement from bath could be the guiding principle for coherence protection with reduced overhead when the bath 
dynamics is more or less understood. In the present case of protecting electron spin from the nuclear spin bath, we need a pulse 
sequence to produce a time where the decoherence from all three sources can be removed. We give here a solution which is a 
two-pulse control. Fig. 14 c) shows that, after a second tt pulse at 3t, the two conjugated paths corresponding to the electron |±) 
states, driven by the intrinsic nuclear interaction, cross again at 4t, coinciding with the secondary spin echo time for the other 
two causes. The electron spin is thus disentangled from entire nuclear spin bath to the leading order at t = 4t (see numerical 
evaluation in Fig.[T4jb)). 

The power of concatenation design of pulse sequences has been shown in the context of dynamical decoupling of quantum 



systems from baths 11331 . Similarly, the control of quantum system for disentanglement from the bath may be enhanced by 
concatenation design. We notice that the pair evolution with the two-pulse control of the electron spin can be constructed 
recursively from the free-induction evolution ^/q , by the concatenation, Uf = Uj_^Uf_^, I = 1,2 I13QII134 I. The process can be 
extended by iteration to any level as shown in Fig. 15 c). Disentanglement from local and non-local pair-dynamics both occur at 
T/ = 2V coinciding with the classic spin echo. The decoherence is reduced by an order of b^r^ at t/ for each additional level of 
concatenation till saturation at the level /o ~ - log2(bT) [130J . Hence, the coherence echo magnitude scales with the echo delay 
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time according to exp (-(T//r/)^^^^) as shown in Fig. ISm), suggesting that short-time decoherence can by arbitrarily suppressed 
with additional levels of concatenation. Numerical calculation further shows that a proper level of concatenation allows the 
protection of electron spin coherence by pulse sequences with interpulse interval as large as ~ yus (see Fig. 15 a)). 

The concatenated control of the decoherence can also be optimized in terms of the number of pulses so that the control errors 
due to imperfections in the controlling pulses is minimized. The invention of the pulse sequences with the minimum possible 
number of pulses to suppress the short-time decoherence to a given order of the pulse delay time is due to Uhrig 11391 in 
considering the qubit decoherence in a non-interacting boson bath. Later, the Uhrig decoherence control was conjectured [ 1401 
and proved [141] to be universal regardless of the bath Hamiltonian. The Uhrig pulse sequences may also be interlaced with the 
concatenated pulse sequences for control not only the pure dephasing but also the longitudinal spin relaxation Ill42ll227ll . All 
these advances help to clear the obstacle of the qubit decoherence in solid environment. 



F. Summary: Quantum dot opportunity 

The three-dimensional confinement in QDs leads to the quantized electronic and excitonic energy levels. As a consequence, 
a QD resembles an atom in terms of the discrete energy levels with long coherence times and well defined optical transition 
selection rules, although it is essentially a mesoscopic system. The atom-like electronic and optical properties have been well es- 
tablished by experiments in the past decades 1 1 991 12QQ1 1205 1 12061 12 1 3 1 1228142301 , including the initial demonstration of quantum 
coherent control [|231h236iI . 

Experiments also showed that, for single electrons in QDs, spin polarization along an external magnetic field can be preserved 
for a sufficiently long time (T^i ~ 20 ms reported for InAs self-assembled dot 11761 ). A major concern has been the inevitable 
cause of transverse decoherence by the lattice nuclear spins in III-V materials. Due to the extremely small energy scales even in 
a strong magnetic field, the nuclear spin bath is of high entropy at experimentally achievable temperature. Ensemble dephasing 
time r* ~ 1 - 10 ns in diff'erent types of QDs II112II202L spin echo decay time Th in the order of yus in gate-defined QDs I108II161I 
and for impurity spins in GaAs 1 181 1, and single spin dephasing time 7^2 ~ 3 yus in self-assembled QDs lil09J have been extracted 
from various experimental approaches, all in agreement with the theoretical analysis of nuclear spin baths 111271 11301 1237L As 
compared with the spin echo decay, FID in single-system dynamics is subject to the additional cause of decoherence from 
the extrinsic mechanism, which is dependent on the external magnetic field. In a moderate field of 1 tesla, T2 ~ 100 ns by 
theory I127II130I , is an order smaller than the spin echo decay time. Polarization of the nuclear spin bath can partially suppress 
this decoherence channel 1 12F, 12061 1238 1. However, a substantial increase of electron spin coherence time would require a 
nuclear polarization over 99%, which can be extremely difficult. Current experimental capability on nuclear spin polarization in 
III-V QDs is in the order of 10-70% l[TT2l[239H243l . 

For spin qubits used as quantum memory, dynamical decoupling schemes can be used to decouple system and environment 
for coherence protection 1 133 . 244 - 2471, which typically requires frequent manipulations of the spin qubits. As compared with 
these schemes designed to deal with a general environment, the available solution to the dynamics of the nuclear spin bath makes 
possible a diff'erent approach aiming at disentanglement of the system from the bath only (their decoupling is a sufficient but 
not necessary condition for disentanglement) |130, 134|. Numerical calculations on realistic QD systems show that designed 
sequence of 7r-pulses on the electron can efficiently preserve its spin coherence up to ~ 100 yus in both single-system and ensemble 
dynamics. As the interval between adjacent pulses can be as long as yus, this disentanglement approach not only substantially 
reduces the overhead but also avoids the problem of unwanted heating of the system from the frequent manipulations required 
by the dynamical decoupling schemes. 
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FIG. 13. (Extracted from Ref. 1 130|) (a) Comparison of the Hahn echo (dashed green line) and the FID (solid red line) signals. The FID 
signal is also shown with the extrinsic hyperfine-mediated pair-flips neglected (dotted blue line), (b) The logarithm plot of (a). The QD is as 
in figure [To] with 5ext = 2 tesla. 
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In the quantum logic control of the spin qubit, is the shortest timescale one encounters. Even with the ultra-fast optical 
manipulation timescale Top ~ 10 ps, T* ~ 1 - 10 ns from FID in a thermal nuclear spin bath is not sufficient to satisfy the current 
fault-tolerant threshold Tcoh/Top ~ 10^ - 10^ P 248H2531 . Considering the efficiency of spin echoes in eliminating the dephasing 
by inhomogeneous broadening, the combination of the desired control action with these coherence protect operations could 
offer a promising route towards fault- tolerant quantum information processing. Efforts are being devoted towards the search of 
universal logic control strategies of coherence protected qubits 12461 . 

An alternative approach is to pre-prepare the nuclear spin bath so that the nuclear field inhomogeneous broadening can be 
squeezed below its thermal value, as suggested by the various nuclear state preparation schemes 1 114lll43H145[|148til51l . The 
resultant enhancement on the T* time can last for seconds or even longer as nuclear spin relaxation is extremely slow. For 
optically controllable electron spin in self-assembled dot, enhancement of T* up to microsecond by nuclear state preparation has 
been achieved experimentally for a spin ensemble 1 146 1, and very recently for a single spin 1 147 |. For electrically controllable 
spin qubits, enhancement of T* of a coupled spin pair to microsecond was also reported in a double-dot configuration 11151 . 

When the inhomogeneous broadening effect is suppressed (by spin echo or bath state preparation), dephasing by nuclear 
interacting dynamics is the limiting factor. Experiments 11091 and theories II 13011 show that, in FID, 72 ~ 0.1 - 10 yus for typical 
self-assembled QDs under a moderate magnetic field, which is sufficiently long to satisfy the fault-tolerant threshold. And the 
spin coherence time can be still elongated further by dynamical decoupling. 

As an alternative, the heavy hole spin in a positively charged QD can also play the role of a qubit carrier. For the p-type 
hole bands, the contact hyperfine interaction vanishes and the hole spin is coupled to nuclear spins through the dipolar hyperfine 
interaction. Theoretical studies shows that the hole-nuclear hyperfine coupling strength is about one order weaker as compared 
with the electron-nuclear hyperfine coupling fSQilSTIl . Furthermore, unlike the isotropic electron-nuclear hyperfine interaction, 
the hole-nuclear coupling is strongly anisotropic. In the absence of heavy-light hole mixing, the hole-nuclear hyperfine interac- 
tion is Ising-like with coupling only between spin components along the growth direction |80, 81 1. Finite hole-mixing effects 
can lead to the coupling between the hole and nuclear spin components perpendicular to the growth direction. In the absence 
of magnetic field, pump-probe and time-resolved photoluminescence experiments on p-doped self-assembled dot have revealed 
hole-spin ensemble dephasing time of 14 ns ISOl . Most significantly, in a magnetic field perpendicular to the growth direction, 
coherent population trapping has been observed for a single ;?-doped dot which suggests a transverse dephasing time longer 
than 100 ns [25 4 J . These experimental findings are consistent with the theoretical studies of the hole-nuclear hyperfine cou- 
pling [sol [811. Since the /7-doped QDs have similar energy level-schemes and transition selection rules to the ^-doped ones, hole 
spins in QDs may be manipulated using optical schemes similar to those for electron spins. For example, high fidelity hole spin 
initialization has already been demonstrated [255 L using similar optical pumping schemes previously adopted for electron spin 
initialization I256M2571 . 

Another kind of solid-state electron spin systems which are under exciting development is NV centers in diamond. An NV 
center in diamond is a defect with a C-C bond substituted with a negatively charged N atom, which has a spin-1 at the ground 
state. NV center spins are a promising candidate for quantum computing for the following virtues 12581 : First, as deep-level 
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FIG. 14. (Extracted from Refs. [1341 and [1301). (a) The electron spin coherence under the control of a short :7r-pulse applied at r = 2 yt/s 
(when the FID signal has vanished), the recoherence at V2t is pronounced while no signal survives at the echo-time 2t. (b) The electron spin 
coherence under the control of a Carr-Purcell pulse sequence. The arrows indicate positions of the tt pulses. The solid blue (dotted red) lines 
are calculated with (without) including the extrinsic hyperfine-mediated pair- flips. The QD is the same as in figure [To| with 5ext = 10 tesla. 
Inhomogeneous broadening is excluded, (c) The projection of the Bloch vector trajectories to the pseudo x-y plane for intrinsically driven 
pair-excitations under the 2-pulse Carr-Purcell control. 
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defects, they have chemical and thermal stability; Second, the spin-orbit coupling in the light C and N atoms is very weak, so the 
spin decoherence by phonon scattering is negUgible even at room temperature |[89l|94l; Third, the natural abundance of isotopes 
with non-zero spin (C-13) is only about 1% and also the hyperfine interaction between the center spins and the bath nuclear 
spins is mostly dipolar which is highly anisotropic and decays rapidly with the distance. Thus the electron spin decoherence 
by the hyperfine coupling is very slow (coherence time > 50 microseconds in natural samples | 89 | and > milliseconds in C-13 
depleted diamond 12591 ): Fourth, the material is optically transparent and the centers are optically active, feasible for optical 
access [W, "921 [93j |95l I2Q4I and coupling with cavities or waveguides f 260H263l . The proposal of quantum computing with 
diamond defects ISTl exploded to a hot research field after the experimental demonstration of electron spin Rabi oscillation 1 90 1 
and two-qubit gates for coupled electron and nuclear spins 1911 in single NV centers. Awchalom group demonstrated coherent 
coupling between a "bright" NV center and a "dark" nitrogen center f93l. Gaebel et al realized strong coupling between an NV 
spin and a nitrogen spin at room temperature 1941 . Lukin and colleagues showed spin echo of an NV spin and observed coherent 
coupling between the electron spin and nuclear spins nearby |97|. In 2007, Lukin group managed to isolate and control an NV 
spin and a strongly coupled nuclear spin |98|. In 2008, Neumann et al claimed multipartite entanglement among C-13 nuclear 
spin near an NV center 12641 . Most recently, Lukin group and Wrachtrup group independently demonstrated readout of an NV 
qubit improved by repetitive retrieval of proximal nuclear spin ancillae | 265 , 266 1. 

In the remaining part of this review, we will discuss the optical manipulation schemes for QD electron spin qubits, and we 
expect that most of them shall be applicable for hole spins and NV centers as well, except the distinctive hyperfine effects on the 
optical processes of the hole spins and of the NV center with its proximity nuclear spin. 



III. PHYSICAL STRUCTURE 
A. Local nodes 

A local node is composed of a few QDs. The specific QD systems of interest include Ini_;cGa;cAs self-assembled QDs and 
GaAs fluctuation QDs. These III-V compound semiconductors have direct band-gaps and thus are suitable for optical control. A 
fluctuation QD is formed by width fluctuation in a narrow quantum well grown with certain procedure (such as interruption for 
introducing interface roughness). This kind of dots has lateral confinement size (> 10 nm) much larger than the growth direction 
size (< 5 nm) and the lateral confinement potential is shallow (usually in the order of a few meV). Thus a fluctuation QD would 
not host many bound electronic states and often just one, which is still subject to ionization due to thermal or optical excitation. 
Nonetheless, the loose lateral confinement makes optical transitions in a fluctuation QD well characterized by selection rules 
resulting from conservation of angular momentum with respect to the growth direction. The large size of fluctuation QDs also 
makes the dipole moment for inter-band transitions to be large and therefore enhances the optical coupling which is useful for 
strong coupling in cavity quantum electrodynamics (QED) and has been utilized in demonstrating optical control of excitonic 
qubits Il231ll232i . A self- assembled QD is formed by the spontaneous nucleation of one material (such as InAs) on the surface of 
a substrate (such as GaAs) which has slightly difl'erent lattice constant. Such a QD is relatively small (with lateral size -10 nm) 
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FIG. 15. (Extracted from Ref. 1 130 |) (a) The electron spin coherence under the /th order concatenation control, / = 0, 4, as functions of the 
pulse delay time r. (b) The logarithm plot of (a). The QD is the same as in Fig. [T0| with 5ext = 10 tesla. Ensemble average is taken at T = 1 K. 
(c) Concatenated sequences of ;7r-pulses flipping the electron spin, represented by vertical bars. 
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and deep confinement potential (in the order of hundreds of meV) which is defined by the off'set between the band-edge of the 
QD material and that of the substrate. As a result, a self-assembled QD could host quite a few stable bound states, providing 
extra flexibility for quantum control. There is no obvious reason to exclude other types of QDs such as those formed in II- VI 
materials and II- VI nanocrystals | 267 1, though they are less comprehensively studied in experiments mostly because of technical 
difficulties such as the requirement of UV lasers, strong charge fluctuation around QD surfaces due to low material mobility, and 
complex defects and impurity centers. The optical control schemes discussed in this review can be applied to the NV centers in 
diamond with some modifications. 

In classical electronic computers, the physical layout of logic circuits is planarly extended. It is not hard to expect the planar 
layout to be used in a quantum computer. In self-assembling growth processes, vertically stacked QDs may be formed to several 
layers | 268 |. The vertical structure, however, is not extendable and furthermore, limits the accessibility by optical pulses. Planar 
distribution is naturally formed both in self-assembled QDs and in fluctuation QDs. Formation of QD molecules or arrays due to 
lateral coupling is possible |195l i269 l. With proper growth art, clusters of QDs with certain patterns may be fabricated 1270- 2721. 
Remarkably, NV centers in diamond may be implanted by ion beams with position precision in the order of 10 nm 1 158> .159l . 
Having such technology within the scope of our consideration, we assume a local node in a quantum computer composed of a 
few (usually fewer than 10) QDs (or impurities), laterally distributed and coupled. The control of electron spins in QDs is to be 
designed after a local node consisting a cluster of QDs has been fabricated and characterized. In this sense, the specific layout 
of a local node, which could be a ring, a line, or any other graphs naturally formed, does not make essential diff'erence. But a 
linearly displaced array would be preferable for its simplicity in coupling and practicality in manufacturing. 

A scalable physical structure of a QD-based quantum computer should have the following features: (1) the QDs are placed in 
an extendable layout; (2) the QDs are connected so that electron spins can be coupled to a common photonic or electronic state; 
and (3) the QDs are individually accessible so that the electron spins are individually controlled. On the one hand, the spatial 
resolution of near-field optical devices is still not high enough to identify each QD in a cluster. On the other hand, within the 
limit of current technologies, it is impossible to control the growth of QDs so that they are almost identical. Diff'erent sizes and 
shapes of the QDs would make the exciton transition energies in diff'erent QDs diff'erent. In this way, the near-field optics and the 
fingerprint transition frequencies of diff'erent QDs may be combined together to individually address each dot or to selectively 
couple a pair of them. The coupling between spins in general is mediated by virtual tunneling between diff'erent QDs which 
may be activated by virtual optical excitation of excitons in the presence of extra electrons which bear the spins (More details 
will be discussed in later sections). For NV centers in diamond, such mechanisms are not yet considered, but other schemes 
may be applicable such as coupling through virtual excitation of cavity modes |[64|, or hyperfine interaction with nuclear spin 
baths dpO I. 



We show an example in Figure 16 a working module of 7 qubits in a linearly displaced array of QDs. With the QDs of 
the size around 50 nm and about 20 nm apart, the 7 dots may be addressed with two micro-lenses attached to, e.g., optical 
fibers, with resolution of about 0.2 jim. The QDs addressed by the same optical fiber are distinguished by their signature 
transition frequencies. Using optical pulses with diff'erent frequencies and polarizations, each dot and each adjacent dot pair 



may be (virtually) excited and various single- and two-qubit gates and initialization could be realized, as illustrated in Fig. 17 



To realize arbitrary single-spin rotation, a magnetic field is applied along the in-plane-direction (denoted as the z-direction). 
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FIG. 16. The physical structure of a local node consisting of seven electron spins in a QD cluster. Each QD is identified by a near-field micro- 
lens as well as by its characteristic transition energies. The numbers are used to label the qubits for factorizing 15 with the seed number a = 4 
(top row) and a = 13 (bottom row), respectively. The numbers in the parentheses label the qubits at the end of computation (see Sec. |VIII A| 
for details). 
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We assume that only neighboring QDs are coupled, for simplicity. Coupling between spins in farther separated dots is to be 
accomplished by recursively using the nearest neighbor coupling, which increases the number of gates by an amount in the order 
of the number of dots between the two ends. To minimize the number of operations between separated qubits, the quantum 
algorithms are compiled by using optimum labeling of the qubits. Taking into account the size of the micro-photonic structures 
(including cavities and waveguides connecting local nodes), a working module of about 10 QDs could occupy an area of about 
10 //mxlO yum, so a quantum chip of size 10 cmxlO cm can in theory accommodate 10^ qubits. 

For the purpose of addressing and optical control of a single QD, several near-field optics technologies have been available, 
such as micro-optical masking and micro-fiber optics |273| . For a cluster of QDs to function under optical control, however, there 
are still non-trivial technical challenges, including at least fabrication of QD clusters with energy levels and inter-dot couplings 
falling in the desired parameter ranges, assembling of micro-lenses on the surface of QD clusters, and designing, shaping, and 
controlling complex laser pulses. The micro-lens technology has been widely used in the digital imaging industry (micro-lenses 
of similar sizes are routinely used to focus light into individual pixels in commercial digital cameras). 



B. Distributed structure 



The electron spin qubits are distinguished by the diff'erent optical transition frequencies of the host QDs. Thus, the dot 
density shall be rather dilute so that each laser spot contains only a small number of QDs (~ 0(10), considering the typical 
inhomogeneous broadening of the excitonic linewidth of ~ 10-100 meV, and the requirement of ~ 0.1 - 1 meV frequency 
separation for optically addressing individual QDs. Therefore, the optical approach predetermines that a local node can only 
have a limited number of qubits. In order to scale up, a distributed architecture could be a solution fT,'214]. In such structures, 
clusters of QD electron spin qubits form quantum nodes where logical operations can be performed locally, and connections 
between clusters are through quantum channels in which the flying qubits take information from one place to another. 

The single photon wavepacket is an ideal candidate as the carrier of flying qubits, being widely used in quantum cryptog- 
raphy 12751 and linear optics quantum computation 12761 . The qubit can be encoded in the photon-number subspace II277II or 
polarization subspace | 276 |. While single photon propagation in free space is un-channeled and inefficient, optical waveguides 
in semiconductors and optical fibers provide directional channels. 

In the distributed architecture of optically controlled spin quantum computation, flying photons in waveguides/fibers are 
responsible for integrating the distributed stationary spin clusters into a globally functioning quantum computer. This requires 
quantum interfacing between single electron spins in QDs and single photons in waveguide. As mentioned in Sec. |Il| QD 
electron spins interact with optical fields via the intermediate states of trions. Such interface at the single photon level requires 
strong light-matter interactions. As a QD has a fixed optical transition dipole moment which is limited by its size, one way to 
have such strong light-matter interaction is to confine photons in optical cavity structures with small volumes. 

Micro-cavities can be realized in a number of ways in semiconductor structures. We list below the essential properties of a 
few representative ones: 



1. Microdisk - See Figs. [18] and [19] for schematic illustrations of the geometry. Light is confined by total internal reflection 
in the inner wall and the confined modes are known as the whispering gallery modes. Quality factor Q ~ 0(10"^) in III-V 
materials and ~ O(IO^) in polymer; mode volume V ~ 6(Ao/n)^ where n is the refractive index of the material and Aq 
the wavelength of cavity mode [278] . Strong coupling regime for GaAs fluctuation QD embedded in micro-disk structure 



ORKKY 
Raman / AC Stark 




Cooling 



FIG. 17. The energy diagram of various optical processes for quantum gates and initialization: "Raman" for single-spin control or optical 
pumping and measurement of a single spin, "ORKKY" for two-spin control, "AC Stark" for transient shift of energy levels to relize selective 
resonance, and "Cooling" for initialization of spins via phonon bath (phonon emission as wavy lines). 
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Quantum Dot 



FIG. 18. Left: coupled cavity and waveguide structure formed by point and line defects in 2D photonic bandgap crystals. Right: waveguide 
coupled micro-disk cavity etched on chip. Layers of QDs can be embedded in the matrix slab where the cavities are formed. 




Nanoparticle 



FIG. 19. Left: fiber coupled micro-disc cavity. Middle: fiber coupled micro-sphere cavity. Right: fiber coupled micro-toroid cavity. 
Semiconductor nanoparticles, e.g. diamond with NV centers, can be adsorbed onto the latter two types of cavities. 



similar to the illustration in Fig. 19 has been achieved 12791 . Because of the large dipole moment of the fluctuation dot, 



the measured cavity-dot coupling constant gcav -0.2 meV signifies strong coupling. 

2. Defect cavity in 2D Photonic Crystal - Two dimensional photonic bandgap crystals are an ideal structure to form a cavity 
resonator 12801 12811 . Propagation of light in the plane has a forbidden bandgap for carefully designed periodical arrays of 
air-holes drilled on the 2-D slab. As shown in Fig. [18] by forming a point defect in the 2D array of air holes, light can be 
almost perfectly confined in the plane of the slab if its frequency lies in the forbidden bandgap. The vertical confinement, 
achieved by total internal reflection at the semiconductor-air interfaces, is imperfect, in that light with small in-plane 
wavevectors can leak out of the top and bottom. Vertical leakage can be greatly suppressed by proper engineering of the 
defect 12821 [2831. 2 ~ 6 x 10^ and V ~ 1.2(^o/^)^ ~ 0.072 jim^ have been achieved f282','5S51. Theoretically analysis 
shows that Q-factors greater than 2 x 10^ are realizable by optimizing the structure | 283 |. The matrix of a 2D photonic 
crystal can be either silicon or III-V compounds 12841 . Strong coupling with single self-assembled InAs QDs has been 
demonstrated II284II285L where gcav ~ 0.1 meV. 

3. Micro-pillar - Light is vertically confined by distributed Bragg ref elector mirrors and horizontally by total internal reflec- 
tion. Q ~ 10"^ and V ~ jdw? have been achieved 1286, 2871 . Strong coupling with single self- assembled InAs QDs has 
been demonstrated 1286112881 . where gcav -0.1 meV. 

4. Epitaxial cavity - Vertical confinement is by distributed Bragg refelector mirrors and horizontal confinement by thick- 
ness variations, similar to the confinement principle of the fluctuation QD 12891 . 2 = 3x10^, y~ jivo? have been 
achieved I289J. 



5. Silicon microsphere - WGMs are confined by total internal reflection. Q exceeding 10^ and V ~ 1(9 [iw? have been 
achieved |290, 291 1. Nanocrystals (such as CdSe nanocrystals |292] and diamond nano-crystals with NV centers [260] ) 
deposited on the surface are usually used for coupling with the cavity photons. 



6. Microtoroid - See Fig. 19 for an illustration of the geometry. WGMs are confined by total internal reflection, g ~ 10 
is achieved with principal diameter Z) ~ 100 yum and the minor diameter d ~ jim 12931 . Theoretical analysis shows the 
possibility of realizing micro-toroid with Q exceeding 10^ and V ~ 0(10) jdvo? 12941 . 
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The high quaUty factor allows photons to be confined for a sufficiently long time inside a cavity and the small mode volume 
makes possible a large intra-cavity electromagnetic field from a single photon. Both features are critical for the strong dot-cavity 
coupling, defined by the criteria that a single cavity photon can induce > In Rabi oscillation of an excitonic transition in a single 
dot within the lifetimes of the cavity photon and the exciton. Strong coupling has already been realized in several cavity-dot 
systems 1 260. 279, 284TI286I [288, 292 1. With micro-cavity being a playground for strong coupling of a single photon with a 
single dot, the interface between a flying photon qubit in a waveguide and a stationary spin qubit in a QD can be achieved by 
evanescent coupling between an optical waveguide/fiber and a micro-cavity containing the QD (see Figs. [18] and [19]). A single 
photon wavepacket propagating in the waveguide can excite a cavity photon which then significantly influences the QD spin 



dynamics, e.g. through the optical Raman process via the trion states (see Sec. |IIB| ) 11191 . In addition to the exciting advances 
in improving the Q-f actor and reducing the mode volume of semiconductor micro-cavities for strong coupling, many other 
key ingredients towards the construction and control of such a dot-cavity- waveguide coupled structure have been progressively 
achieved in laboratories. These include the high-efficiency coupling between micro-cavities and optical waveguides/fibers 12951 - 
I3Q1L precise control of intra-cavity location of QDs [ 285l 13021 1303IK fine-tuning of cavity modes into resonance with a given 
QD transition 12851 13021 1304H306L electrically controllable charging of intra-cavity QDs 13071 . and coherent resonant driving 
of QD excitonic transitions inside a cavity by an external laser 12301 13081 . 



With the possibility of spin-photon interfacing in dot-cavity- waveguide structures, we are able to outline the construction of a 
distributed architecture for scalable quantum information processing in an integrated semiconductor platform composed of QDs, 
optical micro-cavities, and optical waveguides/fibers (see Fig. 20 ) | 309 |. In this network structure, a local node contains a limited 
number of charged QDs, distinguishable by their optical frequencies. The stationary spin qubits form a basis for the quantum 
memory and quantum logic modules. Optical waveguides/fibers connect distributed nodes with single photon wavepackets as 
the ffying qubits. Micro-cavities off'er the playground for the strong interaction between the two types of qubits. Control of a 
local node will be the focus of Sec. IV Control of the interfacing between single spins and single photons will be the focus of 
section r 



VII With the recent progresses in coupling NV centers in diamond with photons in cavities and waveguides Il260ti263 



similar distributed quantum computing with NV centers is also foreseeable. 
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FIG. 20. Distributed quantum information processing in integrated semiconductor structures, (a) Schematics of a distributed quantum com- 
puter where communications between computation modules are mediated by single photons in optical waveguides/fibers, (b) The spin-based 
computation modules on a chip, controlled by ultrafast optics. 
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IV. ONE- AND TWO-QUBIT OPERATIONS 

It has been established that universal quantum computation can be accomplished by a set of single-qubit gates and one kind of 
entanglement gates such as CNOT, VSWAP, or controlled phase-shift gate 1117211183111841 . While there are many alternatives, 
we follow a more traditional and scalable approach of designing a system which has negligible interaction between qubits when 
the system is idling, and creating transient interaction between two qubits during the operation. 

General speaking, there are two strategies for controlling electron spins in QDs, namely, the direct and the indirect control. As 
the terminology suggests, the direct control manipulates the electron spins directly which may be realized by rotation of a single 
spin with an AC magnetic field from, e.g., a microwave or a pulse DC magnetic field ll66l[88ll90l[97ll98lll61L and by coupling 
two spins via exchange interaction mediated by virtual tunneling between QDs switched by, e.g., gate voltages 1631 1651 182L 
The direct control schemes are more applicable to electrically defined QD systems made possible by the feasibility of in-situ 
electrical gates. There is no fundamental obstacle to integrate electrical gates to self-assembled or fluctuation QD systems and 
to apply the direct control of electron spins. But the much smaller size of these QDs as compared with the splitting-gate defined 
ones and the much stronger confinement (or much less tunneling probability) cause non-trivial technical problems. The indirect 
control schemes are based on indirect coupling between electron spin states mediated by virtual excitation of auxiliary energy 
levels which are modified by (or conditioned on) the states of the spins f64l [68ti79lll811[T82 l. Usually, such intermediate states 
are excitons which may be excited by optical pulses. The optical control of electron spins is limited to direct-gap semiconductors 
but otherwise have a great deal of merits particularly due to the energy scale cascade. As bandgaps in semiconductors are usually 
larger by orders of magnitude than the electron spin Zeeman energy, the control of the electron spins could be made much faster 
by using excitation cross the large bandgap. For comparison, electrical control of spins on the nanosecond scale is already 
the state-of-the-art technique 11081 1161L while optical manipulation of electron spins can be completed on the picosecond 
timescale ll78l[79l fT8n[T82ll . 

In essence, all schemes of optical control of the electron spin states (641 168ti76l 11731 are realized by the Raman processes 
where the virtual excitation of excitons plays the central role. In one-qubit operations, the intermediate states are excitons 
with one excess electron, i.e., the trion states 168117211^. In two-qubit operations |64l[68M7Tl[73^-76, 173], the intermediate 
states are excitons with two excess electrons. The eff'ective exchange interaction between two spins may be induced by the virtual 
excitation of one exciton tunneling back and forth between two dots, which is similar to the RKKY process and is dubbed optical 
RKKY interaction to indicate the role of the optical excitation [73 -75 1. Alternatively, even if there is no tunneling between QDs, 
the interplay between Coulomb interaction and the spin blocking eff'ect would induce eff'ective interaction between two separated 
spins when two excitons are (virtually) excited in two dots L68i-i71t .173 J . The spin interaction may also be mediated by virtual 
photon exchange during the virtual excitation of excitons (6411761 . 

More details about the optical control of one-spin and two-spin gates follow. Since the essential physics of various optical 
control schemes is similar, we will focus on two specific examples, namely, Raman control of single spins, and optical RKKY 
control of two spins, as an illustration of the physics and the operation conditions. 

A. Single-spin rotation by Raman process 

Arbitrary single-spin rotation can be realized via adiabatic Raman processes mediated by trion states (72l 11691 . We will 
illustrate the basic idea of such a scheme and discuss the non-adiabatic generalization and certain limits of the method. It is 
worth pointing out that the stimulated Raman adiabatic passage (STIRAP) |310| involving only dark states could not realize 
an arbitrary rotation but only a spin flip from a known initial state. The adiabatic Raman processes in arbitrary spin rotation 
involves both dark states and bright states |72|. 

The control of an electron spin relies critically on the spin states of the excited excitons which are determined by the optical 
selection rules in the semiconductor QDs. In the III-V semiconductors of zinc-blende crystal structure, the selection rules of 
optical excitations at the band edge are well defined by the angular momentum conservation in terms of the electron spin in the 
conduction band and the hole spin in the valence band, which have angular momentum ^ = 1/2 and / = 3/2, respectively. In 
QDs of large lateral sizes and strong confinement in the growth direction (defined as z-axis), the hole states will be split into two 
sets of degenerate states, with spin states = ±3/2 and ±1/2, respectively, designated as heavy and light holes, respectively, 
according to their eff'ective mass along the z-axis. When the lateral sizes of the QD are much larger than the confinement size in 
the growth direction, the mixing of diff'erent angular momentum states by the lateral confinement is small (STI . Thus the optical 
excitation is restricted by the angular momentum conservation along the z-axis. Now if the controlling optical pulses are applied 
normal to the sample surface, the conservation of the angular momentum about the growth direction makes it impossible to flip 
the electron spin along z-axis and thus impossible to complete an arbitrary quantum operation, unless the light beam is incident 
with an angle | 311 1 or the symmetry is broken by a magnetic field with a non-zero in-plane component. Since in the near-field 
optics, the incident light is usually normal to the surface, we need a static in-plane magnetic field applied (whose direction is 
defined as x-axis). Under the strong magnetic field, the electron spin states are split into two Zeeman levels | + x) and | - x) with 
energy ±0Jc/2, respectively. We use these two states as the basis |0) and |1) of a qubit. The Zeeman splitting ojc is in the order 



24 



of 0.1-1 meV under a magnetic field of a few Tesla. In GaAs fluctuation QDs, the hole states are still near-degenerate, since 
the large heavy-Hght hole splitting (Ahi is tens of meV) makes the hole spin splitting ~ col /^w negligible even under a field as 
strong as a few Tesla, so the hole states can still be defined by the magnetic quantum number as | ± 3/2), which will also be 
denoted by hollow arrows as | |i) and | ff), respectively. In this case, the optical transitions can be separated by the selection rules 
even when the energy difference due to the electron spin splitting is relatively small. Note that the near-degeneracy of the two 
trion states is not a necessary condition in the schemes of single spin rotation to be discussed below. Actually, in self-assembled 
dots where hole states can be split by an in-plane magnetic field due to the large heavy-light hole mixing, the splitting between 
the trion states can be exploited to separate desired transitions from unwanted ones by energy diff'erence. 

To illustrate the essential physics underlying the control process, we ignore the high-lying excited states and model the system 
by four states consisting of the two split electron spin states and two trion states |^+) which are formed by two electrons in the 
singlet state | tsl) and one hole in the spin state | ff) or | respectively. According to the angular momentum conservation about 
the z-axis, the selection rules of the optical excitation, as depicted in Fig. 21 (a), is such that a light with circular polarization cr- 
will induce the transition from the electron states | -h 1/2) or | - 1/2) to the trion states \t+), respectively. 

So, under the excitation of an optical pulse with cr^ polarization, the dynamics is reduced to a Raman rotation in the A-type 



three-level system [as shown in Fig. 21 (b)], which is governed by the Hamiltonian 



H = cocs, - Q*(OI + 1/2X^+1 - Q(0k+)(+l/2| - AluKu 



(8) 



where Sx/y/z is the electron spin operator along the x, y, or z direction, respectively, Q(0 is Rabi frequency of the laser pulse in the 
rotating frame, and A is the detuning of the laser relative to the trion state |^+). To eliminate the dynamic phase associated with 
Zeeman splitting, the quantum operation should be done in the frame rotating together with the free spins under the magnetic 
field. By the transformation S = exp (icOctSx) , the Hamiltonian in the precessing frame is written in the matrix form 



H 






-Q(0^"^"^^'/V V2 -Q(0^^'"^^'/V V2 



-A 



(9) 



in the basis of _ x), ^"^"^^^l + x), and |^+). 

For large detuning (|A| much greater than the bandwidth of the optical pulse and the Rabi frequency), the so-called adiabatic 
approximation is justified and thus under a standard canonical transformation, the off'-diagonal terms in the Hamiltonian between 
the electron spin states and the trion state are eliminated up to the second order of the Rabi frequency. So the transformed effective 
Hamiltonian is approximated as 



eff 



|Q(0lV(2A) 
-|Q(0P^"^"^^7(2A) 




-|Q(0P^^^"^^7(2A) 
|Q(0lV(2A) 







-A - |Q(0|VA 



(10) 



This is equivalent to a magnetic field with strength |Q(0P/(2gyUfiA) precessing in the x-y plane with the angular frequency ojc 

of the electron energy |D(0P/(2A) contributes only a trivial global phase-shift and 



(the time-dependent optical Stark shift 13111 
can be ignored). In cases that the optical pulse is much shorter than the spin precession period, the effective magnetic field 




FIG. 21. (a) Optical selection rules for the electron-trion transitions, (b) The Raman process in the A-type three-level system formed by the 
electron spin states and the trion states connected by a cr+ -polarized laser pulse. 
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becomes an instantaneous pulse which can be controlled in the femtosecond timescales (c.f. the GHz limit on the control of 
external magnetic field). For t ^ nnlicoc) in is an integer, the magnetic field pulse is eff'ectively along the z-axis, and for time 
around {n + l/2)7i/oJc, the instantaneous rotation of the spin is eff'ectively along the j-axis. Thus we have two SU(2) generators 
which can be combined to complete an arbitrary rotation of the electron spin. To be specific, a spin rotation with angle y along 
an axis defined by the Euler angles (a,l3), denoted as R(a,l3, y), can be realized by at most 3 elementary rotations along different 
axes in the x-y plane, for 

(71 71 71 \ (71 \ (71 71 \ 

where the three elementary rotations on the righthand side of the equation can be performed in turn at ^ = -7i/(2ajc) -\- j3/ajc, 
jS/oJc, and -\-7i/<jJc + P/oJc^ which can be completed within half the precession period of the electron spin. For an electron with 
Zeeman splitting of 1 meV, the duration of an arbitrary operation of a single spin is less than 2 ps. 

To complete a finite rotation of the spin at a time much shorter than the precession period, the Rabi frequency should be 
large, which, however, could demolish the condition for the adiabatic approximation. For instance, if the Zeeman splitting is 
0.1 meV and the detuning is 1 meV, to complete a tt rotation within a duration one tenth of the precession period, the Rabi 
frequency is required to be about 3 meV, greater than the detuning, which makes the adiabatic approximation unjustified. In 
fact, to validate the adiabatic condition for instantaneous operations, it is required that ojc ^ A, which means slow operations 
for reasonably large detuning (A < 10 meV). Alternatively, by shaping the control laser pulse, the rotation can also be operated 
in the non-adiabatic regime or in a non-instantaneous manner. No matter how intense and how fast the optical pulse could be, 
the spin rotation in general cases, however, is still limited by the precession period, as can be seen from the equation for the spin 
polarization as 

dtisz) = ojcisy) - dtpu (12) 

where pt is the population of the trion state. For a complete operation, the residue population of the trion should be zero, so the 
change of the spin momentum along the z-direction is |^<^^)| = \(jJc ^^{Sy)dt\^ < ojct, which cannot be faster than the precession 
under the static external magnetic field. 

In summary, under a moderate external magnetic field (< 10 T and 6g>c ~ 1 meV), an arbitrary spin rotation can be ac- 
complished well within 10 ps by up to three ultrashort optical pulses (with simple shape and large detuning) or by one pulse 
(with engineered shape). Remarkably, control of single electron spins 11791 12361 [31 2>| and hole spins fW\ in QDs in picosecond 
timescales have been recently realized in experiments. Using the optical control, Yamamoto's group has demonstrated spin echo 
for an ensemble of impurities in GaAs 11811 and for a single QD spin 13131 . Greilich et al have realized optical rotation of an 
ensemble of QD spins along arbitrary axes 11821 . 



B. Two-qubit gates by optical RKKY interaction 

To implement two-qubit quantum gates, optically induced RKKY interaction (ORKKY) between electrons doped in QDs has 
been proposed to couple two spins, laterally via continuum excitons 1731 |74l or vertically via discrete states |75|. The ORKKY 
interaction mediated by continuum excitons ||73l|74l, due to the extension nature of the continuum states, is less controllable 
in selectively coupling certain spins. The ORKKY interaction between QDs vertically stacked |75|, on the other hand, is not 
applicable to a scalable system with a planar layout. For the sake of scalability, we consider to employ the discrete excited states 
in laterally coupled QDs to induce the ORKKY interaction, which, as described below, can be controlled to selectively couple 
spins in designated adjacent QDs. The ORKKY interaction is by nature a Raman process with the ground states formed by the 
two-electron spin states and the intermediate states by excitons charged with two excess electrons. 

The physical process of the optically induced RKKY interaction via discrete intermediate states is depicted by the Feynmann 



diagram in Fig. 22 (a). The qubits under controls are the spins of two electrons in the ground states of two neighboring QDs, 
denoted as Si and S2. When the electrons are in the ground states, the inter-dot tunneling is negligible and therefore the two spins 
have no direct exchange interaction. The basic elements of the ORKKY process for mediated interaction between the two qubits 
can be described as: 

1. an incident optical pulse excites a direct electron-hole pair into the excited electron and hole levels in, say, QD 1, denoted 
as 

2. the optically excited electron will interact with the electron spin Si with the strong exchange interaction (Ai ~ 5 meV in a 
typical In As QD), while the electron-hole exchange interaction is negligible in comparison; 

3. with a strong quasi-cw optical field applied, the excited electron levels in the two QDs can be tuned into resonance by 
the optical Stark eff'ect, and thus the electron in the excited state \ei) can resonantly tunnel into the second QD with the 
tunneling rate in the order of 10 meV for two QDs separated by 15 nm, while the hole tunneling can be neglected due to 
the stronger confinement and off'-resonance condition; 
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4. having tunnelled into the excited level of the second QD, the optically excited electron can exchange spins with the excess 
electron in the second QD via the strong exchange interaction; 

5. after the exchange interaction with the second qubit electron, the electron in the excited state can tunnel back into the first 
QD; 

6. the electron back to the first QD can be recombined with the hole by emitting a photon back into the optical pulse, leaving 
the two excess electrons an eff'ect of indirect spin exchange. Since the laser frequencies can be adjusted to selectively 
excite the exciton in one QD and to selectively shift the level in another QD by the AC Stark eff'ect, a pair of adjacent QDs 
can be coupled on demand. 

Since the Coulomb and the exchange interactions could be strong for electrons in the discrete states, as compared with the 
optical Rabi frequency, it is better to treat the interactions non-perturbatively by first exactly diagonalizing the states of charged 
excitons, and use the eigenstates as the basis for calculating the optically induced indirect exchange interaction. In fact, when the 
intermediate states are discrete eigenstates of the Coulomb interaction, the ORKKY process becomes equivalent to the Raman 
process in the multi-level systems with the two-spin states as the ground states and the charged exciton states as the intermediate 
states. Then the optical excitation can be treated as a perturbation, similar to the case of the single- spin rotation, and the adiabatic 
approximation may also be adopted when the detuning is large. 



As shown in Fig. 22 (b), we assume the following characteristic energies for two neighboring QDs of diff'erent sizes, which 
can be realized by properly adjusting the QD sizes, the separation, and the composite x in Ini-^^Ga^^As in typical cases under 
consideration: The diff'erence in transition energy between the neighboring dots is dn ~10 meV, the distance between the 
ground and excited electron states £^oi ~ 50 meV, the tunneling strength between neighboring excited states ~ 10 meV, 
and the Coulomb exchange energy between electrons /I1/2 ~ 5 meV. When the electron-hole pair is excited into the excited 
single-particle levels in the smaller QD, the diff'erence between the binding energy of the direct exciton (when both the electron 
and the hole are in the smaller dot) and that of the indirect exciton (when the electron has tunnelled into the large dot), which is 
-10 meV, can already compensate most of the energy diff'erence (^12) between the excited single-particle levels in the two dots. 
And if desired, a further optical pumping can be used to fine-tune the energy levels in the dots so that near-resonance tunneling 
can occur between the excited levels for electrons in the two dots. The ground single-particle states for electrons and the hole 
states can be virtually taken as localized due to the relatively high barrier. 

To illustrate the basic idea of using the discrete exciton states in laterally coupled QDs for ORKKY interaction between the 
two qubit spins, we assume that the hole is (virtually) created by the optical pulse only at the first excited level in the smaller QD 
and the dynamics is essentially determined by the interaction between the two qubit spins and the optically excited electron. The 
role of the hole is then just to impose the renormalization of the electron energies. As mentioned above, both the tunneling and 
Coulomb exchange are in the strong coupling regime, they should be considered non-perturbatively. So we can first diagonalize 
the Hamiltonian with the exchange energy and electron tunneling fully included. As the tunneling is spin-independent, we 
consider only the bond state of the QD molecule and treat the anti-bond state as in far off'-resonance (since it is about 10 meV 
above). Thus the relevant spins are qubit 1 (si), qubit 2 (S2), and the spin of the bond molecular state 3 (S3). The Hamiltonian of 
these three spins is 

H = 2/liSi • S3 -\- 2/I2S2 • S3 -h 2ojc(s\ + 4 + 4)' 
where Ai and A2 are the exchange energy in each dot. For simplicity, we assume Ai = A2 = A. The eigenstates and the 




FIG. 22. (a) Feynman diagram for optical RKKY interaction. The straight arrows are electron or holes lines, the wavy arrows are photon lines, 
the tunneling is the cross vertices, the Coulomb exchange is the scattering vertices, and the AC Stark field is the dotted lines with cross, (b) 
Energy level schematics for two-qubit gates. The two-spin states are split by the magnetic field, and the exciton charged with two electrons 
(the hole is understood) is split both by the Coulomb exchange energy and by the magnetic field. The dotted arrows are optical excitation for 
a SWAP gate, which together with the excitation represented by dashed arrows can also accomplish the phase gate. 
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corresponding eigen energies can be worked out as 
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as schematically shown in Fig. [22] (b), where the up/down arrows indicate in turn the spin states of qubit 1, qubit 2, and the 
electron created by optical excitation, quantized along the external magnetic field direction. Notice A - 5 meV» 2ojc ~ 1 meV. 
The two-qubit gates can just be realized by the Raman passages between the two- spin ground states mediated by the charged 
excitons formed by one hole plus three electrons in the spin states shown above. To be specific, we discuss two frequently used 
two-qubit gates as below. 



1. SWAP and ^/SmP gates 



A SWAP gate interchanges the states of two qubits. Its matrix form is shown in Fig. 45 In the Raman process, it just flips the 



two states: | tvl) and | The intermediate state connecting these two states by Raman process is 15", |, |) = -^(| tvlT) - 1 vlTT)) 
(henceforth the hole state in the exciton has been omitted for the sake of simplicity), which is optically coupled only to the 
two-spin singlet state -^(| tvl) - I vlT))- By performing a Inn Rabi rotation between the intermediated exciton state and the 
singlet ground state, a pure phase-shift will be induced to the singlet state, whose value is determined by the detuning. When 
such a phase-shift is controlled to be tt (which can be obtained by, e.g., setting the rotation angle to be 27i and the detuning is 
zero), the two states | ti) and | i^) are just flipped and the SWAP gate is realized. In practice, the optical pulse can be tuned 
off -resonance from the transition to suppress the spontaneous decay. With proper polarization, the optical pulse doesn't excite 
the nearby state \S,^,-^} and all other transitions are at separated by an energy of A. So this operation can be accomplished in 
a period of time ~ 10 x (InA'^) ~ 10 ps. Considering the optical Stark pulse used for tuning the resonant tunneling, a two-qubit 
SWAP gate in such a scheme would require two optical pulses of duration of about 10 ps. 

By choosing a proper detuning, the phase-shift can also be controlled to be ;r/2, and then the SWAP gate is performed halfway, 
or the VSWAP gate is realized, which has the matrix form 
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(15) 



The VSWAP gate can generate entanglement between the two spins and can be used to realize the controlled phase gate and the 
CNOT gate. 



2. Controlled phase gate 



Another two-qubit gate which can be used in lieu of the CNOT gate for universal quantum computation is the controlled phase 
gate which induces a phase-shift of the target qubit depending on the state of the controlling qubit. The general matrix form of 
the controlled phase gate is 
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(16) 



where the subscript of the phase-shift indicates the corresponding qubit states, Rx(0) is the single qubit phase-shift of 6, and 
= 011 -h - - 0^1 . When (p = 71, the phase gate can be transformed into the CNOT gate by applying certain single-spin 
gates, which are usually assumed much easier than the two-qubit gates. 
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As shown in Fig. 22 (b), to obtain the phase-shift, an optical pulse is applied to couple the state | H) to the exciton state 
|P, |, -|) = I iii). The pulse shifts the states by optical Stark effect and induces the phase shift (pn. This pulse will also couples 

the electron state (| U> + I U))/ V2 to the exciton state |P, |, -^) = (| iU) + I TU) + I UT))/ V3, inducing a phase shift to 
(I vlT) + 1 Tvl))/ V2 and thus the rotation between | i^) and | t>l). To obtain a pure phase gate, another pulse (with the same energy 
and polarization as in the SWAP gate) coupling the states (| U> - I U))/ and 15", ^, |) = (| UT) - I UW/ ^ can be used to 
compensate the rotation. Finally, the conditional phase-shift is = 0|| - - (p^i = Similar to the case of the SWAP 

gate, a phase-shift gate can be realized with three pulses (including the optical Stark pulse) of duration of about 10 ps. 



C. Issues to be considered 



In realistic cases, there are several issues which could degrade the fidelity of the quantum gates. 

The first one is the relaxation of the intermediate state by spontaneous emission. In the stimulated Raman adiabatic pas- 



sage I310L the spin states can be flipped without populating the exciton state, which, however, cannot perform a general quan- 
tum gate since such a passage depends on the initial state of the spin. To have a general quantum gate which transforms a 
spin independent of the initial state, both the dark and bright states should be employed. To suppress the spontaneous emis- 
sion and other scattering processes, it is preferable to have large detuning so as to minimize the population of the intermediate 
state |[68l[72l [T69l[T73ll . 

Another important eff'ect aff'ecting the fidelity is the imperfection of the selection rules and the hole-mixing CTl • In realistic 
cases, the QDs would never have perfectly symmetric shape and thus the lateral confinement could cause mixing between states 
of diff'erent angular momentum (such as the heavy-light hole mixing), which is worse for the excited states involved in the two- 
qubit gates. Such eff'ects, however, only induce systematic errors or unwanted dynamics to the quantum gates designed for ideal 
condition. In principle, the shapes, polarizations, and timings of the controlling laser pulses can always be readjusted once the 
realistic system parameters have been measured. It could also be possible to design pulses of certain robustness against small 
deviations in the system parameters I ll20ll314ll , and a scheme of using chirped pulse to implement quantum gates robust against 
the mixing eff'ect has been proposed ITTI. 

For QDs with reffection symmetry with respect to the growth plane, the imperfect selection rules may be tolerated by re- 
designing the polarization of the control light. In general cases, especially for small QDs which have irregular shapes, the hole 
mixing may be used as a resource for quantum control. In previous discussions in Sec. IV A a static magnetic field applied along 
a direction other than the growth direction has been required to break the rotation symmetry so that an arbitrary rotation of a 
single spin is possible. When the conservation of angular moment with respect to the growth direction is not perfect and thus the 
"forbidden" transitions would be made partially "allowed" due to the hole mixing, an eff'ective magnetic field along an arbitrary 
direction for the spin could be induced by a properly polarized light beam through the AC Stark eff'ect and an arbitrary rotation 
of a single spin could be realized even without a external magnetic field applied. More discussions about the hole mixing and its 
eff'ects on the optical control are given in Appendix [A[ 

The third problem with the scalability of the quantum computation is the complexity in the energy level structure of multi- 
dot systems. The analysis and characterization of the many levels with a number of excess electrons require much eff'ort and 
furthermore 1315143191 , the optical pulses applied to a desired transition would inevitably aff'ect the other transitions in QDs 
nearby, making the pulse design very demanding when the system becomes large. In femtosecond chemistry, learning algorithms 
haven been developed to design sophisticated pulses for controlling the complex atomic and molecular dynamics Il320ti322l . We 
expect the quantum learning algorithm be a powerful tool to deal with the design complexity in multi-dot systems. 

With recent experimental progresses demonstrating the feasibility of optical control of single spins, the optical control of two 
spins for implementing two-qubit quantum gates is an immediate milestone for future experiments. Indeed, the recent systematic 
investigation of the optical transitions in coupled QD structures [268. .315. ,317r.319J has laid a cornerstone for this target. 



V. QUBIT INITIALIZATION 



A rapid and continuous supply of refreshed qubits is one criteria for scalable quantum computation |2 |. Such a requirement is 
not only a prerequisite for the initialization of a quantum computer, but also a key element for quantum error correction where 
errors are continuously generated during the operation in noisy environments and by imperfect control. The initial preparation 
of a quantum computer could be done slowly, e.g., by simply cooling the system to very low temperature. For quantum error 
correction, however, rapid reset of qubits is crucial to recycle the spoiled qubits, otherwise a (infinitely) large number of fresh 
qubits should be prepared and preserved before a quantum computation commences so that the erroneous qubits could be 
replaced. The dynamical recycling strategy is more economical than a static supply of many qubits which deteriorate. It is 
desirable that the machine be as small as possible, with operation cost as trade-off'. 
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The essence of qubit initialization is preparing a pure quantum state out of a mixed one. It amounts to cooling a qubit (ideally) 
to absolutely zero temperature. Thus the key physical process is dumping entropy to the environment. The aim is to build a 
quantum refrigerator in analogy to a reverse Camot cycle. In general, the cycle consists of the following steps 11201 12561 12571 
[323 - 326 1 : (1) pumping of the system to an excited state; and (2) relaxation of the excited system with entropy dumped into 
the environment. To initialize the spin qubit in an ultrafast timescale, a quantum channel capable of dumping entropy rapidly is 
required. 

In a QD, the only available thermal baths for dumping the entropy of an electron spin are the nuclear spins, the host lattice (the 
phonon bath), and the electromagnetic environment (the photon bath). The coupling between the electron spin and the nuclear 
spins is very weak (with the rate in the order of 10"^ sec"^), so only the phonon or photon bath could be used as entropy drain. 
The direct coupling of an electron spin with either the lattice vibration or the electromagnetic modes is known to be very weak. 
The solution is to transfer the spin state into orbital states which couple to the phonon or photon bath strongly. The entropy 
dissipation by rapid photon channel can be realized by optically pumping the spin states to a trion state and coupling the trion to 
a photon in a strongly coupled QD-microcavity- waveguide structure 1 120|. Or alternatively, excited trion states may provide an 
efficient entropy channel realized by the rapid phonon emission in QDs. Below we discuss these two possibilities. 



A. Initialization by entropy dumping to photon baths 

1. Optical pumping 



The idea of initializing a spin by optical pumping is illustrated in Fig. [23j An optical light brings one of the two spin states 
into a trion state, and then the trion state relaxes to either spin state by spontaneous emission. After sufficient cycling of the 
pumping process, the electron spin will be in the state that is not coupled to the trion state by the pump light, which has been 
demonstrated recently in experiments for single electron spins in QDs I256II257I . The scheme has also been applied to initialize 
single heavy-hole spins in QDs |255 1. With essentially the same physics, optical pumping has also been used to initialize spins 
of NV centers in diamond 1901 l95l . 

To selectively excite only one electron spin state, one can choose the light polarization so that one transition is "forbidden" 
due to the selection rules. The selection rules are based on some symmetry, such as the (approximate) rotational symmetry about 
the growth direction. To exploit the selection rules resulting from the rotational symmetry, one can use a circularly polarized 



light. A (T+ polarized light, e.g., couples the spin state | + 1/2) to the trion state |^+) (see Fig. 23 ). Then the trion state will relax 
to both spin states by spontaneous emission. Now we have a dilemma: The selection rule prevents the spin state | - 1/2) to be 
reached. But remember the selection rules we have are approximate and the "forbidden" transitions are actually slightly allowed. 
According to the discussions in Appendix [a| the trion state \t+) has a small probability (~ f'^) to relax to the "forbidden" state 
I - 1/2). Also, due to the approximate selection rules, the state \t-) may be slightly excited, and the spin state | - 1/2) can be 
excited by the (T+ -polarized light to the two trion states, which degrades the initialization efficiency. Furthermore, a magnetic 
field can be applied along the growth direction (i.e., in the Faraday configuration), so that all the transitions except the desired 



ones as depicted in Fig. 23 are off'-resonance from the pump light. Thus the pump process is characterized by several rates: The 
excitation rate G, the spontaneous emission rates F and f = Ff'^ (to the spin up and spin down states, respectively), and the 
longitudinal spin relaxation rate ji = T~^. The rate equations for the trion population pt and the electron populations in the two 
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FIG. 23. Optical pumping of a spin in the Faraday configuration. The solid line indicates the optical excitation, and the dotted lines indicate 
the spontaneous emission. The polarization and the relative dipole matrix element for each transition are indicated. The two spins states are 
connected by a spin-flip process with rate yi . 
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spin states p+ are established as 



p+ = -yip+ +7iP-+ Tpt - Gp+ 
p- = -yip- -\-yip+ -^tpt. 



(17a) 
(17b) 
(17c) 



The quahty of the initiaUzation may be quantified by two factors: the saturation time Ts and the saturation spin polarization Pg. 
In a typical GaAs or InAs QD at low temperature and under a moderately strong magnetic field, F ~ 10^ s"\ yi ~ 10^ s"^ (see 
Sec. lie for details), and ~ 1% (for a unstrained dot |81 1, the value may be increased by strain tSOl ). Assuming the pumping 
rate G » F, we obtain by the rate equations (note the trion population will eventually becomes spin up population when the light 
is switched off) 



Ts ~ 10/f ~ 0.1 ms, 
Ps=P--p^-pt^l-2y,/Y-9^%. 



(18a) 
(18b) 



Such a high degree of electron spin polarization by optical pumping has been experimentally demonstrated 12561 . The initializa- 
tion by optical pumping in the Faraday configuration, however, is rather slow and the saturation polarization is limited by the spin 
flip rate relative to the "forbidden" spontaneous emission rate. Furthermore, the energy cost of the pump light is considerable. 
To see the energy cost of the optical pumping in the Faraday configuration, we notice that for just one useful photon emitted 
(which results in the target spin state), the number of photon wasted (by spontaneous emission resulting in the original state) is 
^'-2 ^ 400 - only 0.25% energy of the pump light has been eff'ective. The limiting factor is the small dipole moment for the 
"forbidden" transition. 

The solution is straightforward. One can work in the Voigt configuration in which a magnetic field is applied in plane to 
quantized the electron spins to be | ± x) 111811257113251 . As shown in Fig. 24 now both spin eigenstates are connected to a 
trion state by a (T+ -polarized excitation with the same dipole matrix element. The efficiency of optical pumping in the Voigt 
conffguration can be discussed similarly to that in the Faraday conffguration. But now there is another limiting factor of the 
saturation polarization: the backward optical excitation process from the target spin state to the trion state which, though off- 
resonance, is not negligible because of the large dipole matrix element. The generation rate of the trion by the off'-resonance 
excitation is G = GF^/ (f^ -h cofj where coc is the angular Lamor frequency of the electron spin in the magnetic ffeld. The rate 
equations for the trion population Pt and the electron populations in the two spin states P+ are established as 

Pt = -2Tpt + Gp++Gp., (19a) 
p+ = -yip+ +yiP-- Gp+ -h Tpt, (19b) 
p- = -yip- +yip+- Gp- + Ypt. (19c) 

lO^F ~ 10^^ s"^), the efliiciency of the initialization by 



Under the condition that G » F » G » yi (e.g., G = lOF and cDc = 
optical pumping in the Faraday conffguration is characterized by 

Ts ~ 10/F ~ 10 ns, 
Ps=p--p+^l-G/T 



99.9%. 



(20a) 
(20b) 



Owing to using the allowed transition in the Voigt conffguration instead of the forbidden transition in the Faraday conffguration, 
the spin initialization is faster by orders of magnitude. Also, the saturating polarization is much closer to unity since now the 
limiting factor is the off'-resonance transition probability relative to the resonant one, instead of the spin-ffip rate relative to 




FIG. 24. The same as Fig. [23] except that the setup is in the Voigt configuration. 
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the trion recombination rate due to the forbidden transition in the Faraday configuration. The off'-resonance transition could be 
suppressed simply by enlarging the electron Zeeman splitting. There are calculations 13251 and experiments [257,1 demonstrating 
efficient optical pumping of single electron spins in the Voigt configuration. 



2. Single-shot initialization with cavity enhancement 



In the previous discussion on spin initialization by optical pumping, we have seen that it is crucial to have a rapid entropy 
dumping channel to have efficient qubit cooling. The cooling duration of 10 ns is acceptable in many cases, but it is still highly 
desirable to have even faster spin initialization so that the qubit refreshing rate could catch up with the quantum gate and the 
error generation rate. Suppose the error rate per qubit per operation is about the quantum error correction threshold 10"^ and 
about 10 physical qubits are used to code one logical qubit to incorporate the error correction within the quantum logic. As we 
have discussed, the optical control for a simple quantum gate should take about 10 ps. That means we should have a qubit reset 
rate well above 1 qubit per nanosecond to avoid quantum computation being held up due to the lack of refreshed qubits. To 
have such ultrafast spin initialization, the entropy dumping channel should be specially engineered. One possible approach is to 
enhance the coupling to the photon bath by increasing the local density of states of photon modes. This enhancement is possible 
by putting a QD in the proximity of a microcavity 1 120|. As has been discussed earlier, a micro-photonic structure is in any case 
needed to form a scalable large structure of QDs for distributed quantum computation. Here we discuss how the in-situ cavity 
QED may be used to selectively enhance the photon emission for entropy dumping and hence for ultrafast spin initialization, 
following the procedure of Ref. 1 120|. 



The basic structure is depicted in Fig. 25 To be specific, we consider a micro-ring coupled to a QD. Such a structure may 
be constructed by etching a semiconductor surface where QDs are located. Structures in photonic crystals should be ideal 
alternatives (see Sec. |IIIB| ). The specific structure, however, is not crucial to the physics to be discussed below. The attached 



micro-cavity would strongly modify the electromagnetic vacuum in the vicinity of the QD. The coupling to the cavity mode 
(which is taken as a whispering gallery mode) is realized due to the overlap between the QD and the evanescent wave of the 
cavity mode. A waveguide coupled to the cavity serves as a quantum channel for a cavity photon to escape rapidly to the 
environment. For the purpose of spin cooling, such a directional waveguide is not necessary, and actually one could as well use 
a cavity with large leakage (bad cavity situation). Since the emitted photon carries certain information about the qubit, a guided 
channel with cavity also enables such information to be retrieved, either for quantum measurement (as will discussed later) and 
for quantum error diagnosis. The incorporation of a waveguide along a cavity then involves extra designing and fabricating cost. 

We consider the Voigt configuration. The spin eigenstates under a static magnetic field in the x direction are denoted |±). The 
spin states may be flipped to the two degenerate trion (exciton plus an electron) states |r=F) or |r±) by an X- or F-polarized 
tipping pulse, respectively. The trion states are, by design, ofl'-resonance from the cavity modes. We assume such ofl'-resonance 
condition for several considerations: (1) It avoids the cavity-induced optical decoherence during quantum operations of the spin. 
(2) A local node usually consists of a number of QDs to a cavity, so it is unlikely to have all the QDs are in resonance with a 
cavity mode. (3) The resonance coupling may be realized by transient control via optical Stark efl'ect which provides flexibility 
for selectively initialize a spin in a QD in the cluster. The relevant cavity mode is denoted by |C). The evanescent wave of the 
cavity mode is designed to be X-polarized in the vicinity of the nanodot, so that when brought within resonance, the trion states 
|r±) and the cavity states |=F, C) are coupled into two split trion-polariton states, respectively. This provides a fast decay of the 
trion to a spin state by emitting a photon into the quantum channel. The pump light is 7-polarized so that a Raman pathway is 
formed from the spin state |+) to the trion state \T-\-) (by the optical pumping), to the cavity state |-, C) (by dot-cavity coupling), 
and to the spin state |-) ( by spontaneous photon emission into the waveguide). 



The optical cycle (similar to a reverse Carnot cycle) for cooling a spin qubit is illustrated in Fig. 26 1 120|. Without loss of 



generality, we assume that the electron before optical pumping is in an unpolarized state, i.e., p(-oo) = 0.5|-)(-| -\- 0.5|-h)(-h|. A 




FIG. 25. Schematics of the dot-cavity- waveguide coupling structure. 



32 



cooling cycle consists of four basic steps: 

1. An X-polarized AC Stark pulse is adiabatically switched on, bringing the states and |-,C) into resonance by AC 
Stark effect. 

2. A F-polarized pump pulse flips the spin up state |+) to the polariton states formed by and |-, C). 

3. The polariton states relax to the spin down state |-) rapidly by emitting a photon into the waveguide, dumping the spin 
entropy to the environment. 

4. The AC Stark pulse is adiabatically switched off*. No photon-generation or spin-flip would take place if the initial spin 
state is |-). 

Ideally, after one cooling cycle, the spin is fully polarized with the entropy mapped into the quantum channel, and the final 
density matrix becomes |-)<-| (0.5|0)(0| -h 0.5|1)(1|) , where \n) is the ^-photon waveguide state. This is an idealization of the 
single-shot initialization of a spin qubit in a QD. In reality, the single-shot initialization is subject to errors due to the spontaneous 
emission of photons into free space, by which the trion state may relax to either spin state depending on the polarization of the 
emitted free-space photon, while the cavity photon couples only the transition to the target spin state. 
The whole system under the optical control is described by the Hamiltonian 

H = ^ca^a ± y l±><±l + ^t\T±){T ± \ + {g,^,\T±){^a + H.c) 

+ [xt{t)^t ^Xpit^p] • ex (|7^±)<=F| + rca^) + H.c 

+ [xtm +Xpit)^p] ' er|r±)<±| + H.c, (21) 

where a is the cavity mode annihilation operator, Qc is the cavity mode frequency, is the bare trion state frequency, ojc is the 
electron Zeeman splitting, gcav is the coupling between the cavity and the QD, Xp is the Rabi frequency of AC Stark pulses with 
polarizations e^, Xt is the Rabi frequency of the pump pulse with polarization e^, and rc is the strength of the direct coupling 
between the light pulse and the cavity mode (relative to the coupling to the trion state). 

The dynamics in the cavity-dot system is rather complicated and some attention should be paid to the designing of the 
controlling pulses. To bring the trion state into the resonance with the cavity mode and to maintain the resonance, the AC Stark 
pulse is designed to have an almost square profile. Also, the switch-on and off' of the pulse should be made smooth enough to 
avoid non-adiabatic excitation of the trion states from the target spin state. We choose the pulse to be of the profile 

Xp(t) = Xpe-'""^' [^rf{cTp(t -h))- erf [cTpit - ^2))] , (22) 



[see Fig. [28] (a)]. As the AC Stark pulse maintains the resonant cavity-dot coupling which facilitates the photon escape to the 
quantum channel, the trion state relaxes very fast (on the time-scale of g~^^ and y"^ ~ 10 ps). The ffipping pulse should be 
a TT-pulse for Rabi rotation between and \T-\-}. Due to the dynamical nature of the states (dressed by the AC Stark pulse) 
and the rather small polariton splitting (~ 0.1 meV), a perfect 7r-rotation requires an extremely long pulse. In principle, a 
full excitation of the spin up state to the polariton states could be made much faster by pulse shaping. One of such pulse 
shaping is to use geometrical control of the transition which is robust against uncertainty of the polariton frequencies. The 
geometrical control is realized by using a a chirped pulse as XtiO = ;^?^"^'^^^^"'^'^sech (crf(t - tt)) with the phase sweeping rate 
(pit) = -(Tc tanh ((T^(^ - tt)) 13141 . The frequency of the pulse now will sweep from above to (T^ below and the sweeping 




FIG. 26. Basic optical processes for initializing a spin by controlled cavity QED. The dotted, solid, and wavy arrows represent the AC Stark 
pulse (X-polarized), the tipping pulse (7-polarized), and the spontaneous emission, respectively. 
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range [Q^ - cr^^t + crj covers both of the trion-polariton states. The initial spin state |+) will be brought adiabatically into 
a superposition of the two polariton states, which relaxes rapidly to the target spin state |-). Such a geometrical flip can also 
tolerate to some degree laser fluctuations and uncertainty in the dipole moment, transition energy, and selection rules. 
The cooling process is simulated by numerically solving the master equation of the dot-cavity system 

dtp = -i [fi.p] - ^-^^LaP - ^ X -^I^X^^'I^' 

where LbP = 2dpd^ - o^op - pb^b denotes a Lindblad super-operator, y is the cavity- waveguide escape rate, y' is the cavity- 
free-space loss rate, F is the trion decay rate due to spontaneous emission into free-space. The multi-photon cavity states 
were included in the numerical calculation, as they renormalize the AC Stark shift (the real excitation of multi-photon states 
is negligible due to the ofl'-resonance condition). Inclusion of up to 3 -photon states was found sufl&cient to obtain converged 
results. 

We test the cooling efliiciency with a set of realistic parameters fT20'l. The Zeeman splitting coc - ^ meV, y - 0.2 meV, 
y' - 0.045 yueV (corresponding to an intrinsic Q-factor ~ 3 x 10^), the dot-cavity coupling gcav = 0.1 meV, the cavity-trion 
detuning Qc - - = 0.5 meV, F = 1 yueV, and tq - 0.3. For the parameters given above, the trion state and the 
cavity state |-, C) are brought into resonance when the AC Stark pulse strength i^Xp) is maintained at 1.21 meV. Maintaining 
the resonance for h-h - 70 ps is found sufficient for the total dissipation of the photon. The spectral width of the AC Stark 
shift pulse {p-p - 0.354 meV) is set much smaller than the detuning (Q^ + ^l/2 - = 5.5 meV), so that the excitation due to 
non-adiabatic switch-on and off' is negligible. The ffipping pump pulse has a frequency sweep range of cTc - 0.4 meV, strength 
Xt - 0.2 meV, and duration l/cr^ = 6.58 ps. The spin state is ffipped to the polariton states with negligible error. 

Figure [28] (b) shows that a single cooling cycle completed within 80 ps produces an almost 100% polarized spin from a 
maximally mixed state. The density matrix at the end of the cycle is p = 0.9945|-)(-| -h 0.0040|-h)<-h| -h perr, where pen- is the 
probability (^^0.15%) of the system remaining in the trion states which results mainly from the non-adiabatic switching of the 
AC Stark pulse. The extra error (^^ 0.4%) comes mainly from the decay of the trion with photon emission into free space. 



B. Initialization by entropy dump to phonon baths 

The phonon bath in a QD is often taken as a source of the qubit decoherence, but it can also be used as as a resource for ultrafast 
cooling. The electronic energy levels in QDs are discrete, it has been argued that the phonon emission would be much suppressed 
due to the lack of available ffnal states fulfflling the energy conservation | 327 1. This so-called phonon bottleneck eff'ect, because 
of its importance in QD lasers and detectors, has been extensively studied both in experiments and in theories. Nonetheless, 
experiments have established various mechanisms for rapid relaxation of electrons from exited states to the ground states, such 
as the Auger-process and the multi-phonon process 13281 . The observed relaxation time varies from tens of picoseconds to a 
few picoseconds 1 3291 13301 . The spin relaxation, especially for holes due to the large spin-orbit coupling, always accompanies 
with the energy relaxation, for the mixing in the excited states is much stronger than in the ground states. Thus we can use the 
phonon emission as fast entropy dumping channel for an electron spin in a QD. 



A cycle of the spin initialization is illustrated in the inset of Fig. 29 



1. A circularly polarized laser pulse resonantly pumps the electron states to an excited trion state. 

2. The excited trion state relaxes rapidly to the ground trion states | tvlH) and | Tvlli)- During the relaxation, the hole spin is 
not conserved due to the strong heavy-light hole mixing in the excited state, and without loss of generality we assume the 
rate is the same for the relaxation to the two ground trion states. 



(meV) 


6 - 




4 - 








2 - 


xl 




0) 

c 


- 




CD 


-2 - 


1 1 




1 2 1 

pump strength (meV) 



2 



FIG. 27. (Extracted from Ref. 1 120 |) A quantum version of the reverse Carnot cycle for spin qubit initialization in a QD. The grey curves are 
the energies of different states versus the Rabi frequency of the AC Stark pulse, in the rotating frame. 
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3. To deplete the population of one of the ground trion states, a circularly polarized laser pulse is applied to resonantly couple, 
e.g., the trion state | tvlli) to the electron spin state | t)- 

4. After a period of pumping and relaxation, the population will be accumulated to the trion state | TvlH). and, thus, an 
ultrafast 7r-pulse can be applied to flip the trion state into a pure electron spin state. 

The most time-consuming step of this laser cooling process is the pump and relaxation processes which is limited by the carrier 
relaxation rate and the pump strength. Care should be taken to avoid the population being trapped in some dark state which could 
be possible in the strongly driven multi-level system (similar to the electromagnetically induced transparency efl'ect). Fig. [29] 
shows a typical spin initialization cycle in which the spin is pumped from a fully mixed state to a nearly fully polarized state 
in 120 ps (the polarization at the end of the cycle is about 97%), where p^it) is the population stored in the ground trion state 
I TvlH). and P(t) is the polarization of the electron spin. The pump strength for the excited and ground trion states has been set to 
be 2 meV and 0.5 meV, respectively, the carrier relaxation time is 2 ps, and the spontaneous emission time is 1 ns. The 7r-pulse 
flipping the trion state into the spin state is chirped for optimized performance. 

The phonon baths ofl'er an alternative solution of ultrafast spin initialization with comparable speed to a photon bath enhanced 
by cavity QED. The real excitation of the excited electron and hole states, however, could cause some complications. One such 
issue is the trapping of electrons and holes in the dark states, which cause the loss of a spin qubit. And in fluctuation QDs, it is 
not likely the excited state energy is high enough for LO phonon emission. 



VI. QUANTUM NON-DEMOLITION MEASUREMENT OF SPIN QUBITS 

The result of a quantum algorithm is obtained by quantum measurement of specific qubits. To diagnose the errors generated 
in quantum error correction, quantum measurement is also needed. 

In current single-spin experiments in QDs I79l 12561 f273l 13261 133 1H333II and in diamond color centers | 88J, the measurement 
is usually achieved by cycling read of the spins, which amounts to time-ensemble measurement. Remarkably, eflficiency en- 
hancement in measurement of NV center spins in diamond has been made using ancillary nuclear spin qubits I266II334I , but the 
single-shot measurement still remains illusive. 

Eventually, single-shot measurement of quantum registers is demanded for scalable quantum computation for large-scale 
problems. In Appendix [Cj we show that quantum non-demolition (QND) measurement, even with certain errors, may be con- 
verted to a single-shot measurement and thus is scalable. To realize single-shot measurement, the crucial issue is how to enhance 
the coupling between the probe (photons, e.g.) and the spin qubits. Below we provide a possible solution of using cavity QED, 
which may be implemented in-situ in a quantum network (discussed in Sec |VII| ). The discussion is based on QD cavity-QED 
systems, but may be readily extended to NV centers where cavity-QED have been also demonstrated 1.260 J . 
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FIG. 28. (Extracted from Ref. 11201 ) (a) The Rabi frequencies of the AC Stark pulse and the tipping pulse (amplified by a factor 5), and the 
sweeping frequency of the tipping pulse, (b) Probabilities of spin down and up. Different steps of the cooling cycle, indicated by (D-®, are 
distinguished by shadowed areas. 
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A. Scalability of quantum measurement in quantum computation 

In quantum computation (such as in Shor's algorithm IH 11721 ). the quantum register could be in a superposition state right 
before the measurement |^) = 2x ^xk), where \x} = |xi, X2, . . . , xn) is a computational basis state for the A/^-qubit register. The 
measurement should be in the computational basis which returns any \x) contained in the superposition and the computation 
result is derived from the measured x. If an algorithm requires the measurement of the wavefunction Cx, i.e., the tomography 
of the quantum state, it would be an analog computation instead of a digital one, and much worse, it would not be scalable as 
the number of measurements would increase exponentially with the number of qubits in the register to be measured L335 J . Such 
measurements are ensemble measurements. 

A point we would like to put forward here about the scalability of quantum measurement in quantum computation is that 
an ensemble measurement is not scalable in the sense that the size of the ensemble would increases exponentially with the 
problem size (defined as the number of qubits registering the computation result involving a quantum measurement) I12QI . 
It has been well-known that an ensemble quantum computation is not scalable if the qubits cannot be initialized to a pure 
state 1336 1. The issue of scalability associated with the ensemble measurement |120|, however, has received less attention in 
spite of its importance in quantum computation. This problem is briefly explained below, and more detailed discussion is given 
in Appendix [C] 

First we notice that the uncorrected ensemble measurement cannot be used to read out the quantum register in general 
algorithms (especially for those terminating in superposition states such as Shor's algorithm). In an uncorrelated measurement, 
the spins are measured independently. Thus, in general which basis states are in the superposition cannot be deduced from the 
measured result. For example, the superposition state |000) + |111) will give the same uncorrelated measurement result as the 
state |000) + 101 1) + 1101) + 11 10). In some algorithms such as Shor's algorithm for factorization, the number of basis states in the 
superposition may increase exponentially with the number of qubits, thus the number of possible superposition states yielding 
the same uncorrelated measurement result would increase exponentially with the number of qubits measured. So the quantum 
measurement has to be a correlated one (e.g., the photon counting should be in coincidence) to be scalable. If a coincidence 
measurement is a destructive one, the procedure has to be run from the very beginning in each repetition, and the superposition 
state can collapse into any possible which is in general difl'erent from one cycle to another. To have a certain x to be measured 
at least twice for the sake of confidence, the number of repetitions to be performed should be in the order of the number of basis 
states in the superposition, which again could be an exponential function of the problem size. In conclusion, the enhancement of 
signal-background contrast by ensemble quantum computation is not a scalable solution. 

In realistic cases, the signal of a single- shot measurement of a single quantum object is usually too weak to be distinguished 
from noise. Thus the signals are to be amplified either by simultaneously measuring a large number of identical "quantum 
computers" running the same quantum program or by repeating the quantum program under identical conditions for a large 
number of times. This renders the quantum computation to be an ensemble one and hence not scalable. 

Using ensemble measurement as a solution to the detection efficiency problem has been applied in various systems, including 
nuclear spins in liquid-phase NMR |336| or solid-state NMR 1621 . Here we consider quantum computing with optically con- 
trolled spins. The spin-dependent absorption (or other optical methods such as spin-dependent scattering and Faraday rotation) 
may be used to measure a single spin in a QD, in which a probe pulse resonant with the spin-trion transition in a QD measures 
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FIG. 29. The numerical simulation of a typical initialization cycle using phonon emission from excited trion states. Inset: Schematics for the 
spin initialization process. 
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the spin in the computational basis by detecting whether or not a single photon of the probe pulse has been absorbed. To measure 
a register with many qubits, the probe should be composed of many pulses in which each pulse addresses a QD individually 
by spatial and spectral resolution. In practice, the interaction between the probe pulse and a single spin is very weak, and also 
photon collection and detection efficiency is less than one, so the information obtained about a spin is on the average much less 
than one bit. The probe measurement has to be repeated for many times to accumulate statistical confidence in the measured 
result (or signal-to-noise ration). Considering the fact that the probe process is destructive to the spin state (since a trion state 
excited by the probe pulse may return to either spin state regardless of the original spin state), the repetition has to be run from 
the very beginning of the quantum computation which prepares the quantum computer to the same superposition state. This 
makes the quantum computation an ensemble one, which in general is not scalable. 

Nonetheless, the ensemble measurement could still be useful in demonstrating quantum algorithms for small-size problems. 
For instance, as will be shown later, to demonstrate Shor's algorithm for factorizing 15 by optical control of electron spins in 
QDs, the program would be completed in less than a few ns, and the number of computational basis states in the superposition is 
less than 10. A commercial Ti-Sapphire pulse laser with repetition rate of about 100 MHz could be used to carry out an ensemble 
of repeated running and measurement in a reasonably short time in a pump-probe configuration, where the initialization and gate 
control are viewed as a single complex pump pulse and the probe pulse is composed of many frequency components (and 
detected in multi-channels.) 



B. Quantum non-demolition measurement via cavity quantum electrodynamics 



In essence, the realization of an efficient single- shot quantum measurement involves two key elements: rapid quantum state 
entanglement of the target qubit (here a spin in a QD) with a detectable information carrier (such as a photon), and efficient 
and faithful collection of the information carrier (such as by a photon detector with a high efficiency and a low dark count rate). 
Errors in the measurement, due to the efficiency or dark count problems, e.g., can not be fully eliminated. Thus cycling of 
single-shot measurement is required to accumulate statistical confidence. As discussed in Appendix [Cj the cycling of single- 
shot measurement is scalable when the measurement is a quantum non-demolition (QND) one, i.e., the qubit state after the 
measurement is (ideally) an eigenstate in the measurement basis corresponding to the measurement output. 

The QND measurement being rapid is an essential element in the following sense: The cycling of measurement has to be 
completed in a time much shorter than the qubit is significantly disturbed by the environment. So we should look for a rapid 
quantum information transfer between the spin qubit and a medium to be detected. An ideal medium is photons. The cavity- 
enhanced entropy dumping in the ultrafast initialization is indeed a rapid information transfer (but there the information has 
been viewed as noise). Thus the ultrafast spin initialization and rapid QND may be integrated in the same micro-photonic 
structure |120|. 

To ensure efficient detection of the transferred quantum information, the photons escaped from the cavity should be directed 
into a quantum channel, unlike in the initialization process where it does not matter which direction the dumped entropy flows. 
To realize such a directional information flow, a waveguide could be fabricated in the proximity of the cavity, and the waveguide 
may be coupled to a fiber which conducts the photon to a detector I295H301I . 

The optical control of the cavity-QED for a rapid measurement is similar to that for the ultrafast initialization 1 120|. But to 
enable measurement cycling, the measurement should be non-destructive, i.e., the spin basis state should return back to its initial 
state after a cycle of measurement, with close to unity probability. Thus we need to switch the polarizations of the tipping and 
the AC Stark pulses from (7, X) to (X, 7), respectively. The energy diagram with optical transitions indicated is shown in Fig. 30 
A measurement cycle may be processed in four basic steps (see Fig. [31]): 
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FIG. 30. Basic optical processes for QND measurement of a spin by controlled cavity QED. The dotted, solid, and wavy arrows represent the 
AC Stark pulse (F-polarized), the tipping pulse (X-polarized), and the spontaneous emission, respectively. 
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1. An X-polarized tipping pulse flips the spin state |+) to the trion state 

2. A F-polarized AC Stark pulse adiabatically switched on drives the trion state into resonance with the cavity state |+, C). 

3. The trion state resonantly tunnels into the cavity state and relaxes rapidly back to the spin state |+), leaving a photon 
emitted into the quantum channel. 

4. The AC Stark pulse is adiabatically switched off. 

Suppose that the spin state to be measured is Qf|+)+/5|-) and the channel is initially in the vacuum state |0). An ideal measurement 
process will transform the system into the entangled state Qf|+)|l)+yS|-)|0), so that the detection of the photon projects the electron 
into a spin basis state, providing a QND measurement of the spin. Note that the pulse timing for measurement is difl'erent from 



that for cooling [cf. Fig [28] (a) and Fig. [32] (a)]. In the measurement cycle, the flipping pulse need not be chirped, since here 
the Rabi flop occurs between stationary energy levels and the transition between the spin state and the trion is well separated in 
frequency from the cavity mode. Instead, a simple Gaussian n pulse ;^^(0 = Xte~^'^^~^'^ l^-i^tt used. The AC Stark pulse 

is chosen F-polarized to avoid direct excitation of the cavity mode. 



Numerically simulation is done to check the efficiency of the measurement cycle 11201 . Since the initialization and the 
QND measurement are supposed to be operated via the same photonic structure, the simulation is done with the same physical 
structure as that in Fig. [28] The number of photons emitted into the waveguide is calculated with dtU = y{a^a). The tipping and 
the AC Stark pulses are set such that l/cr^ = 2.19 ps, Xt = 0.192 meV, = - cTp = 0.707 meV, 2xp = 2.08 meV, 

Qt-^^l/'^-^p = 5.5 meV, and the duration of the pump pulse t2-ti = 50 ps. After a single cycle of measurement, an initial state 
Po = results in the final state pi = 0.0161 + 0.9824|+)(+| +perr with the number of photon emitted into the waveguide 



n = 0.9806 [see Fig.|32](b)], while an initial state po = |-)<-| results in the final state pi = 0.9955|-)(-| + 0.0040|+)(+| + pen- 
with n = 0.0015 (not shown). The photon emitted into the waveguide can be detected with high efficiency. If the detector has a 
zero dark-count rate and an efficiency of 50% [337 J , the POVM (positive operator- valued measures II172L see Appendix [C] for 
more discussions) for the measurement process can be defined as 

P. = 0.9992|-)<-| -h 0.5097|-h)<-h|, (24a) 
P+ = = 0.0008|-)<-| 0.4903|-h)(-h|, (24b) 

for a non-click or click event, respectively. Within 5 measurement cycles, e.g, the spin state is measured with accuracy higher 
than 97%, and the back-action noise to the spin is less than 10%, while the time duration is less than 0.4 ns, much shorter than 
the spin decoherence time. 

The QND measurement of single spins thus can be completed within 100 ps, and the high-efficiency of collecting and detecting 
photons propagating in waveguides enables a near unity accuracy by repeating the single- shot measurement for only a few 
times. The measurement is still subject to the problem of less than unity efficiency as well as back-action noise. As discussed 
in Appendix [C] quantum gates and error tolerating coding can be combined to achieve a sufficiently faithful measurement of a 
qubit without rewinding a quantum computing program. 



VII. NETWORKING LOCAL NODES 



In section IIIB we have introduced the dot-cavity- waveguide coupled structure for distributed quantum information process- 
ing. Quantum nodes are formed by clusters of singly charged QDs with electron spins as carriers of stationery qubits. Single 
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FIG. 31. (Extracted from Ref 1 120 |) Detailed optical process for the measurement cycle. The grey curves are the energies of different states 
versus the Rabi frequency of the AC Stark pulse, in the rotating frame. 
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photons in optical waveguides or fibers can transport quantum information between distant nodes. Interfacing between single 
spins in QDs and single photons is made possible by the strong photon confinement in solid state micro-cavities. 

The separation between quantum nodes on a single chip which can range from ~ jim to ~ cm allows parallel optical control, 
and intra-chip communication is realized by optical waveguides. Inter-chip communication is possible by wiring chips together 
with optical fibers 1338 -3401 . For intra-chip communication and short distance inter-chip communication, the decoherence of 
the photon qubit is negligible I341H343I . Thus, it is possible and also highly desirable to perform inter-node operations in a 
deterministic way (to be contrasted with most quantum cryptography and linear optics quantum computation schemes based on 
projective measurement which renders the logical controls probabilistic). The key component is a quantum interface that allows 
the deterministic state transfer between spin and photon qubit. 

The prototype quantum interface for this purpose was proposed by Cirac et al fTJT] . It is composed of a cavity coupled to a 



three-level A system, illustrated in Fig. 33 The two ground states, \g} and \e}, of the three level system form the stationary qubit. 
State 1^) is coupled to the intermediate state \t) by the cavity mode and \e) to \t) by the external laser field. Direct excitation 
of cavity by the external laser is assumed absent. Through the imperfect mirror, the cavity is coupled to the electromagnetic 
continuum which forms a photonic channel. A Raman path from |^) to |g) through the intermediate state \t} is thus formed. If 
the three level system is initially in state \e), an external laser pulse can bring it to state \t) by a tt rotation which relaxes to state 
\g) by spontaneous emission of a cavity photon. The cavity photon then goes into the photonic channel forming a single photon 
wavepacket. If the three level system is initially in state \g}, it will remain in this state provided the cavity is in its vacuum. The 
quantum state carried by the three-level system is thus mapped into the photon number subspace of the outgoing photon, 

{Cg\g) + Ce\e)) Ivac) ^ \g) [c,|vac) + CJ^out)] • (25) 

where \a) denotes a single photon wavepacket in the photonic channel and |vac) the channel vacuum. This process forms the 
basis for the sending function (i.e. the mapping from a stationery qubit to a flying qubit) of a quantum node. The receiving 
function is the mapping from a flying qubit to a stationary qubit and can be realized as the time reversal of a sending process, 

\g} (C,|vac) + C,|a,n)) ^ (Cglg) + Ce\e}) |vac). (26) 



With the output of the sending node directed as the input of the receiving node (see Fig. [33]), transfer of qubits between two 
distant nodes can be performed. 

This quantum network by cavity QED was originally proposed for quantum computation with atomic systems II277II . Critical 
experimental steps towards realizing such a quantum interface in atom-cavity QED systems have been progressively demon- 
strated 1 153, 154, 344 1, including the initial demonstration of reversible state transfer between photons and atoms 13451 . Via 
similar cavity-assisted Raman processes, schemes for mapping between motional states of single trapped atoms 13461 or collec- 
tive excitation of atomic ensembles li347J and the quantum states of single photons are also proposed. 




FIG. 32. (Extracted from Ref. II120I ) (a) The Rabi frequencies of the AC Stark pulse and the tipping pulse (amplified by a factor 10). (b) 
Probabilities of spin down and up, and the number of waveguide photons, for a spin initially polarized up. Different steps of the cooling cycle, 
indicated by ®-(4), are distinguished by shadowed areas. 
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For quantum computation with QD spins (63], deterministic quantum network control is indeed possible in the distributed 
structure discussed previously in section IIIB| In subsection VII A we will first show how to realize such a prototype quantum 
interface between a QD electron spin and a photon mode in a waveguide/fiber. In subsection ] VII B| we will describe the exact 
solution to the interface dynamics between a spin qubit and a photon qubit in the prototype quantum node in the most general 
scenario. Control schemes based on the exact solution form the basis of a variety of inter-node operations in quantum networks 
as discussed in Sec. VII C The issue of unavoidable inhomogeneity of solid-state quantum nodes is properly resolved with the 
exact interface solution. In subsection |VII D| we study the eff'ects of various sources of errors on the interface operations. A 



summary and outlook are given in subsection 



vira 



A. Dot- Cavity- Waveguide structure as spin-photon interface 



Here we show how to implement a dot-cavity-waveguide coupled structure as a deterministic quantum interface for single 
spins in QDs and single photons in waveguides. 

The coupled structure can be realized with any of the micro-cavity systems discussed in section [niB] In Fig.[34ja), a QD 
sitting in the evanescence field of whispering gallery mode of a ring cavity is illustrated as an example. There are two essential 
requirements for such a coupled structure to be an efficient interface. First, the dot-cavity coupling must have a large Purcell 
factor, so that the QD optical transitions are dominantly coupled to the cavity field. Second, the leakage of the cavity photon 
into free space should be small as compared with the tunneling into the attached waveguide. If both conditions are fulfilled, the 
evolution through the desired quantum pathway, i.e. QD trion cavity photon waveguide photon, occurs on a much faster 
timescale than the leakage out of it 111911120 '!. Lowering of the cavity Q-factor due to the coupling to the waveguide is part 
of the process and has no deleterious eff'ects on the quantum operation. The cavity with the reduced 2-factor is not required to 
be in the strong coupling regime II 191 13481 13491 . This will become clear when we discuss the control of this interface in the 
following subsections. 

The qubit is represented by the two spin states \x-) and \x-\-) and spin manipulations are mediated by the two trion ground 
states and \T-\-} (see section II B ). In the convention of the prototype quantum interface, the two spin states are also denoted 
here as \g) and \e) respectively, and the two trion states \T-\-) and as \t) and \t), respectively. While \g) and \e) have energies 
ojg and oje in a static magnetic field normal to the optical axis of the dot (see Fig. 34 a)), \t) and \t) typically have a much smaller 
energy splitting (ojt ^ ojf) in GaAs fluctuation dot because of its negligible in-plane g-factor of the heavy holes I2QQII35QL but a 
splitting comparable to the electron's in self-assembled InAs dot | 257 |. 

By design, the selected cavity mode of frequency ojc is X-polarized in the vicinity of the QD and a 7-polarized control laser 
of central frequency oji and complex Rabi frequency Q(0 directly couples to the QD transitions I23QI 13081 . Therefore, by 
the selection rules shown in Fig. 34 the cavity field couples only to the straight transitions \g} \t} and \e} \t) , and the 
controlling laser couples only to the cross transitions \g) \t) and \e) \t). The laser light and cavity mode satisfy the two- 
photon resonance condition: <jJl-\- oje = coc -\- oJg. The Raman detuning A = ojt - ojl - oJe is also much smaller than the electron 
Zeeman splitting oje - oJg. Thus, by the Zeeman splitting and the selection rules, the trion state \t) is off-resonance to the laser 
light and the cavity mode [shown by dashed lines in the Fig. 34 c)]. 

At a sending node, the Raman process consists in first the laser field exciting the spin state \e, 0) to the trion state |r, 0), then 
the trion state resonantly coupled to the cavity state \g, 1) which finally is rotated to the spin state \g, 0) forming a photon wave 
packet in the waveguide (here and 1 denote the number of photons in the single cavity mode). The receiving mode is just the 
time-reversed process. Undesirable dynamics involving the state \t) is eliminated by making the Zeeman splitting sufl&ciently 
larger than the cavity-dot coupling and the Rabi frequency. The resultant optical process is the cavity-assisted Raman process in 
a A-type three-level system as required by the prototype quantum interface (see Fig. [33}. Numerical calculations including the 




FIG. 33. Illustration of state transfer in a quantum network. The node is composed of a cavity coupled to a three-level A system. The two 
ground states \g) and \e) of the three level systems form the Hilbert space for the stationary qubits. State \g) is coupled to the intermediate \t) by 
the cavity mode with strength gcav and \e) to \t) by a classical light with Rabi frequency Q(t). Direct excitation of cavities by the classical light 
is assumed absent. The cavities themselves are coupled to the outside continuum which forms a photonic channel. Two nodes are connected 
by the photonic channel in the following way: the output of node 1 is directed to node 2 as its input and vice versa. 
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non-resonance transitions and realistic decoherence have been performed and high fideUty of desired operations at the quantum 
interface is demonstrated (see subsection ] VII D| for details). 



B. Control of spin-photon interface 

In the previous subsection, we have shown that the dot-cavity-waveguide coupled structure could be an efficient solid state 
realization of the prototype quantum interface. In order for such quantum interface to be suitable for a quantum network, proper 
control schemes are required. 

The difficulty in realizing the network lies in the receiving end. Instead of being trapped at the quantum node, the single photon 
pulse can be reflected by the cavity unless the pulse shape of the classic control laser matches the single photon wavepacket 
exactly, known as quantum impedance matching |347, 351 1. One way to deal with this requirement was provided in 12771 . The 
central idea is that if a laser pulse can be found for generation of an outgoing photon wavepacket with time reversal symmetry, by 
using this laser control at the sending node and its time reversal at the receiving node, the time reversal symmetry will guarantee 
the photon wavepacket to be completely trapped at the receiving node. A solution for such a laser control pulse was given 
inlETTl. 

The time reversal symmetric control scheme requires the sending and receiving quantum nodes to be identical in terms of 
optical transition frequencies and strength of coupling between the components. Unfortunately, such requirement is practically 
impossible to fulfill in solid state systems since the fabricated QDs and micro-cavities naturally have shape variations and size 
fluctuations, leading to inhomogeneity in optical frequencies. The cavity field has a highly non-uniform profile and thus cavity- 
dot coupling diff'er from dot to dot. For realization of a quantum network with solid state structures, control schemes beyond the 
time symmetric one is required. In the sending function, it is straightforward to solve for the outgoing single photon wavepacket 
a(t) if we know the form of the laser control pulse 0^(t). The key problem is the inverse functional relation Q(0 = 9^~^[(^(t)], 
i.e., given an arbitrary single photon wavepacket a(t), the exact laser pulse that can generate this single photon wavepacket at a 
sending node, or completely trap it at a receiving node. Knowing this relation allows the sending and receiving operations to be 
separately addressed and hence the construction of a quantum network with heterogeneous quantum nodes. 

Solutions to D(0 = [(^(t)] were first given in the adiabatic approximations 1347113511 . Exact solutions in the non-adiabatic 
form were later found by the authors I119ll348ll , as briefly described below. 

The Hamiltonian including the interaction between the single-mode cavity and the three-level system and the channel contin- 




FIG. 34. (a) A high- 2 micro-ring coupling a 'tapered' waveguide and a singly charged QD. (b) Optical selection rules in the basis where spins 
are quantize along the field (x) direction. The hollow two-headed arrow denotes X-polarized light and the solid arrow denotes F-polarized 
light, (c) The level diagram and optical process. In \s,n), s = g,ej ox t denotes an electronic state in the dot and n denotes the number 
of photons in the single cavity mode. Straight, curved, and wavy arrows represent the laser excitation, dot-cavity coupling, and cavity-fiber 
tunneling, respectively. The resonant and off'-resonance processes are represented by solid and dashed lines, respectively. 
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uum IS, 

H — (jJcCi^ a + cOf \t) {t\ + (jJe \e) {e\ + I dojojbl^bt^ + gcav (i 10 (^1 ^ + H.c.) 

Jo 

+ - [/D (0 ^-^"^^^ 10 {e\ + T^.c] j doj (/ ^/y/2nbia + H.c.) , (27) 

where Z?^ is the annihilation operator for the mode of frequency oj in the channel continuum and a is the annihilation operator for 
the cavity mode. The energy of state \g) is set as zero. The \g) \t) transition is coupled to the cavity mode with strength gcav 
The \e} \t} transition is coupled by the external control laser of time-dependent Rabi frequency Q(0 and central frequency ojl. 
The coupling of the cavity mode to the channel continuum is assumed constant: ^y/2n. An ideal situation is assumed neglecting 
photon leakage into free space through intermediate state \t} or the cavity sidewall. 

We note that the system described by this Hamiltonian, under the laser excitation and the cavity-dot and cavity-channel 
interaction, has two invariant Hilbert subspaces, with the basis {\g, 0) |vac)} and {\e, 0) |vac) , \t, 0) |vac) , |g, 1) |vac) , \g, 0) |6i;)} 
(where in l^-,^), s = g,e or t denotes the state of the three-level system, n denotes the number of photons in the cavity mode, 
and \oj} denotes the one-photon Fock state of the channel mode of frequency oj). So the evolution of the system can be generally 
described by the state Cg \g, 0) |vac) -h Ce (0) in the interaction picture, where 

1^^ (0) = Pe (0 k, 0) Ivac) + A it) \U 0) Ivac) 

+ Pc (0 |g, 1) Ivac) + doj^a^ it) |g, 0) \oj) . (28) 
Jo 

The time evolution of the amplitudes in the interaction picture is described by the following Schrodinger equations, 

= g,,W^^'-^'^% + ^V., (29b) 

0c = -8ca.e-^^'"-'''^% - ^|^ain (0 - (29c) 

= -fev^"^"^"'""^^^A - ^|yc^out (0 + \Pc, (29d) 

where (0 = / dcoa^ito) e'''^'^''^'''^ I ^|2n with ^ and a^out (0 = / dcoa^iti) e~'^'^~'^''>^ I ^|27: with h -hoo can be 
regarded as the incoming and outgoing wavepacket of the photon in the quantum channel, respectively. From Eqs. \29c) and 
P9d| ), we note that evolution of Pc(J) is simply an instantaneous map of the difference between the input and output field in the 
photonic channel 

^|^Pc (t) = Qfout (0 - C^in (t) . (30) 

jStit) is also readily expressed in terms of ain(t) and aout(0 as, 

_^iKco,-co,)t^ (31) 



-(dput - ain)l ^|y - (Q^in(0 + Q^out(O) Vr/2 ,(^^_^^)^ 



8^ 



cav 



So as the amplitude of PeO)^ 



and the phase. 



|argOS,)=l^S,r2(M*-M*) 



Finally, from Eq. (29b ), the complex Rabi frequency of the laser pulse Q{t) can be expressed in terms of the amplitudes that have 
been solved above, 

n(t) = 2 ^'~^''''^ . (34) 
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Thus the desired operation, with ai^(t) and CKoutCO arbitrarily specified, can be generated on demand as long as the normalization 
condition of the wavefunction is not violated 

^ {\Pef + \Ptf + \Pcf) = |t^m(Ol' - koutCOP . (35) 

The functions of this quantum interface can be classified into three types: 

1. If there is no incoming photon, the quantum interface generates an outgoing photon wavepacket of a specified shape; 

2. If there is an incoming photon wavepacket of a specified shape, it is completely trapped by the quantum interface so that 
there is non outgoing field; 

3. If there is an incoming photon wavepacket of a specified shape, the quantum interface generates an outgoing photon 
wavepacket of another specified shape - a controlled scattering process. 

The first two types of control form the basis for the quantum network operation. With control of the type III, the quantum 
interface can act as a controllable scatter or pulse shaper for single photon wavepackets. This control can also be considered 
as the combination of consecutive controls of type II and I. In the following subsection, we will discuss in more details how to 
implement the first two types of controls for a quantum network. 



C. Inter-node operations in a quantum network 



The sending node of a quantum network is operated with control of type I. The initial conditions are: ain(t) = 0, jSdto) = 0, 
jSeito) = 1 and jStito) = 0. The integral form of Eq. ([35]) becomes 



= 1 - sin^ f l^out (r)f dr - - l^eavl"' \$c + 7fic/2f , (36) 

where o^out is the normalized wavepacket of the emitted photon, and sin^ 6 is the average photon number. For a photon number 
and a pulse shape arbitrarily specified, the amplitude of the cavity mode is determined by Eq. (30) as jSdt) = ^Jy&out(t) sin 6. If 
we pose the problem of finding the laser control pulse to produce a specified shape of the outgoing photon wavepacket, the fact 
that the right-hand side of Eq. ([36]) is positive requires the specified output pulse be sufficiently smooth, i.e., the pulse generation 
process be slower than the cavity-channel tunneling and the dot-cavity coupling rate (with time scales and g~^^, respectively). 
At the remote future time, ti -hoo, the photon emission process is completed, i.Q.,j3c(ti) = 0c(ti) = 0, so /3e(ti) = e^^ cos 6 with 
the controllable phase given by Eq. ([33]). The general form of the photon generation process can be expressed as 



a(0 



(C,|g) + C,k))0|vac) 

Cg\g) |vac) -h Ce [e'"^ cos e\e) |vac) -h sin e\g) l^out)] • (37) 
A full Raman process corresponds to ^ = 7r/2 and Pdh) = 0, where Eq. ( [37] ) is reduced to, 

{Cg\g} + Ce\e}) Ivac) ^ \g} (c,|vac) + C,|^out)) , (38) 

which results in the mapping of the stationary qubit onto the flying qubit. If initially the three level system is entirely in state 
\e), this mapping operation can function as a deterministic generation of a single-photon wavepacket with any desired shape 
Q'out(0- If the Raman cycle is controlled to be partially completed (0 < 7r/2), the state initially in \e) |vac) is transformed into 
an entangled state of the stationary spin and the flying photon 

\e} |vac) — > e'^P cos e\e} |vac) -h sin 0\g) l^out). (39) 

The entanglement entropy E = - cos^ log2 cos^ 6 - sin^ log2 sin^ 6 can be set at any value between and 1 depending on the 
rotating angle 0. 

The receiving node is operated with control of type II, typically as a full Raman cycle, in the quantum network scheme. With 
the three level system initially on state \g} and the incoming photon Qlvac) -h Celai^t)}, the mapping transformation is, 

\g} (C,|vac) + C,|a,n)) ^ (Cg\g} + Ce\e}) |vac). (40) 
As in the sending process, the incoming photon pulse ain(t) can be arbitrarily specified, provided that it is smooth enough. 
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By combining the sending and receiving processes, the transfer of a qubit from one node to another can be easily implemented, 
with the outgoing photon from the sending node directed as the incoming photon for the receiving node. As the photonic channel 
is linear, two state-transfer operations with opposite directions can be performed in parallel, and qubits at the two nodes will 
be swapped. Swap operations can only be performed between nodes separated with sufficiently large distance so that photon 
traveling time in the channel is longer than the interface operation time. 

If the operation at the sending node has been designed to produce an entangled state of the stationary and the flying qubit, the 
mapping process at the receiving node will just produce an entangled state of the two nodes by the transformation, 

\e}i\g}2 |vac) e''^cose\e)i\g}2 |vac) + sm0\g}i\g}2 l^out) 



[e'^cose\eh\g}2 + sin%)ik)2] |vac). (41) 



Non-local entanglement can thus be generated deterministically in the quantum network. 

With the exact solutions for the interface dynamics, the sending and receiving functions can be separately addressed. This 
enables the construction of a quantum network with heterogeneous quantum nodes, which is essential for solid state realization. 
To illustrate this, we give below an exemplary control strategy. The control laser pulse at a sending node can have a general 
shape Qi(0^~'^^^ with the pulse area satisfying the specified rotation angle 6. The outgoing single photon wavepacket has a 
definite shape c^l^^it) (in the rotating frame defined by frequency col), determined by the laser control Qi(0 and the parameters 
of the sending node gl^^, y\ and Raman detuning Ai = co] - - ojj - col. Then we pose the problem of finding the 

optical control ^2(0^"^^^^ to trap the single photon wavepacket £? (0 = ^^^^(O^^^^^"^^^^ at the receiving node which may have 



a diff'erent set of parameters g^av 72, ojf, co^ and co^. From Eq. ( |36[ ) and the discussion follows, it implies that the tolerance of 
the node inhomogeneity is determined by the node bandwidth, i.e. the dot-cavity coupling gcav and dot-waveguide tunneling y. 
Thus, a large dot-cavity coupling gcav is essential. This tolerance was discussed more explicitly by Fattal et al in Ref. 13491 . 

D. Operations with imperfections 

In this subsection, we discuss the eff'ects of various imperfections that may occur in a realistic quantum network and the 
corresponding mitigation. 

1. Intrinsic photon leakage into free space 

The desired quantum evolution is through the pathway of trion - cavity photon - waveguide photon. The main causes of photon 
leakage out of this quantum pathway is through the trion decay by spontaneous emission and the cavity mode leakage other than 
the tunneling into the waveguide. The waveguide and fiber loss is negligible on the distance- scale of relevance for intra-chip 
communication and for inter-chip communications if the chips are distributed in a spatial range of < 10 cm. As long as the 
photon leakage rate is much smaller than the bandwidth of the desired quantum pathway (determined by the dot-cavity coupling 
gcav and the cavity- waveguide tunneling y), high fidelity operations can be expected. Typical trion decay rates in self-assembled 
QDs are Y ~ yueV [|229[ 12301 , and the intrinsic loss rate of a high- 2 cavity (i.e. excluding coupling to the dot and the waveguide) 
can be potentially achieved at yo ^ 0-1 y^eV (corresponding to a Q-factor ~ 10^) |282, 283 1. The state-of-the-art dot-cavity 
coupling constant achieved is gcav = 0.1 meV 12791 1284 , 286], while the cavity-waveguide tunneling rate is controlled in design 
by the gap distance. Thus the bandwidth can be two orders larger than the leakage rate. 

The simulation result of mapping a spin state to a flying photon wavepacket with the shape targeted as an asymmetric su- 



perposition of two sech-functions as a^ouT^(^) ~ sech(7^/6 -h 5) -h 0.5sech(7^/6 - 5) is shown in Fig.l35| The trion decay rate, 
based on experiment |229 |, is set at F = 3yueV, and the intrinsic cavity loss rate is assumed to be yo = 0-1 y^eV |283 |. The 
cavity-fiber tunneling rate is chosen to be y = 0.2 meV and the dot-cavity coupling constant gcav = 0.1 meV | 284|. The fidelity 
of the single photon generation KcKou^Vout)! ~ 0.9912. Because of the non-adiabatic optical pumping and dot-cavity coupling, 
the whole mapping process can be completed within 300 ps. The simulation of the photon absorption process shows an over- 
all fidelity greater than 0.99 as well. With the above chosen parameters, the cavity mode broadening due to the coupling to 
waveguide is actually larger than the dot-cavity coupling. The system is therefore not in the strong coupling regime by the usual 
definition 12791 12841 1286J . High fidelity is nonetheless guaranteed by the large Purcell factor. 



Unwanted coupling to energy levels beyond the 3-level A system 



In the energy level structure of the dot-cavity coupled system, non-resonance coupling to other energy levels could lead to 
excitation out of the 3-level A subspace and AC Stark shift of the qubit states of interest (see Fig. 34 c)). These eff'ects have 
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been included in the numerical simulation given in Fig. 35 As shown in Fig. 35 c), AC stark shift induces a deterministic phase 
drift between \g) and \e) . As the two excitation pathways starting respectively from \g) and \e) are independent of each other 
(see Fig.l34tc)), this phase drift is independent of the coefficients Cg and Therefore, it can be compensated by a single-qubit 
operation irrespective of the quantum state being mapped. Leakage out of the qubit subspace by the non-resonance excitation to 
multi-photon states can be greatly suppressed if the Zeeman splitting is much larger than the Rabi frequency and the cavity-dot 
coupling. For InAs self-assembled QDs, Zeeman splitting ~ meV can be achieved in a moderately strong magnetic field (~ 10 T) 
due to the large g-f actor of InAs materials. 



3. Unknown parameter offsets 



Solutions to the laser pulses for desired controls of the quantum interface are based on the knowledge of the coupling strength 
^cav,7 and Q(0. But in practice, there could be various unknown errors on parameters due to imperfect characterization of the 
system. The robustness of the control schemes in presence of unknown system parameter errors is thus a critical feature. In 
Table. |lj we list the eff'ects of unknown off'sets from the assumed values of various parameters on the fidelity of two typical 
quantum network operations: (i) entangling two quantum nodes into state e^^\g)i\e)2 + k)il^)2; (ii) transfer of the state \g) -h \e) 
between two nodes. In both cases, the target shape of the involved single photon wavepacket is sech(^) and the design of 
the control laser pulses uses the assumed parameter values. The system shows a surprising robustness: 10% unknown errors 
on gcav,T or 1^(01 only reduce the fidelity by less than 1%. This intrinsic robustness against unknown parameter errors paves 
the way for learning studies of the system parameters by trial and error L320...322. 352.K and classical feedback controls in the 
quantum network 13531 13541 . 



4. Laser intensity fluctuations 



The |Q(OI error studied in Table. [l|is a global one on the amplitude, e.g. induced by the QD being slightly out of focus from 
the laser control field or an unknown off'set from the assumed value of the dot dipole moment. The actual control laser pulse 
may have shape deviations from the desired ones, e.g. due to temporal fluctuations in laser amplitude. The efl'ects of this error 



source have been studied in | 348 |, where the control scheme is found to be immune against fast fluctuations (see Fig. 36). This 
robustness is due to the finite bandwidth of the quantum interface determined by the cavity-dot coupling strength and the cavity- 
waveguide tunneling rate. Any temporal fluctuations in the control field with the frequency higher than the interface bandwidth 
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FIG. 35. (Extracted from Ref.| 119|) Generation of a single photon wavepacket with an asymmetric double sech shape (see text), (a) Real 
part of the dimensionless amplitude of the simulated single photon wavepacket (solid line) as a function of the dimensionless time ytl2. The 
deviation from the target shape (dashed line) is not visible, (b) Imaginary part of the simulated single photon wavepacket (solid line) and that 
of the target one (dashed line), (c) Phase drift of the component \g, 0>. (d) Rabi frequency of the control laser pulse. 
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are effectively averaged out. The time independent amplitude error listed in Table. |T| may also be considered as a special shape 
error which is actually the worst scenario for the fidelity. 

5. Laser phase fluctuations 

The complex Rabi frequency Q(0 can also have phase uncertainties due to laser phase drift which is typically slower than 
the time scale of interface operation. Assume that the laser control field at the moment of sending operation has an unknown 
phase of (fi and, hence, the Rabi frequency is Di (0 e^'^^ . From the form of the coupling term that involves the laser field in the 



Hamiltonian Eq. (27 ), the unknown phase factor can be absorbed by redefining the state \e)i = e ^'^^ \e)i, so that, 

i [i^i(t)e'^^e-'^''\t)n{e\ + H.c] ^ ^ [/Qi(0^-^""^^|Oii<^l + H.c] (42) 

and we can make the same transform at the receiving node where the control laser has the unknown phase (p2 at the moment the 
single photon arrives. Starting from a general state {Cg\g)i + Ce\e)i^ \g)2 |vac) = Cg\g}i\g}2 |vac) + Cee^'^^\e)i\g)2 |vac), 
two-node operations in the presence of laser phase uncertainties can be generally expressed as, 

Cg\g}i\g}2 Ivac) + Cee^^^\eh\g}2 |vac) 

C^|g)i|g)2 0|vac) 

+ Cec'^^ p"^ cos^|^)i|g)2 Ivac) + sme\g)i\g)2 l^out)] 

""^'^ [C,\g}i\g}2 + Cee^^' {e^^cose\eh\g}2 + sin %)i 1^)2)] |vac). (43) 

The final state is equivalent to [C^|g)i|g)2 + Ceie^"^ cos^|^)i|g)2 + e^'^^e'^'^^ sin ^|g)i 1^)2)] |vac). If the control laser fields at the 
two nodes is phase locked so that there is a certain relative phase between Qi(t) and fl2(t - t), the two-node operation is well 
protected from laser phase fluctuations. 

6. Deterministic phase and shape variations in photon propagation 

Unknown off'sets from the assumed waveguide/fiber dispersion relation or nonlinearity of the dispersion relation can cause 
phase and shape variations of the single photon wavepacket during the propagation. At low temperature where thermal fluc- 
tuations are suppressed, such variations are deterministic and can thus be incorporated in the design of the receiving node 
control. Close-loop adaptive feedback control 13531 13541 or quantum learning algorithms 13201 13221 13521 can be implemented 
to pre-characterize such variations. 

7. Loss and indeterministic fluctuation in photon propagation 

For inter-node operations between well separated nodes (distance > m), photon losses and indeterministic fluctuations during 
the propagation in the fiber could be non-negligible. Error correction schemes dealing with such propagation loss has been 
proposed using auxiliary stationary qubits in the quantum node 13551 . The idea of quantum repeaters with nested purification 
schemes | 356| might also be incorporated into the quantum network design for protection against the photon propagation loss 
and indeterministic fluctuations. 

E. Summary: Coherent quantum manipulation by remote control 

In this section, we have discussed how to unite clusters of QD spins into a network for distributed quantum information 
processing. The advances in the fabrication of coupled structures of QDs, semiconductor micro-cavities, and optical waveguides 



TABLE L Effect of errors in coupling parameters on the fidelity of entanglement and of state transfer. 





no error 


10% g error 


10% 7 error 


10% Q(0 error 


Entanglement 


0.9912 


0.9872 


0.9894 


0.9862 


Transfer 


0.9901 


0.9870 


0.9891 


0.9879 
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make possible a high efficiency quantum interface between stationary spin qubits and flying photon qubits, where the latter 
can be used for communication between seperated clusters. As indicated in the DiVincenzo criteria |2|, such capability greatly 
enhances the chance towards the construction of fault- tolerant scalable quantum computers in these systems. 

Counterintuitively, real excitations of the continuum modes in the quantum channel act as a conduit for mediating coherent 
operations rather than the cause of dissipation. A more abstract picture for the network structure here is the 3 -level A system 
coupled to a one-dimensional continuum modes, formed by the cavity-waveguide coupled structure, considered as a whole 
electromagnetic continuum. The important function of the cavity is in providing a local resonance in spectrum with large 
spectral weight near the frequency of the QD optical transition, which results in a large Purcell factor (i.e. the ratio between the 
spectral weight of the cavity-waveguide modes and the free-space modes). The unitary evolution dominates over the irreversible 
processes within the coupled system of the QD and the cavity-waveguide continuum. Temporal shaping of the laser control 
pulse provides sufficient freedom to control such evolution. While the relevant part of the spectrum of the cavity-waveguide 
continuum is of a simple Lorentz shape with width y <^ ojc, this Markovian condition is not a mandatory requirement. Provided 
that the 3 -level A system is coupled to a continuum where a spectral weight peak (irrespective of shape) results in a high Purcell 
factor at the relevant optical frequency, control schemes for high efficiency network operation is possible 13571 . This opens 
up possibilities for networking localized stationary qubit by continuum modes without cavity QED and in the non-Markovian 
regime. 



VIII. CHALLENGES 

We have discussed all the necessary elements for implementing scalable quantum computation with electron spins in QDs 
under optical control. To reach the larger goals, there remain many obstacles. Here we present an overview of what technologies 
may be at the top of the required list to accomplish the goals. To have a concrete idea of the challenges for short- to mid-term 
pursuit, we will also give an estimate of the resources for a benchmark task: factorization of 15 with Shor's algorithm. 



A. Resource estimate for Shor's algorithm for factorizing 15 

Undoubtedly, Shor's algorithm for factorizing integers is the most important example demonstrating the superpower of quan- 
tum computation [Ti l 17211 . It gives a solution with time consumption only polynomially increasing with the problem size (the 
bit length of the number to be factorized) and thus ofl'ers an exponential speedup over all known classical counterparts. Histori- 
cally, Shor's algorithm has stimulated the exploding enthusiasm in quantum computation by showing its computation power. An 
efficient factorization scheme can be used to break the public-key encryptions such as the RSA protocol which is widely used 
in internet communication. The factorization of the first "non-trivial" number - 15 by Shor's algorithm has been realized with 
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FIG. 36. (Extracted from Ref. |348 |) Transfer of state |g> + \e) between two nodes in presence of amplitude fluctuations of the control lasers, 
(a) The desired control laser pulse by design at the sending node (solid curve) and the actually applied control pulse with slow fluctuations 
(dashed-dotted curve) and fast fluctuations (light curve). The control laser pulse at the receiving node has a similar error applied, (b) The 
generated intermediate single photon wavepacket. The solid curve shows the target shape of the single photon if the ideal control laser pulse is 
used (solid curve in (a)). Deviation of the generated single photon wavepacket is invisible in the case of fast laser fluctuations. The fidelity of 
the transfer is 0.9912 with slow fluctuations and 0.9922 with fast fluctuations. The parameters used are: y = 0.2meV, gcav = 7/2, = 3//eV, 
and intrinsic cavity leakage 0.05 yt/eV. 
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the liquid NMR spectroscopy fT|, which serves as a benchmark for quantum computation in other systems. Here we give an 
estimate of the resources required to accompHsh such a milestone in the optically controlled spin-based QD system. 

The key step of Shor's algorithm for factorizing a number N is to find the smallest nonzero r satisfying = I mod(N) 
(i.e., the order of a with respect to N), where the seed a is co-prime to N [i.e. the maximum common factor of a and N, 
denoted as (a,N) is 1] and can be randomly selected. By the fact that the order r is the period of the modular exponentiation 
function f(x) = mod(A^), it suffices to search for the period or to measure the frequency of the function f(x), which can be 
accomplished by the quantum Fourier transformation (QFT). The process of finding the period can be expressed as dEHl 

_^ 2"-l 2«-l 

|0") lO'^) Yj I-^> ^ l^"") ^ X '-^^ ^ '^"^ mod(N)} 

2«-l 2«-l 

^ YjYj ^^''''^'^^ ® l^"" mod{N)), (44) 

where |0"^^^) denotes a n(m)-qubit register set at the zero state, m = [log2 A'^] + 1 is the bit length of A^, x is a binary number, and 
the normalization constants for the states have been omitted. Both the QFT and the modular exponentiation can be carried out 
with the number of elementary quantum gates polynomially increasing with the problem size (m). A more detailed review of the 



order-finding algorithm is given in Appendix C 1 



As f(x) is a periodic function with period r, the population of the basis state of the first register at the end of the program 
will be non-zero only for y ^ 2"c/r (c = 0, 1, . . . , r - 1), for the QFT will transform a period function into sharp peaks at the 
multiples of the frequency 2"/r. Thus, a quantum measurement of the first register in the computational basis will result one of 
states in the superposition, thus yielding the number 2^c/r from which the order r can be derived and the number N factorized. 
The order r generally is of the order of magnitude of A/^ ~ 2'^, so the number of basis states in the superposition is ~ 2^, an 
exponential function of the problem size, which explains why a QND measurement is required to make Shor's algorithm scalable 



for factorizing large numbers, as discussed in Sec. VI 



The flow-chart of Shor's algorithm for order finding is shown in Fig. [37] The first QFT subroutine can be simplified to a series 
of single-qubit Hadamard gates (see Appendix |C 1| for the definition of the elementary quantum gates) since the initial state is 
set to be |0"). The second register is initialized to be |1) to facilitate the modular exponentiation. In principle, the length of the 
first register n should be large enough to reproduce the real number 1/r with sufficient eff'ective bits. To factorize 15, it turns out 
that the orders of all co-primes of 15 are either 2 (for a=4, 11, 14) or 4 (for a=2, 7, 8, 13), both of which are a factor of 2^, so 
three-qubit register should be enough to resolve 2"/r. Including the second register, 7 qubits are sufficient to demonstrate the 
algorithm for the first "non-trivial" target: 15. 



Figure[38]and|39]give the specific quantum circuits constructed with the elementary quantum gates (as defined in Appendix[D|) 
to implement Shor's algorithm for factorizing 15 with the seed number selected as a = 4 (an easy example) and a = 13 (a difficult 
example). Compilation and optimization have been performed to reduce the number of physical operations: 

1 . To minimize the number of SWAP gates for coupling remote qubits in the linear configuration, the labelling of the qubits 
in the register has been optimized, and at the end of the program the qubits are not swapped back to their original positions 



but re-labelled instead, as indicated in Fig. 16 (a). 



2. The quantum gates are omitted if they act on the second register after the last gate controlled by the first register, since 
such gates have no effect on the QFT of the first register. 

3. Whenever possible, the sequential single-bit operations on the same qubit are combined into one single-qubit gate (in 
experiment, all one-bit gates between two controlled gates can be realized by a single-spin rotation). 

The numbers of gates and pulses and the running time needed in different cases are estimated with the following considerations: 




FIG. 37. Logic flow-chart of Shor's algorithm for order searching. 
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TABLE IL The numbers of gates and pulses and the time required to factorize 15 by Shor's quantum algorithm with two typical examples of 
seed numbers a for - \ mod(15). 



a 


# 1-bit gates 


#SWAP 


# phase gate 


# pulses 


time-cost 


a=4 


4 


1 


3 


44 


0.9 ns 


a=13 


20 


8 


15 


142 


1.3 ns 



1. All sequential single-qubit gates on the same qubit are combined as one single-spin rotation which can be completed with 
three pulses within 10 ps. 



2. Each SWAP gate can be completed by two pulses within 10 ps. 



3. Each controlled phase gate can be completed by three pulses in 10 ps. 



4. Each qubit can be initialized by three pulses within 120 ps. 



5. To reduce the complexity of pulse designing, all operations within the local module should be performed serially (to 
factorize larger numbers, a large-scale quantum network of distributed modules will be needed and in that case operations 
could be performed in parallel in separated modules). 



6. The time consumption for quantum measurement is not counted. 



The timescales of the single- and two-qubit gates are chosen in accordance with the discussions in Sec. IV While there are so 
far no experimental demonstrations of optically controlled two-qubit gates, the estimated timescales of the single-qubit gates are 
similar to the experimentally realized gates l78| [79lfT8nil82 1. The estimated resource requirements are summarized in Table [11] 
Alternatively, the controlled phase gates can be replaced by the VSWAP gates, and the resource requirement is expected to be 
similar. 
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FIG. 38. (a) Quantum circuit for searching r for - \ mod(15), where the subprograms for the modular exponentiation and the QFT are 
blocked together as indicated, (b) The simplified circuit with only local two-qubit gates. 
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B. Technologies most needed 



1. Complex multi-pulse optics 



As listed in Table |ll| to factorize the small number 15, the number of optical pulses needed is in the order of 100. It is would 
be extremely difficult to synchronize so many laser pulses and to stabilize the relative phases, but all the pulses may be viewed 
as a single one with complex frequency and phase design. The complex pulse may be formed from a single laser source by pulse 
shaping techniques, such as acousto-optical modulation 13581 13591 . Design of such a complex pulse is also very challenging. 
Learning or genetic algorithms may be developed to deal with the problem, with lessons learned from controlling complex 
chemical processes by ultrafast optical pulses l320H322 i. The measurement involves many probe pulses interacting with many 
QDs. The signals are to be analyzed through frequency multiplexer as well as phase modulation with heterodyne detection | 331J. 
Thus multi-dimensional spectroscopies I360H3631 are desirable for characterizing the systems and for implementing a small- 
scale benchmark demonstration. Eventually, the measurement would have to be done with an efficient quantum channel like a 
cavity- waveguide structure. 



2. System fabrication and characterization 



The design of the quantum computer and its operation follow the sequence: The physical system is constructed first, with 
certain uncontrollability. Then it is characterized. The optical control will be designed according to the system parameters. In 
such a procedure, we do not require an ultimate control of the system fabrication (hardware) but defer the difficult work to the 
control design (software) stage. For example, we do not require, and actually do not desire that all the QDs are almost identical. 
But of course, the system should be fabricated fulfilling certain conditions. Basically, we need a system made of QD clusters, 
which should meet the following requirements: 

1. QDs are only locally coupled. For a system with coupling between remotely separated dots, the design of the controlling 
pulses would require overhead increasing exponentially with the number of qubits, since essentially the control design 
amounts to solving the Schodinger equation of the whole system. 

2. The coupling between diff'erent QDs in a local node should be weak. Otherwise, tunneling between diff'erent dots would 
make the local node rather a large QD molecule instead of a cluster of individual dots, making even a single-qubit gate as 
complex as a control of all the qubits in the local node. 

3. Even though we do not require all the QDs to be nearly identical, the size of diff'erent QDs in a cluster should be in 



|o)m]_ 

|0)^hH^ 
|0)HhH- 



145 



-e- 



90j-|Hfr)" 
QFT I 



13" mod 15 



(a) 




FIG. 39. (a) Quantum circuit for searching r for 13^ = 1 mod(15). The faded components can be removed without changing the computing 
result, (b) The simplified circuit with all gates experimentally realizable. 
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a relatively small range. Otherwise, the optical control, especially the energy shift by the AC Stark effect, would be 
extremely difficult. 

4. The number of active electrons in each QD should be controlled to be one, either by doping or by gate voltage control. 

Once a QD system has been constructed, the characterization is not any less demanding. The parameters to be determined 
include electron doping level, electron spin g-factors, frequencies of the ground and excited state transitions, selection rules 
or relative dipole matrix elements for different transitions, tunneling rates and exchange interaction of electrons in different 
levels and different dots, etc. Identification of the exciton, bi-exciton, or trion transitions in two-dot systems have already been 
demonstrated in two vertically coupled dots with the help of varying the gate voltage 11951 . A pair of laterally coupled dots 
remains a challenge. The characterization of a cluster of, for example, seven QDs, is a tremendously difficult target for the 
current experiment capability. In the long run, we would expect no full characterization of the system be required. The solution 
may lie in, again, the learning algorithm, by which, the controlling optical pulses are to be designed, self-adaptive to the fidelity 
of a certain set of quantum gates. 

3. Nano-photonics 

The requirement of nano-optics is two-folded: First, the control of a local node requires near-field optical addressing. Sec- 
ond, the ultrafast fast initialization and QND measurement of a qubit, and coupling between different nodes need cavities and 
waveguides fabricated in-situ with the QDs. 

To address individually each cluster of QDs, a micro-lens and micro-fiber may be used. Each QD within a cluster is dis- 
tinguished by its fingerprint transition frequencies. The spatial resolution required is given by the distances between clusters 
instead of the size of the clusters. Using high-index material for the micro-lens, spatial resolution -0.1 yum may be achievable. In 
cases where the spatial resolution is not enough to single out dot clusters, an alternative solution would be further pulse shaping 
(probably with a learning algorithm as well) to eliminate the coupling between different clusters covered by one micro-lens. 

To connect QDs with cavities and waveguides, two structures look promising. One is obtained by etching the surface where 
the QDs are grown |278|. Electron-beam lithography and chemical etching (sometimes plus some annealing) have already 
produced high-quality microcavities and waveguides on semiconductor surfaces. Even strong coupling between a QD and 
microcavity in such a system has been demonstrated 12791 . Photonic crystals are another promising possibility. Point defects 
in photonic crystals can be made into nano-cavities with Q-factor > 10^ and effective volume less than a half- wavelength 
cubed 12821 1283li . Strong coupling between a QD and a photonic crystal nano-cavity has also been demonstrated 12841 12851 . 
Quantum electrodynamics of single QDs in nano-cavity in photonic crystals can be engineered |281 , 304, 305 1. Waveguides in 
photonic crystals can be made by line defects which may be coupled to remote nodes by optical fibers | 301 1. How to combine 
these photonic structures with the QD clusters, especially, how to assemble them in proper layouts and positions, would demand 
a great deal of progress in sample processing technologies. 

We would like to point out that the photonic structure fabricated on the QD system may also be used to individually address 
each local node. Thus, we eventually may need no micro-lens or micro-fibers to attach each dot cluster. 

A wilder conjecture is that the lasers be integrated in the photonic structure. Micro- or nano-lasers made of photonic crystal 
cavities have actually come into being. The remaining problem would be to make the laser emit into a waveguide, and to tailor 
it into wanted shapes and sequences. How is one to control the laser on-chip then? Electrical gates may be used. So we come 
around a full circle and find a point where different quantum computation strategies may be synthesized to achieve a common 
goal. 

IX. CONCLUSION 

We have discussed various aspects of a scalable scheme of quantum computation based on optical control of electron spins in 
semiconductor QDs. To implement such a scheme, a number of outstanding challenges remain to be overcome. 

We would like to add some remarks on two different philosophies in implementing digital computing or more generally 
automatic reasoning, which may give us some inspiration in the journey to realize quantum computation. One is related to the 
von Neumann structure of computers and the other is related to the Turing ACE structure, both of which are based on Turing's 
insightful view of programs and data as essentially the same for a universal computer. Turing's design is deeply rooted on his 
finding of universal computers and thus has a hardware of minimal instructions with complex functions to be implemented by 
software programming. The von Neumann structure, while keeps as well the important role of programming, tries to maximize 
the usage of hardware design to implement a large number of mathematical functions which would otherwise be solved simply 
by programming. We follow Turing's perception in describing the blueprint of a quantum computer based on semiconductor 
QDs and optical control of spins in them, simply for one reasons: In implementing quantum computation, the hardware part is far 
more formidable than the design of optimal control. Thus we propose no need of perfectly controlled arrays of almost identical 
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QDs but a control scheme programmed after a physical structure has been constructed and characterized. The randomness in the 
system synthesis is not to be eliminated but to be utilized. The requirement for a functional physical block is relatively simple: 
sufficient coupling for a universal set of quantum gates to be programmed. Finally, we note that Turing's philosophy is becoming 
more and more important nowdays in large-scale classical computers which tend towards the RISC (reduced instruction set 
computer) architecture. 

When will a quantum computer come into being of practical usage? We do not know. But a hint may lies in the comparison 
between our present situation and the situation something 60 years ago when engineers were working hard to maintain thousands 
of vacuum tubes functioning together for a while before one or another went wrong. Scaling up of a quantum computer may be 
not as rapid as classical computers have done, but just be aware that adding one functioning qubit supposedly doubles the power 
of a quantum computer, which is worth 18 hardworking months in the sense defined by Moore's law. 
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Appendix A: Hole-mixing and selection rules in a quantum dot 



In the III-V bulk material, the top of valance band occurs at the T point of the Brillouin zone (see Fig. 40). The spin-orbit 
coupling of the /7-orbit splits the p-type valance band to a quadruplet with symmetry (or a total angular momentum J = 3/2) 
and a doublet with symmetry (or a total angular momentum / = 1/2). The T-l band lies much lower in energy than the 
r-8 bands. From now on, we focus on the bands derived from the J = 3/2 quadruplet, which are described by the Luttinger 
Hamiltonian in the framework of eff'ective-mass theory 13641 



Hr = 



2mo 



( Pi Q R 

Pi R 

7?* P2 -Q 

7?* -2* Pi J 



expressed in a matrix form in the basis {|/^ = 3/2), = 1/2), = -1/2), = -3/2)}, with 

Pi = (n - 272) kl + (n + 72) {kl + k^y) , 
Pi = (n + 272) kl + (71 - 72) {kl + k^y) , 
Q = -2^/3y3k,{k,-iky), 
R = -^[y2{kl-k^y)-2iy3k,ky], 



(Al) 



(A2a) 
(A2b) 
(A2c) 
(A2d) 




conduction band 



heavy hole band 
light hole band 
split-off band 



FIG. 40. Bulk band structures of direct bandgap III-V semiconductors. 
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where 71,2,3 denote the Luttinger coefficients. The = ±3/2) and = ±1/2) bands are usually referred as the heavy hole and 
light hole bands respectively by their different effective mass in the z direction. 

Hole mixing effects are caused by the QD confinement potential Vh{x,y,z). Here for simplicity, the growth direction z is 
assumed to be along the [001] direction. For other growth directions, the discussions below can be generalized straightforwardly, 
and the essential conclusions would not be changed. Strong quantum confinement in the growth (z) direction lifts the four-fold 



degeneracy at Y point (see Fig. 41 a)). The light hole states actually have a larger in-plane effective mass and therefore their 



dispersion tend to cross into the heavy hole bands at finite in-plane wavevector (see Fig. |4TJb)). Spin-orbit coupling thus 



results in anti-crossing and the mixture of the heavy-hole and light-hole bands (see Fig. 41 c)) | 364 |. 



Let us consider first the zeroth-order approximation for the hole-mixing. When the confinement size along the z-direction 
is much smaller than the lateral size, we have {k^) » (^^), {ky). If the off-diagonal couplings Q and R are neglected, the 
Hamiltonian is diagonal in the basis of the angular momentum quantized along the growth direction. The heavy and light holes, 
with kinetic energy Pi and P2, are characterized by their angular momentum states as | ± 3/2) and | ± 1/2), respectively. The 
heavy hole (HH) and the light hole (LH) are separated by an energy IS^hl = '^72{k^}/^o |364|. The hole-mixing are induced 
by Q and R terms. The Q terms couple | ± 3/2) to | ± 1/2), and the R terms couple | ± 3/2) to | =F 1/2). As we usually have 
» (2), (P), the zeroth-order approximation and the corresponding selection rules determined by the angular momentum 
conservation are often adequate to understand the optical transitions. 

For small QDs, the mixing may be important. Now let us consider the HH-LH mixing in different situations. 

First we consider a confinement potential with rotational symmetry about the growth direction. If we assume 72 =73, the 
Hamiltonian has the rotational symmetry. The HH and LH states coupled by the Q and R terms must have different orbital 
angular momentum. So the light-hole components mixed into, e.g., the heavy hole ground state are not optical active. In this 
case, the mixing has no effect on the optical transitions but a reduction of the dipole matrix element. The problem comes from 
the fact that in reality we usually do not have a cylindrical V(x, y, z) or 72 ^ 73 (since the lattice has no spherical symmetry). 

Let us still make the assumption that the confinement potential V has inversion symmetry in the z-direction, which is reason- 
ably fulfilled in fluctuation QDs. Since the Q terms are linear in the k^, they couple HH and LH states with different parity which 
are separated by large energies due to the strong confinement in the z-direction. In the even parity ground state, the hole mixing 
effect caused by the Q terms is negligible, especially in optical transitions. Thus, we need only to consider the coupling induced 
by the R terms. Thus the optically active components of the HH ground states are 

|/f±) = |±3/2) + 7/| + l/2), (A3) 

where rj ~ L^/L^^y (with and L^^y being the vertical and lateral confinement sizes, respectively) and normalization is under- 
stood. The dipole matrix elements between the electron spin states | ± 1 /2) and the trion states |^±) which contains the HH in 
|^±) are (in arbitrary units) 

<r±|d-o-±|± 1/2) = 1, (A4a) 
{t±\d'0-^\±l/2} = T]\ (A4b) 
(r=F|d-o-±|±l/2) = 0, (A4c) 
<r=F|d-o-^|±l/2) = 0, (A4d) 



where cr+ is the circular polarization of a light normally propagating and 7/' ~ rj. The selection rules are shown in Fig. 42 The 
two subspaces {| + 1/2), \t-\-)} and {| - 1/2), |^-)} are disconnected under optical coupling, which means the electron spin along 
the growth direction is still conserved unless an external magnetic field is applied along an in-plane direction. With a magnetic 
field along the x-direction, we can choose a light with the polarization, e.g., 

&+ = 0-+ - ri'cT-, (A5) 




FIG. 41. Dispersion relation of the = 3/2) and 1/^ = 1/2) valence bands in a quantum well }36A\. (a) Energy shift due to the confinement 
in the growth direction, (b) In-plane dispersion of the two valance bands when the "off-diagonal" terms are neglected (without mixing), (c) 
Level anti-crossing when "off-diagonal" terms are included. 
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which couples the electron spin eigen states | ± x) to the common trion state |^+) with renormalized dipole matrix elements 
(l ±T]'^y The theories for the case without hole mixing can be generalized simply by replacing the dipole matrix elements there 
with the renormalized. 

In systems without the inversion symmetry, which is often the case for InAs self-assembled QDs, the mixing would make the 
two ground HH states to be 

m) = I + 3/2) + f I + 1/2) + 7/1 - 1/2) + ^1 - 3/2), (A6) 
\H-} = I - 3/2) + f I - 1/2) + 7/1 + 1/2) + ^1 + 3/2), (A7) 

where t] and f are small numbers (< 10% for typical self-assembled QDs), and ^ ~ t]^. Now the dipole matrix elements are 

<r±|d-o-±|± 1/2) = 1, (A8a) 

{t±\d-cr^\±l/2) = T]\ (A8b) 

{t + \d-cr^\±l/2) = C, (A8c) 

{t + \d'0-J±l/2} = ^\ (A8d) 

with rj' ~ 7/, f ' ~ f , and ~ ^, all being small numbers. Interestingly, all the states now are optically connected. Thus it is 
possible to perform an arbitrary spin rotation by optical control without applying an external magnetic field or tilting the light 
beam from the normal direction. (For a QD without a symmetry axis, the normal direction is not special). For example, if the 
light polarization is chosen to be X = o-+ the dipole matrix elements between the spin states | ± x) and the trion states 

\t±x} = \t-\-} ± \t-} are 

jc|d • XK jc) = 1 -h 7/' -h -h ^\ (A9a) 
{t -h jc|d • X| - jc) = 0, (A9b) 
{t - x\d • X| -h jc) = 0, (A9c) 
{t-x\d'X\-x} = l-\-T]' - (A9d) 

Thus the AC Stark eff'ect induces an eff'ective magnetic field along the x-direction, which, as compared the eff'ective magnetic 
field along the growth direction induced by a circularly polarized light, is reduced by a factor ~ ' -h ^0- To realize an arbitrary 
rotation without a static magnetic field, the optical field would need to be much stronger than with a static magnetic field, unless 
the QD lateral size is comparable to the vertical size (such as in the case nanocrystals or spherical QDs formed by chemical 
deposition). 

The heavy-light hole mixing coefl&cients f and r] can be extracted using the polarization dependent absorption spectroscopy as 



described in Ref. 13651 . For typical InAs self-assembled dots studied by Steel's group, f is negligible and t] varies in the range 
~ 0.1 - 0.2 for diff'erent dots l[T?7ll257l . 



Appendix B: Theory of electron spin decoherence by interacting nuclear spins 

In this appendix, we give the details on the theory of electron spin decoherence by interacting nuclear spins in a strong 
magnetic field (> 1 Tesla), which was first formulated in Ref 1 127 |. 

The mesoscopic system consists of an electron with a spin vector Se and N nuclear spins, J„, with Zeeman energies and 
(jjn under a magnetic field ^ext. respectively, where n denotes both positions and isotope types (e.g. ^^As, ^^Ga and ^^Ga in 
GaAs). The interaction can be separated as "diagonal" terms which involve only the spin vector components along the field (z) 
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FIG. 42. Selection rules of optical transitions in a doped QD with inversion symmetry along the growth direction. The solid lines denote the 
allowed transitions, and the dashed lines denotes the "forbidden" transitions which have a relative dipole matrix element t]' . 
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direction and "off-diagonal" terms which involve spin flips. Because the electron Zeeman energy is much larger than the strength 
of the hyperfine interaction, the ofl'-diagonal term is eliminated by a standard canonical transformation, with the second-order 



correction left as the hyperfine-mediated nuclear interaction I127L called the extrinsic interaction in Section II D 1 For the same 
reason, the ofl'-diagonal part of the nuclear interaction contributes only when the terms conserve the Zeeman energies (so-called 
secular terms in the NMR terminology). Hence, the non-secular terms are negligible. The total reduced Hamiltonian is obtained 
for the limit of long longitudinal electron spin relaxation time (Ti oo), 

^red = 4 + + X l±)^^<±l, (Bl) 

+ 

with = ^e^e' = oJnJn, and the interaction terms, 

= ±Ha +Hb + Hd± He, (B2) 



given by, 



= Y;^J:J- = J]'a„,./:4, (B3a) 



4Qe 



Hb = Y,^nJlr,n^ (B3b) 

nt-m 

Hd = YPn,mJlP,n, (B3c) 

n<m 

HE = Y,^J'n^Yu^nJl, (B3d) 



where |±) are the eigenstates of ^g, the summation with a prime runs over only the homo-nuclear pairs, the subscript A denotes 
the extrinsic hyperfine mediated interaction, B the off'-diagonal part of the intrinsic nuclear interaction, D the diagonal part of 
the intrinsic interaction, and E the diagonal part of the contact electron-nuclear hyperfine interaction. The hyperfine energy, 
determined by the electron wavefunction, has a typical energy scale En ~ a„ ~ 10^ s"^ for a dot with about 10^ nuclei 13661 . 
The sum, J{ = Yjn^nA^ the hyperfine constant depending only on the material. The intrinsic nuclear spin-spin interaction has 
the near-neighbor coupling strength Bn,m ~ ^«,m ~ 10^ s"^ The extrinsic hyperfine-mediated interaction, which is unrestricted 
in range within the QD and associated with opposite signs for opposite electron spin states, has an energy scale dependent on 
the field strength, An,m ~ 1-10 s"^ for field ~ 40-1 T. Thus, the intrinsic interaction dominates local pair-interactions, while 
non-local pairs are driven by the extrinsic mechanism. 



a. Formal theory of de coherence in a nuclear spin bath 



After the initialization step of the electron spin qubit, the electron-nuclear spin system is prepared in a product state with the 
nuclear spins in a thermal state with temperature T, described by the density matrix 

p(O)=p^(O)0p^. (B4) 

The time evolution of the reduced density matrix of the electron spin, 

p\t) = TrNP(0, (B5) 



A'' 



+ 



1/2) 



— 1-1/2) 



FIG. 43. The relative dipole matrix elements in a doped QD with irregular shape but a large lateral to vertical size ratio. The light polarization 
is assumed to cr+. The solid lines denote the allowed transitions, and the dashed lines denotes the "forbidden" transitions which have relative 
dipole matrix elements ri\ ^' , and ^' . 
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is obtained by tracing over the nuclear spins, may be expressed in the form, 

pIM = ^,,v-,;v' {t)pl,y (0) (B6) 

where p^^ = (yu|p^|y), and |yu), |y) g {|+), |-)}. The superoperator or correlation function JL,_i^y-,_i' y can be expressed in terms of 
the evolution operator and contains the information on the electron spin relaxation and decoherence. 

The Hamiltonians of Eq. (Bl ) for the Ti ^ oo limit conserves the electron quantum number: [H.Sl^ = 0. Hence, the 
correlation function has following properties, 

■Cfi,v;fi\y(t) = £fi,vO)^fi,fi'^vy, (B7a) 
J:^At) = 1' (B7b) 

£^,-(t) = i;;+(o, (B7c) 

and the specific expression for the free-induction decay (FID), 

£^,.(t) = e-'^'^'Tv^ [p'^e^'^'^e-'^^'] , (B8) 

which can be straightforwardly extended to dynamics under pulse control. 

The ensemble of nuclear spins, at temperature T > ojn ^ ^n,m^ ^n,m, J^n,m, En, may be approximated by the density matrix, 

J 

where \J'} = \jn}, jn being the quantum number of nuclear spin n in the magnetic field direction. Pj is the thermal 
distribution factor. While single-system dynamics (i.e., with the nuclear bath initially in a pure state \J'}) could be the ultimate 
aim for quantum applications, we note that all experiments to date are performed under the ensemble scenario, either in a spatial 
ensemble of many dots 11091 11771 12Q5I or in a single dot with repeated measurements II1121 11611 12Q2I , where the statistical 
average over the initial configurations is needed. 

The correlation function X+ -(0 can then be generally expressed as, 

i;,,_(0 = Pje-'^^^'^ \{J-(t)\J^(t)}\ (BIO) 

J 

In FID, \J'-(t)} = e'^^'^lJ'} and (pj(t) = (Q^ + 8j)t where Sj = 2« jn^n is the contribution to the electron zeeman splitting 
from the Overhauser field in the nuclear configuration 



b. Pair correlation approximation and pseudo-spin picture 

The solution to the single-system evolution \J'-(t)} is key to both single spin decoherence and ensemble decoherence behaviors 
under free induction decay and pulse controls. Due to the slowness of the nuclear spin interacting dynamics, this evolution is 
well described by the pair-correlation approximation for the nuclear spin bath I127II129J . Within a time t much smaller than the 
inverse nuclear interaction strength, the total number of pair-flip excitations A^flip is much smaller than the number of nuclei N. 

The probability of having pair-flips correlated can be estimated to be Pcorr ~ 1 - e'^^^'p^^ (q being the number of homo-nuclear 
nearest neighbors), which is negligible in the relevant timescale of electron spin decoherence 11271 11291 113011134112201 . Thus, 
the pair-flips as elementary excitations from the initial state can be treated as independent of each other, with a relative error 
^ ^ ^corr- Then the single-system dynamics \J'-(t)} can be described by the excitation of pair-correlations as non-interacting 
quasi-particles from the "vacuum" state \ J'), driven by the "low-energy" efl'ective Hamiltonian, 

H'j = Z'^^^Z^t-^k/2, (BID 

k k 

which has been written in such a way that the pair-correlations are interpreted as 1/2-pseudo-spins, represented by the Pauli 
matrix &k, with k labeling all possible flip-pairs 11341 . The time evolution from the initial state \ J') can be viewed as the rotation 
of the pseudo-spins, initially all polarized along the -\-z pseudo-axis: | T^t), under the efl'ective pseudo-magnetic field, 

= (±2Ak + 2Bk, 0, Dk ± Ek\ (B12) 

where, for the electron spin state |±), ±Ak and are the pair-flip transition amplitudes contributed by the extrinsic nuclear inter- 
action Ha and the intrinsic nuclear interaction Hb, respectively, and Dk and ±Ek are the energy cost of the pair-flip contributed 
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by the diagonal nuclear coupling Hd and the hyperfine interaction He, respectively. The decoherence then can be analytically 
derived as 

X;_(0 = n |(^,-(0l^,"(0)| ^ n e<l\ (B13) 

k k 

where 10^^(0) are the two conjugated states of pseudo-spin k at time t conditioned on the electron spin state |±). In FID, 

= e-'^k^\ Ty^); while with a n pulse to flip the electron at ^ = r, \il/^it > r)) = e-'^k^'-^^e-'^k^\ Ty^). ^1 = 1- |<<A^I<Apr 
possesses a simple geometrical interpretation: the squared distance between the two conjugate pseudo-spin states on the Bloch 
sphere, which quantifies the entanglement between the electron spin and the pseudo-spin. 

A couple of justified simplifications can provide an understanding of the eff'ects of various mechanisms on the spin deco- 
herence. First, the energy cost by the diagonal nuclear coupling {Dk) can be neglected as it is by three orders of magnitude 
smaller than that by hyperfine interaction {Ek). Second, for near-neighbor pair-flips, the intrinsic nuclear interaction is much 
stronger than the hyperfine mediated one for the field strength under consideration. Third, for non-local pair-flips, the intrinsic 
interaction is negligible due to its finite-range characteristic. Thus we can separate the flip-pairs into two subsets, Ka, which 
contains 0(N^) non-local flip-pairs driven by the efl'ective pseudo-magnetic field ^ (±2Ak, 0, ±Ek), and g Kb, which contains 
0(N) near-neighbor flip-pairs driven by ^ (2Bk, 0, ±Ek). The conjugate pseudo-spins will precess along opposite directions 
in the non-local subset Ka, and symmetrically with respect to the y-z plane in the near-neighbor subset Kb. The decoherence can 
be readily grouped by the two difl'erent mechanisms as 

keKe keKA 

where hk = |h*| and sinc(x) = sm(x)/x. 

In III-V QDs, because of the large number of nuclear spins, the lost of electron spin coherence is much faster than the build- 
up of higher-order nuclear spin correlations. The decoherence is therefore well described by the pair-correlation approximation 



as given in Eq. ( |B14| ). In other systems such as Si or diamond NV centers with a dilute nuclear spin bath, corrections from 
higher-order nuclear spin correlations will become important 1129112201422211 . 

Lattice distortion can result in local electric field gradients, inducing the quadrupole interaction for nuclear spins with moment 
greater than 1/2. Recent experimental works have indeed demonstrated signatures of quadrupolar interactions for nuclear spins in 
InAs self-assembled dot I367[ |368 l. The quadrupolar interaction can be well incorporated in the theory described in this appendix 



as contributions to energy cost for nuclear pair-flips [i.e. the Dk term in Eq. ( B12[ )] when reliable parameter is extracted from 



experiments. We expect that quadrupolar interaction does not afl'ect electron spin free induction decay and Hahn echo decay 
where the Dk term is unimportant, but may afl'ect Carr-Purcell echoes and spin echoes by complex pulse sequences when the Dk 
term plays a non-negligible role |[T27l[T3Ql[T34ll22T1l . 

Appendix C: Quantum measurement in Shor's algorithm 

We show by explicitly examining Shor's algorithm the crucial role of quantum non-demolition measurement for quantum 
computation to be scalable. 

1. Order finding for Shor's algorithm 

For the reader's convenience, we give a brief review of the algorithm for finding the order of a number, the core subroutine for 
Shor's algorithm. For a comprehensive description of the algorithm, the reader is referred to Ref. 11721 . The order r of a number 
X with respect to a number N is defined by the relation x^ = 1 mod N. For 2^"^ < N < 2^ and 1 < x < N, the task is to find r 
with resources at most polynomial to L. 

The observation is that f(n) = x^ mod N is periodic with period r, i.e., f(n -h r) = f(n), indicating a quantum Fourier 
transformation (QFT) may be used to find the spectrum of this function and thus to find r. The algorithm is outlined as follow. 

1. Two registers with t and L qubits, respectively, are zeroed initially, and thus the initial state is, 

|00 • • • 0)|00 • • • 0). 

2. The QFT is applied to the first register to obtain the state, 

2^-1 



1 



E 

where n = bib2 " - bt is 3. binary number with /?/ = or 1. 



n}\00 • • • 0), 
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3. With the advantage of quantum paralleHsm, one evaluation of the function f(n) is added to the second register to reach the 
state 



By f(n) = fin + r), the state can be rewritten as, 



1 ^'"^ 

— y |^)|x"modA^). 



r-l [(2^-^)/r] 



k^O /=0 

where [n/m] denotes the greatest integer not greater than n/m. Such an expression suggests a solution of the spectrum by 
QFT. 

4. After an inverse QFT applied to the first register, the state becomes. 

If 2^ happens to be an integer multiple of r, the terminating state is just, 

-ZZ^"™|2'./r>|/>, 

only for j = 2\slr) will the amplitude be nonzero. Generally, r would not divide 2\ but if 2^ is much larger than r, the spectrum 
after the QFT will be distribution composed of peaks around 2\s/r) for (0 < s < r). The larger is the first register, the sharper 
will the peaks be. The probability of the state of the first register being away from l^sjr by a distance 2^ is calculated to be less 
than 1 / (iP^^ - 4^, so with probability greater than 1-1/ (l^^^ - 4^, the fraction sjr can be determined up to the first t - p bits 
by a measurement of the first t- p bits of the first register. If ^ - is chosen to be 2L + 1 = N, r can be determined from the first 
bits of the binary fraction number s/r, i.e., [2^ s/r], by continued fraction. 

2. Issues with the measurement 

The key feature of Shor's algorithm is that, though the terminating state is a superposition of many computational basis states 
Yjx M (where \x) = \b1b2" - bt) with bi = or 1), it is not necessary to know all the amplitudes to solve the problem. Actually, 
an ideal measurement on the computational basis will project the superposition state into an arbitrary basis state which has 
a nonzero amplitude, and with high probability, the fraction sir can be determined up to 2L + 1 bits. If the measurement is 
performed in a single shot, the register may be read out bit by bit, and the superposition state will collapse in a cascade manner, 
so the resources required by the whole readout step is 0(L) and thus the measurement is scalable. The cascading readout can 
be illustrated by the example of reading the state |000) + |010) + |110) + |111), in which the state collapse may follow the steps 
shown in Fig. [44] Only single-bit measurements are required to have a N-hii superposition state collapsed into a basis state 
and only bits of classical memory are needed to record the measurement result. So a single-shot measurement on a single 
quantum register needs less resources than polynomial to the problem size in Shor's algorithm. 

3. Ensemble measurement 

Now we will show, that, in Shor's algorithm, an ensemble measurement requires resources increasing exponentially with the 
size of the problem. The terminating state of Shor's algorithm can be written as. 



^ r-l r-l 



where sir denotes the first N binary bits of sir, and R{s, k) denotes the state of the last p bits of the first register which are not 
accurate in describing sir. Only the first N bits of the first register need to be measured. The detected probability is uniformly 

distributed among the r states |^/r^, which are spaced almost equally by the distance 27 
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a. Correlated measurement 



Suppose the first N qubits of the first register are measured with coincidence counting. Each basis state \b1b2 • • • bf^) could 
lead to a click in a corresponding channel. To accumulate confidence in an ensemble measurement, a channel should get at least 
two clicks. In each single-shot trial of the measurement ensemble, the state would collapse into diff'erent states. Each output 
has to be recorded and stored before one of them is confirmed. So the number of recording channels and the size of classical 
memories used for data storage scale as 2^, exponentially increasing with the problem size in Shor's algorithm. This requires 
exponentially increasing physical resources such as spectral resolution or spectrometer bandwidth in spectroscopy. 

If the number of available channels and classical memory registers are limited by C. We may randomly or uniformly choose 
C x's from the 2^ possible numbers. What is the chance that we can receive a signal at any one of these C channels? As the 
probability of having a signal in an arbitrary channel is ~ 2"^, the probability of having signal in any of the C channels is 
Pc = 2~^C. So, for any finite probability Pc, the number of channels needs to be C ~ 2^Pc, increasing exponentially with the 
problem size. 



b. Uncorrelated measurement 



In an uncorrelated measurement of an A^-qubit register, the ensemble is first divided into 4^(N/g) portions (g is a fixed small 
integer), and every 4^ portions can be used for 4^ independent g-qubit measurements to obtain the density matrix of these g 
qubits. Once the reduced density matrices for all the g-qubit subsystems have been obtained. One could apply some classical 
algorithm to re-construct the density matrix of the N qubits. The re-construction, of course, could not be certain, but the 
uncertainty nonetheless does not exclude the possibility that one might search for a correct result from all possible states which 
give the reduced matrix elements. The number of measurements, the number of recording channels, and the size of storage, all 
these resources scale only linearly with the problem size. 

However, there is a fundamental problem underlying the uncorrelated measurement strategy: to yield an A/^-qubit output, 
essentially N bits of information have been generated. To generate the same amount of information, the g-qubit reduced density 
matrices have to be measured with - ln(N/g) bit accuracy. By the rule of thumb in experimental physics, measuring any 
physical quantity with A/^-qubit precision would require resources scaling as 2^. Only in exceptional cases, may the register state 
be derived from the knowledge of the reduced density matrices. For instance, if all one-qubit density matrices are pure states, 
the register state is obviously the outer product of all of them. But in general, it is difficult to determine the register state from 
the one-qubit density matrices. For example, if all the N qubits are maximally entangled with other qubits, such as in the states, 

10000) 10011) -h 10110) -h 11001) -h 11100) -h 11111), 

and, 

10000) -h 10101) -h 11010) -h 11111), 

the uncorrelated measurement would turn out to be N maximally mixed density matrices in both cases, from which little infor- 
mation can be obtained about the register state. To determine the order r, one has to search from all possibilities, of which the 
number is ~ 2^. 

Below we will show that, for the first N qubits in Shor's algorithm, there are ~ 2^ possible terminating states which would 
produce almost the same one-qubit reduced density matrix. 

o^looo) Q > |ooo) 

^|ooo)+|oioX;^. , , , V 
Q 77I / I / 7-^ 010) — iL^ioio) 



|000) + |010) + |110) + |111) 




1 \l 1 

|llO) + |lll)— ^|110) + |111) 



i^liu) 



FIG. 44. The collapse (quantum jump) of a multi-qubit state under a measurement of the qubits in sequence. The number associated with each 
arrow indicates the output of the measurement. 
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Lemma 1. For a state X/Lq I^^) (k is an odd number greater than 1), the reduced density matrix of any qubit has no 
off-diagonal term in the computational basis. 

Proof If there is an ofF-diagonal term (0|p|l) for the jth qubit, there have to be at least two states \x) and \x') in the 
superposition, which are different only at the jth bit. So we have that |x - = 2""^ is divided by k, which is impossible since k 
is an odd number □ 

Corollary 1. For the state ^J^o \l2^k} (k is an odd number greater than 1 ), the reduced density matrix of any qubit has 

no off-diagonal term in the computational basis, and the reduced density matrices of the last p qubits are all |0)(0|. 

Corollary 2. For the state 2/^o 1^^^^' reduced density matrices of the first n - p qubit are all (|0) + |1)) ^ ((0| + <1|) , 

and the reduced density matrices of the last p qubits are all |0)(0|. 

Lemma 2. For the state ^[lo ^^^^^ 1^^^' if the odd number k < 2^^^^ for a specific number a e (0, 1), the states \x) with xj = 
or 1 have the same probability up to cN-bit precision to occur in the superposition, where c is a constant less than 1. 

Proof The integer numbers x = {xiX2 • • - x^} with xj = or 1 form alternatively 2^ segments with length 2^"A For an 
arbitrarily chosen number ft e (a,l), the segment length 2^"^ is either less or greater than 2^^^^. If 2^"^ > 2^^^^, as the number 
of the multiples of ^ in a segment is greater than 2^"^^^/^. The numbers of multiple k in two neighboring segments differ by at 
most 1 . So the occurring probability of = is different from that of = 1 by at most 1/2(^-^^2. If 2^"^ < 2^^/^ we observe 
that the first x of every k segments is k2^~^, a multiple of k. As k is an odd number, each k segments starting with xj = will be 
followed by k segments starting with Xj=l, and vice versa, until the end of all segments. So, in every 2k segments, the numbers 
of Ik's with Xj = or 1 are exactly the same. As the difference of the occurring numbers in 2 neighboring segments is at most 
1, the difference in k segments is at most (k + l)/2. So, the difference of the occurring probability in all the segments is at most 
(k + 1)2^-^/2^+1 ^ ii2^(i-a/2-m)^ Let c = min(l - a/2 -yS/2,/5/2 - a/2), we have the probability of occurring of xj = Ois the 
same as that of occurring of xj = 1, accurate up to cN bits □ 

From the Lemma 1 and 2, we have directly the following theorem. 

Theorem 1. For all states zEf'^^^^"' \lk} [k is an odd number, and k < 2^"^^ for a specific number a e (0, 1)7, all one -qubit 
reduced density matrices obtained by uncorrelated measurements are the same up to cN significant bits. 

Corollary 3. For all states t^^^-i^/^^/^^)] 

is an odd number, andk < 2*^^^^ for a specific number a e (0, 1)7, all one-qubit 
reduced density matrices obtained by uncorrelated measurements are the same with cN significant bits. 

The theorem and corollary above are consistent with a theorem recently proved by Popescu et al 13691 : Almost all A/^-qubit 
states would give almost the same g-qubit reduced density matrix, as long as N is large and N g. 

The terminating states of the register to be measured in Shor's algorithm have the form of the superposition states in the 
theorem above, with at most one-bit deviation. So, unless at least the density of matrix of one qubit is measured with 0(N) 
effective bits, there are ~ 2^ possible superpositions corresponding to the same set of one-qubit reduced density matrices. On 
the one hand, searching the correct one from all those possibilities needs resources ~ 2^ in all known classical or quantum 
algorithms. On the other hand, determining the density matrix of a qubit with 0(N)-bit accuracy also requires resources ~ 2^. 
So uncorrelated single-qubit ensemble measurement is provably unscalable for Shor's algorithm. Though there is no proof in 
general cases, it would be rather surprising that some ensemble measurement scheme is scalable for Shor's algorithm. 



4. Single-object measurement with error 



In general, scalable quantum computation needs to be performed on a single quantum object (rather than an ensemble) with 
single-shot measurement. In reality, detectors used in the readout procedure have unavoidable inefficiency or errors. Thus, the 
Kraus operators l>172il for a POVM (positive operator- valued measure) of a certain qubit can be written as 

Ao = Vl - J|0)<0| Vr^|l)<l|, (Cla) 
Al = VJ|0)(0|+ V^|1)<1|, (Clb) 

with detector efficiency e and dark count rate d. For a state \i//}, the probability and the resultant state for the output and 1 are, 
respectively, 

Po = <(AKAo|<A)andAo|<A), (C2a) 
Pi = {i//\AlAi\i//}sindAi\i//}, (C2b) 
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Suppose that the detecting error rate at each qubit is greater than a finite number the probabiHty of reading out a A/^-qubit 
register is less than (1 - e)^ , exponentially small with increasing the size of the register. A single-shot measurement then is 
insufficient for scalable quantum computation. Thus, if there exist detector errors as always, the measurement has to be repeated 
for enough times to obtain suflficient confidence in a readout result. Furthermore, if the measurement is destructive, the quantum 
computation has to be rewound from the very beginning of the algorithm, making the measurement equivalent to an ensemble 
measurement. 

In an uncorrected measurement, diflferent qubits of the register are measured and recorded independently, and the error rate 
at each bit is finite, so the density matrices of each qubit of the register can be measured by repeating the quantum computation 
for a number of times proportional to the register size. But the problem here is the same as discussed for the uncorrected 
measurement in the previous section. 

In a correlated measurement, the probability of correctly reading out the projected state is exponentially small [~ (1 - sf\ 
and yet the probability of the terminating state collapsing into the same basis state is also exponentially small (< 1 /r). So before 
a readout result is repeated once for accumulating sufl&cient confidence, the quantum computation has to be repeated a number 
of times which increases exponentially with the problem size. 



5. Quantum non-demolition measurement 



In a QND measurement, the state will remain unchanged after the projection into the measurement basis (which is also the 
computational basis). So the readout can be repeated many times in every qubit to accumulate confidence of the readout result, 
without rewinding the whole algorithm from the beginning. 

Now we calculate the resources required in reading out the state of a A/^-qubit register. If the error rate in reading out each qubit 
by an M-shot QND measurement is sm, the probability of successfully reading out the register is (1 - suf - To have a finite 
success probability, we require ^-a^ = (1 - em)^ > s where 5 is a finite number smaller than 1. When sm is small, ~ e~^^^, so 
we require the error rate of a M-shot QND measurement sm < -(^ /N) In s. 

To obtain the error rate of a M-shot measurement, we define its POVM's. The Kraus operators for the POVM of a M-shot 
measurement giving m photon counts can be derived as, 

vl/2 



= (i-^/^-"^^/'^-/'|0)<oi 



When dM,m < eM,m, it is more probable that the qubit is in the state |1>, and vice versa. As 



dM,m/^M,m - | "j 

I - e 



e{l - d) 



(C4) 



monotonically decreases with m, we can define a mo so that all AM,m<mo indicators of |0> and all AM,m>mo indicators of |1>, 
with the mo given by dM,mJeM,mQ = 1 or, 

Mln^ 

^0 = , , ^ ^M. (C5) 

In f-^ -h In 5 

l-e d 

As dM,m and eM,m as functions of m have peaks at dM and eM, respectively, we have, 

d < a < e. (C6) 

Now the POVM can be calculated from, 

> ^ 2 Al^^AM,m ^ (1 - ^m)|0)<0| + (1 - ^m)|1><1|, 

m<mo 

= A^^^Au^m = dM\0){0\ + eM\m\, 

m>mo 
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where, 



m>mo m>mo 

Ml 



(M-mo)\(mo 
M!(l - 

(M-mo)!(mo-l)! 

(1 - 



— — r r^-^o(l-^r«"^^/^ 
- I)- Ji-d 



aM 
ln(\ - a) 



(1 - Qf)l-^Qf^ 



and similarly, 



\-eM< 



\ aM(\-ef 



(1 - 



(1 - 



M 



(C7) 



(C8) 



The error rate is defined as = max(JM? 1 - ^m)? so the number of required QND measurements per qubit M ~ InA/^, and 
the total number of measurement is proportional to N In So a QND is scalable. In experiment, e and J cannot be determined 
exactly, but fortunately a need not be determined exactly and whenever it is between d and e, the results for du and eu are 
unchanged in the equations above. 

Note that with the aid of entanglement gates and a supply of fresh qubits, a destructive measurement can be converted into 
a QND one. The idea is based on the transformation of a qubit and M zeroed auxiliary qubits into an entangled state by M 
entanglement gates: 

M-l M-1 M-1 

(a|0) +j8|l)) (g) 10) ^ a|0) (g) |0) +yS|l) (g) 

m=0 m=0 m=0 

The M auxiliary qubits are to be read out. As all these qubits are entangled, once one qubit is collapsed into a basis state |0) or 
|1>, all the qubits will be collapsed into the same state. So even a destructive measurement with detecting error can be used to 
read out the qubit. 

How about a QND measurement with back-action noise, i.e., reaction to measurement that disturbs the state after a mea- 
surement cycle? With the idea above for converting destructive measurement into a QND one, we can employ the concept of 
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FIG. 45. Some elementary quantum gates and their matrix representation, including (a) Hadamard (//), (b) rotation of 90 degrees about the 
jc-axis (X2), (c) NOT (A^), (d) SWAP (W), (e) Controlled NOT (C), (f) single-bit phase-shift (S^X and (g) Controlled phase-shift (P^). 
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FIG. 46. Realization of several control gates, including (a) CNOT, (b) remote CNOT, and (c) Toffoli (C2). 



fault-tolerant measurement to deal with this problem. In the so-called fault-tolerant measurement, a qubit is first encoded into a 
stabilizer code, after a single measurement, any back-action noise will be diagnosed and corrected using the error syndrome since 
this noise acts only on a single qubit (by assumption). After the error correction, another measurement would be performed, and 
so on. The fault- tolerant measurement thus allows an imperfect QND measurement to read out the result of a quantum algorithm 
with polynomial resources. 



Appendix D: Elementary quantum gates 



Here we define a series of elementary quantum gates used in the quantum circuits presented in Sec. VIII A The quantum gates 
of interest are defined in Figs. [45] and [46] The Hadamard gate can be realized by the spin rotation operations up to a trivial global 
phase, as the transformation operator for the Hadamard gate can be expressed in terms of spin rotations sls H = e^^/^e^^^z^i^Sy/i^ 
The CNOT gate can be realized by the controlled 7r-phase gate together with two Hadamard gates. The CNOT gate on two 
remote bits can be realized by a local CNOT gate plus some SWAP gates. The doubly controlled NOT (Toff'oli) gate can be 
realized by six CNOT gates plus some single-qubit gates. The controlled 7r-phase gate and the CNOT gate can be realized by 
two VSWAP gates plus some single- spin rotations. As any controlled phase gate can be realized by two CNOT gates plus some 



single-qubit gates, the circuits in Figs. 38 and 39 can be alternatively realized with VSWAP gates as well. 
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